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I. MACROOTLOCUS SYSTEM 



A. INTRODUCTION 

MacRootLocus system is a program written for the Apple Macintosh 
computer. It was designed to allow a user to analyze and design linear feedback 
control systems. It is written in the computer language Turbo Pascal which is the 
native language of the Apple Macintosh and designed with the same 
user-friendliness and standard interface philosophy the Macintosh was designed for. 

The Macintosh gets most of its input from the user through the 'mouse' which 
is a small cigarette pack— sized control that is moved across the table much like a 
pencil across paper to move the cursor or pointer on the computer screen. Rather 
than typing in commands like the IBM, the Macintosh lets you select the function 
you w'ant performed with the mouse. Since the commands are not typed in, the 
commands do not have to be remembered as with other computers. It is the 
user-friendliness that sets MacRootLocus apart from other system analysis 
programs. Prior computer experience is not required to use MacRootLocus. A few 
minutes to learn how to use the mouse and pull down menus is all that is needed. 
All dialogs are easy to use and there are on-line Help menus. So the first— time user 
can get desired results without using trial and error. 

Apple Macintosh models supported by MacRootLocus are the Macintosh SE 
and the Macintosh II system. 

MacRootLocus was tested with several examples and is now available to any 
user on the Naval Postgraduate School Controls Laboratory computers under the 
icon of MacRootLocus system. 
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B. MATHEMATICAL CONCEPT 

The root locus method is a graphical technique for determining the roots of the 
closed— loop characteristic equation of a system as a function of the static gain. 
Consider the general feedback control system, as shown in Figure 1.1. 




Figure 1.1 Closed— loop for Control System 
In order to find the roots of the characteristic equation, it is required that 

1 + G(s)II(s) = 0. (1.1) 

Of course, Eq. (1.1) may be rewritten as 

G(s)H(s) = -l. (1.2) 

Since s is a complex variable, Eq. (1.2) may be rewritten in parametric form 

|G(s)H(s)| = l (1.3) 
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L G(s)H(s) = (2k-l)7T 



(1.4) 



where k is an integer. 

For a specific value of s to be a root of the characteristic equation, it must 
satisfy both Eqs. (1.3) and (1.4). Since the roots are those values of s that satisfy 
both equations, then the root points are points where these curves intersect. 

The original development of the root locus method was concerned with the 
determination of the locus of roots of the characteristic equation as the system gain, 
K, is varied from zero to infinity. It appears that the root locus method is a single 
parameter method; fortunately it can be readily extended to the investigation of two 
or more parameters. 

The characteristic equation of a dynamic system may be written as 

a S + a ^ + ....~l"aS~l"a =0. (I-^) 

n n-l 10 

Clearly, the effect of the coefficient ai may be ascertained from the root locus 
equation 



a S 

1 + 1 . ( 1 . 6 ) 

a S" + a S"-i + . . . . + a S2 + a 

n n-l 2 0 



The parameter of interest, A (in MacRootLous), can be isolated as 



a S" + a + + (a - A)S"''' + AS"'"’'! + . . 

n n-l n-q 

+ aS + aS = 0. (1.7) 

10 
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For example, a third— order equation of interest might be 



+ (3 + A)S2 + 3S + 6 = 0. (1.8) 

In order to ascertain the effect of the parameter A, we isolate the parameter and 
rewrite the equation in root locus form as shown in the following steps: 

S3 + 3S2 + AS2 + 3S + 6 = 0 (1.9) 

A q2 

1 + — = 0 . ( 1 . 10 ) 

S3 + 3s2 + 3S + 6 

Then, to determine the effect of two parameters, we must repeat the root locus 
approach twice. Thus for a characteristic equation with two variable parameters, A 
and B (in Mac RootLocus), 

a S" + a S"-i + + (a - A)S"-‘^ + + . . 

n n-1 n-q 

+ (a -B)S"-' + BS"-'-i + . . + a S + a = 0. (1.11) 

n-r 1 0 

The two variable parameters have been isolated and the effect of A will be 
determined, followed by the determination of the effect of B. For example, for a 
certain third-order characteristic equation with A and B as parameters, we obtain 

S3 + S2 + BS + A = 0. (1.12) 
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In this particular case, the parameters appear as the coefficients of the characteristic 
equation. The effect of varying B from zero to infinity is determined from the root 
locus equation 



1 + 



BS 

+S2 + A 



= 0 . 



(1.13) 



One notes that the denominator of Eq (1.13) is the characteristic equation of the 
system with B = 0. Therefore, one first evaluates the effect of varying A from zero 
to infinity by utilizing the equation 



S3 + S2 + A = 0 



(1.14) 



rewritten as 



1+ = 0 (1.15) 

S2(S + 1) 

where B has been set equal to zero in Eq. (1.12). Then, upon evaluating the effect 
of A, a value of A is selected and used with Eq. (1.13) to evaluate the effect of B. 
This two-step method of evaluating the effect of A and then B may be carried out 
as a two— parameter root locus procedure. First, we obtain a locus of roots as A 
varies, and we select a suitable value of A; the results are satisfactory root locations. 
Then we obtain the root locus for B by noting that the poles of Eq. (1.13) are the 
roots evaluated by the root locus of Eq. (1.15). 

In Mac RootLocus, this design approach is used to calculate the roots for the 

plot. 
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II. OPERATING THE MACROOTLOCUS 



A. BASIC MACINTOSH USE 

This document must serve as a user's manual as well as a technical explanation 
of the program and its capabilities. For this reason, the following brief explanation 
of Macintosh use is included. It is by no means a substitute for the Apple 
Macintosh Users' Manual but it will contain enough information for the beginner to 
be able to use Mac Root Locus. 

1. Basic Operation. 

The Macintosh has a finder, a special application you use to organize and 
manage your document and to start other applications. You use the finder every 
time you start your Macintosh, or whenever you move from one application to 
another. 

The Macintosh screen looks like a light gray desktop, rather than a 
textual list of commands and responses. The desktop simulates the working 
environment. It is initially clean, displaying small graphic images, called icons, 
with short titles directly under them for each disk presently being used and a trash 
can in the lower right corner. An icon is an image representation of an application 
document or a control to a usable function. They offer quick recognition as to the 
type of item they describe and are easier to identify than lists or directories of file 
names with extensions. 

The main interface between the user and Macintosh is the mouse. The 
Macintosh responds instantly to every movement you make with the mouse. You 
can start applications and get documents, work on them, and put them away again 
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just by moving the mouse and pressing the mouse button. 

There are three techniques for the mouse: pointing, clicking and dragging. 
Moving the mouse moves the cursor or pointer on the screen in the same direction. 
Positioning the pointer on an item is called pointing to it. The mouse is used to 
select various items or icons on the screen. You select any item or icon to let 
Macintosh know this is what you want to work on next. You select icons by using a 
technique called clicking. As you click the icon, it becomes highlighted. This 
highlighting shows that you select it. This mouse can also be used to drag across 
something. This means the button is pressed and held while the mouse is moved. 
This action is called dragging. When dragging the mouse across the screen, a 
rectangle is outlined. When the button is released, everything within the rectangle 
is now highlighted and selected. Dragging also refers to moving items on the screen. 
This is done by pointing to an item, pressing and holding the button and moving 
the mouse. This will also move an outline of the icon selected and when the button 
is released, the icon is moved to the new location. 

Whenever you work with Macintosh, you tell it two things: what you 
want to work on and what you want to do. First, you tell the Macintosh what you 
want to work on by selecting it as you have been doing with icons on the desktop. 
Then you tell the Macintosh what you want to do with the selection. You usually 
do this by choosing a command from a menu. 

Along the top of the screen, in the menu bar, are titles of the menus. 
Pressing the mouse button while you are pointing to a menu title causes the title to 
be highlighted and a menu to appear, much like a window blind being pulled down. 
The menu contains commands you can carry out on what you have selected. 
Commands that you cannot use right now appear dimmed in the menu. Wdien you 
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release the mouse button, the menu disappears. 



To choose a command from a menu, you use the same dragging technique 
you used to move icons. As you drag through a menu, each usable command is 
highlighted in turn. If you change your mind about choosing a command, move the 
pointer off the menu and release the mouse button. Nothing is chosen unless you 
release the mouse button while one for the commands is highlighted. You'll follow 
this same pattern whenever you work with the Macintosh; 'select' some information, 
the 'choose' an action for it. For example in MacRootLocus, if you double click the 
mouse on MacRootLocus icon, the icon will be selected and MacRootLocus 
application will start running. You can see the menu bar along the top of the screen 
with the greeting message in the center of the screen. Now, you can select and 
choose the items to use MacRootLocus for your design. 

2. Manipulating Window 

The window is the area that displays information on the desk top. You 
view application documents and folders. Folders are a way of storing and organizing 
things, much like in a file drawer. Double clicking on a folder will open a window 
that displays its contents. Folders can be located in disks or within other folders. 
Usually 8 to 12 windows are the maximum that can be open at any one time. When 
several windows are displayed at once, they will usually overlap. If the window you 
want to view is partially covered by another window, pointing anywhere in your 
window and clicking will select it and make it the active window. The active 
window is always in front of all the others. It is the place you want the next action 
to happen, such as move or select icons or open other folders. 

The active window can be identified by its highlighted title bar with 
narrow horizontal lines on either side of the title. The active window also usually 
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has a close box (to close the window) at the top— left corner; at the top— right corner 
is a zoom box that expands the window until it nearly covers the screen. On the 
bottom-right corner is the size box you use to change the size of a window. As you 
click or drag these boxes, you can get the function that you want. 

3. Handling Input 

Most information in the form of data or text is entered through the 
keyboard. The keyboard includes character keys, numerical keys, direction keys and 
other special keys. The return key tells the computer that the data just typed in 
should be accepted now. In MacRootLocus, the apple key in combination with 
another key is often a shortcut to choose a command from a menu. The tap key is 
used to move between data input points in the dialog for entering data. The 
direction key is important in MacRootLocus. It will be used for some word 
processing applications in a small box of the plot window. This will be further 
explained in the next section. 

When the Macintosh requires information from you, it will display a 
dialog box like the one shown in Fig 2.1. It will tell you exactly what data is needed 
and shows you where to enter it. Data insertion points are small boxes in the dialog 
box that let you type in numbers or text. There are usually several such insertion 
points in each dialog box. The tab key lets you move from one point to another to 
enter data. Usually there will already be data in an insertion point box. This is the 
default data. You can change it if you want but you do not have to. When you 
enter new values in an insertion point box, they will become the new default values. 
Data insertion boxes can also be selected by pointing and clicking with the mouse. 
Clicking the mouse between two characters in the box will let you insert characters 
between them. Double clicking a box will highlight the entire number or an entire 
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Figure 2.1 Sample Dialog 

word if text is entered. You can also select all or portions of numbers or text by 
dragging the mouse across the text you want to select. When all or part of a text is 
selected, it will be highlighted. It can be removed by using the delete or backspace 
key, or it can be replaced by typing in whatever you want. After all data in the 
insertion box is correct, you can enter the data by hitting the return key or by 
clicking the OK button. If you click on Cancel, any changes to the data in the 
insertion boxes will not be saved and the operation which called the dialog box will 
be canceled. If you type in data that does not apply to the insertion box, such as 
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typing in letters in the AMin Gain insertion box of one— parameter plot data dialog, 
the Macintosh beeps. It informs you that you inserted the wrong input and you 
should correct the input that you just typed. 

4. Printing Out 

In order to print out your work, there are a few ways available in 
MacRootLocus. One way is to select the print command in the File menu. This 
allows the user to get a hard copy of any plot displayed by MacRootLocus. 

Another way is to do a screen dump which will print the contents of the 
active window immediately. This is done by holding the command key and then 
typing the number '4'. This is a fast way to get print out of a plot in 
MacRootLocus. You can also create a MacPaint document by pressing the 
command and shift key and typing the number '3'. You can take up to 10 of these 
'snapshots' and can then alter them with MacPaint for transferring to a word 
processor for lab writeup. 

B. MENUS AND DIALOGS 

There are five menus for MacRootLocus: Apple, File, Edit, Plot and Help. 
Actually, the File, Plot and Help menus are used to get the plot. These have to be 
used in order which will be explained later. 

The Apple and Edit menu are not directly used by MacRootLocus, but they are 
used in order to follow the standard Macintosh programming philosophy for 
user-friendliness. This allows the experienced user to easily adapt to a new 
program since many of the operations are already familiar. Also, this is to allow for 
easy interaction with various other programs and desk accessories [Ref. 1]. 

As mentioned earlier, MacRootLocus presents commands in menus you pull 
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down from the menu bar. As soon as you choose the command you want, the dialog 
box appears for each command. This allows you to insert input data. Now more 
details for menus and dialogs in MacRootLocus will be explained. 

1. Apple menu 

The Apple menu is identified by a small apple in the top left corner. It is 
used primarily for desk accessories. It is also used by most programs as a way of 
offering program information. This is usually the first item of the Apple menu. 

There are several accessories for some applications. For example, the 
Mac's clipboard and scrapbook are used to interact with the word processing 
program. It is highly recommended that the user read the Macintosh Users Guide 
as it explains the use of the Mac's Accessories which will be of use but will not be 
discussed here. 

2. File Menu 

MacRootLocus starts with the file menu first. It offers EQ Parameter, 
Get Coeff, Print Screen, Print Window and Quit commands. These selections also 
have a keyboard shortcut by holding down the command key, which has a clover 
leaf symbol on it, and hitting the first letter of the menu item at the same time. 

0 . EQ Parameter 

The EQ Parameter stands for Equation Parameter. As mentioned 
the previous chapter, the MacRootLocus program calculates the roots of the 
characteristic equation for plotting points. It is neccessary to get the degree of 
polynomials and some equation parameters shown in Figure 2.3 since the Laguerre 
algorithm [Ref. 2] is used in order to calculate the roots. 

The degree of the polynomial should be between 1 and 10 in 
MacRootLocus. Then there are the default values for the other parameters. These 
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values avoid the convergence error for almost all polynomials. But if convergence 
error messages appear on the screen, you can change these parameter values. These 
parameters must satisfy the following conditions: 

(1) Initial guess > 0 

(2) Maximum Iteration > 100 

(3) Tolerance > 0 



Characteristic Equation Parameter 
Degree of the polynomial 



InitGuess 

Maniter 

Tolerance 



J 

100 

IE-6 



( ^ 

OK 

V > 

/ — \ 

Cancel 

V 



Figure 2.2 Equation Parameter Dialog 
b. Get Coeff 

After getting the degree of the polynomial, you should insert the 
coefficients of the polynomial. The 'Get Coeff stands for Get Coefficients. The 
dialog for the coefficients is shown in Figure 2.3. This dialog is varied depending on 
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the degree of the polynomial as shown in Figures 2.3 and 2.4. 

The algebraic expression for the coefficients of the characteristic 
equation of the system may have up to two undetermined parameters (A and B). In 
the case of the one— parameter root locus method, you use only one undetermined 
parameter (A). The routine uses standard algebraic, or infix, notation with 
parenthesis allowed. Operators can include +, — , *, /, and * (exponentiation). The 
unary minus sign is allowed. For example, the characteristic equation is 
S3 + (10 + A)S2 + (10 * A + 5000 * B)S + 5000 * A = 0 
Figure 2.3 shows how you insert these coefficients. 



Characteristic Equation Coefficient Data 



* > 

Cancel 

- - ^ 



S**3 



10+fl 



S**l 

10*fl+5000*B 






5000*fl| 



Figure 2.3 Characteristic Equation Coefficient Data Dialog 
Box for Third Degree Polynomial 

If you choose the Get Coeff command without the degree of the polynomial. 
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the message shown in Figure 2.5 appears on the screen. This message tells you that 
degree of the polynomial has not yet been entered. 



Characteristic Equation Coefficient Data 



f .... I ^ 

OK 



Cancel J 



$♦♦7 



S**6 






Figure 2.4 Characteristic Equation Coefficient Data. Dialog 
Box for Seven Degree of Polynomial 





There is no Initial Degree. 







Figure 2.5 Message (1) 



c. Print Screen and Print Window 

These comands allow the user to get a hard copy of any plot 
displayed by MacRootLocus. When you choose the 'Print Screen' command, the 
whole contents of the screen is printed. If you choose the 'Print Window' command, 
the plot that is to be printed should be on the active window of the display. 

Before attempting this method, ensure the printer is properly set up 
for friction feed. 

d. Quit 

The last item under the File menu is Quit. Selection of this will 
cause you to leave MacRootLocus and return to the desk top. 

3. Edit Menu 

Since these operations are not actually used in MacRootLocus, no further 
explanation of the Edit menu will be presented here. Any additional information 
regarding the Edit menu can be found in the Macintosh Users Manual. 

4. Plot Menu 

There are two commands. One Parameter and Two Parameter. These are 
called to display the dialog boxes of Figures 2.6 and 2.8 to insert the plot data. 

These commands are followed after selecting the EQ Parameter and Get 
Coeff commands. If not, the message shown in Figure 2.7 appears to tell you that 
the degree of polynomial and polynomial coefficients should be entered before you 
choose the plot menu. After reading the message, just click once and this box will 
disappear. 

o. One Parameter 

When you choose the 'One Parameter' command, the dialog box 
shown in Figure 2.6 for plotting data of the one— parameter root locus method 
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One Parameter Root Locus Plot Data 
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Figure 2.6 One Parameter Root Locus Plot Data Dialog 



appears. The plot default values are shown in Figure 2.6. They can be changed as 
desired. 

First, the user enters the minimum and the maximum gain values 
into the 'Min Gain' and the 'Max Gain' insertion box. Next, the user selects one of 
two types of interval. Linear and Logarithmic. When you click the radio button, 
the desired type of interval is chosen. For the 'Linear' interval, the gain step size is 
calculated by subtracting the minimum gain from the maximum gain entered in the 
dialog box, and then dividing by the number of points to plot. Using 'Logarithmic 
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interval can best be described as giving equally spaced intervals on a logarithmic 
scale. 

Most MacRootLocus programs calculate the gain intervals using the 
'Linear' interval. This emphasizes gain values that are closer to the max gain. This 
becomes more evident as the maximum gain to the minimum gain ratio increases. 
Using the 'Logarithmic' interval gives more emphasis to the lower gains so more 
continuous loci can be drawn. As a basic rule of thumb, if the maximum gain to 
minimum gain ratio is greater than 100, selecting 'Logarithmic' interval will give a 
more continuous plot. 



There i$ no Initial Degree or 

Characteristic Equation Coefficient. 



Figure 2.7 Message Box (2) 

Next, the scale for the axis will be chosen. For the 'Auto Scale' the 
system calculates the minmum and maximum value of each axis. When the 
'Manual Scale' is chosen, the user should insert the minimum and maximum values 
for each axis. 

The last item, 'Points to Plot' sets the plot resolution. The bigger 
the value you choose, the better resolution plot you get, but the calculation time 
will be longer. 
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Figure 2.8 Two Parameter Root Locus Plot Data Dialog 
b. Two Parameter 

The 'Two Parameter' conmiaiid calls the two-parameter plot dala 
dialog. It is shown in Figure 2.8. The items shown in Figure 2.8 are similar lo 
those shown in Figure 2.6, but several items are different. 

There exists one more undetermined parameter 'B' to be inserted. 
The 'How many loci' item lets you decide how many loci are to be drawn for each 
parameter. The number of loci will be from 1 up to 10 for each parameter. In 
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Figure 2.8, this value is 5. There is no auto— scale for axes. Only the manual scale 
is available. You focus on the interesting area for your design. The last item is the 
marking and justification in order to draw the selected 'A' and 'B' values on the 
plot. There are four radio buttons. Two buttons are chosen each time for each 
parameter, one for position, the other justification to draw. There are sixteen 
combinations available for this work as shown in Table 2.1. Figure 2.10 shows you 
the ninth case in Table 2.1. 



Table 2.1 The Combination for Marking and Justification 
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S : Start Point E : Endpoint R : Right Hand Side L : Left Hand Side 



5. Help menu 

MacRootLocus supports an on-line help menu so that the first— time user 
can get desired results without using trial and error. 

Help is the last item in the menu bar. There are the same item names in 
the menu bar. It makes it easy to look for the item for which the user wants 
information. The contents of each item are the subject of section 2.B. 
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6. Information Box 

Finally, MacRootLocus gives you a convenient way to identify your plot. 
It is a small box in which you can type the information you want to memorize. As 
soon as the plot has been completed, the box appears at the bottom of the plot 
window automatically. 

When the up direction key is pressed twice, the cursor comes out on the 
box. Then you can use the keyboard just like a typewriter. The capacity of the box 
is 160 to 200 letters. If you do not want that box, just click once. It will disappear. 
It is shown in Figure 2.10. 
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Figure 2.10 Sample Plot (1) 
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HI. DETAILED PROCEDURE MODULE DESCRIPTION 



This chapter will basically be programmer's notes covering the significiant 
procedures, functions and libraries to be used in MacRootLocus. MacRootLocus 
follows the Macintosh programming technique to make the program user friendly. 
It is a fairly simple Macintosh application that uses menus, windows, dialog boxes, 
and graphics. Most procedures and functions called in MacRootLocus were 
developed in Turbo Pascal verson 1.0 for Macintosh. 

MacRootLocus consists of one main program, one resource file and six unit.s. 
The main program integrates the resource file and units, then it shows you tlie 
menu bar and the greeting message to start the work. The details for these will now 
be explained. 



A. MAIN PROGRAM 

The Main Program is named MacRootLocus.Pas. This program consists of a 
main body and some precedures to handle the system. 

1. Main Body 

The structure of the main body uses the concept of event driven 
programming. It looks something like this: 



Initialize; 

repeat 

SystemTask; 

if GetNextEvent (everyEvent, theEvent) 
HandleEvent (theEvent) ; 
until Finished; 

Cleanup; 



{ set everything up 1 

{ keep doing the following } 

{ update desk accessories } 

then { if there's an event... ) 
{ ... then handle it } 

{ until user is done 1 
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The program is set up once with the user— defined routine 'Initialize.' It 
then enters a loop that continues until some condition (such as the user selecting 
Quit in a menu) causes it to set the boolean flag 'Finished' to true. Within that 
loop, it performs two major tasks. 

First, it calls 'System Task' (a Toolbox routine), which allows the Mac 
operating system to update any desk accessories that might be in use. Second, it 
calls 'GetNextEvent' (another Toolbox routine) to see if any events have occurred. 
If any have, the highest priority event is returned in the data structure 'the Event.' 
The program then passes the event to 'Handle Event,' which is a user— defined 
routine that handles all the different events that might occur. Such events include 
key— presses, selection of menu items, mouse clicks, and windows being opened, 
closed, or resized. When the program is ready to terminate, it calls the 
user-defined routine clean up. 

2. Handle Event Procedure 

When an event occurs, the operating system creates an event record and 
puts it in a queue, ready for you to handle. To see if there is one waiting, you call 
'GetNextEvent.' a boolean function that returns true if there is an event there for 
you. You give it a mask of the events you are interested in; you can use the 
predefined mask 'Every Event' to look at all events. This event is passed to 'Handle 
Event,' which takes care of it . 'HandleEvent' is just a case statement using the 
'what' field in 'the Event' to determine which of the procedures to call. 

3. DoMouseDown Procedure 

The routine 'DoMouseDown' determines which window' the mouse was in when 
the clicking took place and where exactly it happened. Like 'HandleEvent,' 
'DoMouseDown' is mostly a case statement. 
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4. DoUpdate Procedure 

The Macintosh keeps track of a lot of things for you. For one, it tells you 
when some portion of a window needs to be drawn, because of resizing or removing a 
covering window. This is known as an update event, and it requires special 
handling. To handle an update event, the routine 'DoUpdate' saves the current 
'grafport' into 'SavePort' and makes 'the Window' the current port so that you can 
write to it. 'BeginUpdate' limits all output to the section of 'theWindow' that 
needs updating. You then do whatever redrawing is needed. When you are done, 
'EndUpdate' lifts those limits, and 'SetPort(SavePort)' restores the old 'grafport'. 

5. Dokeypress Procedure 

This routine handles the 'Key down' and the 'Autokey' events. It is a 
check to see if a command— key combination was pressed; if so, it checks if the key is 
a menu command and takes appropriate action. 

6. Handle Menu 

The procedure 'Handle Menu' decodes the mouse position and figures out 
which menu and which item in that menu were selected. It uses a case statement to 
select the action for the appropriate menu; the menu value is the ID assigned when 
the menu is created. The commands in a menu are numbered from the top down, 
with the first command having a value of one. The action itself is usually a second 
case statement, based on the menu item. 

When an item in a menu is selected, the name of that menu (in the menu 
bar at the top of the screen) is highlighted, that is, inverted to white-on-black. 
When you are done processing the menu command, the menu bar is restored to 
normal by calling 'HiliteMenu(O),' which is at the bottom of 'HandleMenu.' 
Another procedure is called before you can leave 'HandleMenu,' 'UpdateMenu,' a 
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local procedure that tests to see if certain items are to be enabled or disabled. To 
enable and disable menu items, the standard Macintosh procedures 'Enableltem' 
and 'Disableltem' are called. 

7. Initialization 

The initialization procedure for the MacRootLocus program includes the 
following structure: Call 'Init' routines, set up menus, set up windows, do other 
graphic initialization and do program— specific initialization. 

There is an 'Init' routine for most of the major managers. The first, and 
most important, is 'InitGraf (thePort). That sets up 'QuickDraw' (which is used 
by just about everything else) and sets up a 'grafport' for the screen. Other 'Init' 
routines are: InitFonts, InitWindow, InitMenus, TEInit, and InitDialogs (NIL). 

Setting up menus involves four steps. First, it defines the menus 
themselves. If a resource file is used, just do a call to 'GetMenu' for each menu 
handle, or even a single call to 'GetNewMBar.' Otherwise, it has to build each 
menu using an initial call to 'NewMenu', followed by a call or calls to Append 
Menu. Second, if it is handling desk accessories, call 'AddResMenu'. Third, add all 
the menus to the menu bar by marking successive calls to 'InsertMenu.' 'Finally, 
call 'DrawMenuBar' to display the menu titles and make them active. 

As with menus, the window initialization takes several steps. If 
MacRootLocus needs a window at start up, create it using either'GetNewwindow' 
(reading in from resources) or Newwindow (building it in place). Having created 
the window, make it the current 'grafport' by calling 'SetPort', then make it the 
active window by calling 'Select Window.' 

The program— specific initialization should probably come here. All of 
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the default values for the dialogs are defined and the array vectors are initialized 
here. 

B. GLOBALVAR UNIT 

This unit declares all of the whole variables to be used in MacRootLocus 
except a few of local variables. This unit is called by the main program and all 
units. It defines constant, data type, and variables. 

C. MAKEROOT UNIT 

The unit MakeRoot is very important in MacRootLocus. This unit provides 
several procedures and functions to find the roots and the array vectors for plot. 
Since the characteristic equation is derived, the program must be able to parse the 
user's characteristic polynomial coefficient equations in order to understand the 
relations and be able to iteratively substitute in values for undetermined 
parameters: A (for one-parameter root locus method) or A and B (for 
two— parameter root locus method). 

This unit has two main procedures, the Get Root 1 procedure for the 
one— parameter root locus method and the Get Root 2 procedure for two— parameter 
root locus method. The simplified algorithm for the two— parameter root locus 
method is outlined in Figure 3.1. 

There are several procedures to perform this algorithm. These will be 
explained in the following section except for the Rootfinder unit. The InfixtoPolish 
and the ComputePolish are especially interesting. 

1. InfixtoPolish and ComputePolish Procedure 

The coefficients of the polynomial equation are written in algebraic, or 
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Given a system's characteristic polynomial : 



C.E = a S" + a S"'* + .... + aS + a = 0 

n n-1 1 0 

where an, an-i, etc. are algebraic expressions in A and B. Also A 
is to be stepped from the minimum A to the maximum A and B varied 
from the minimum B to the maximum B and the reverse is processed. 

Then 

SET A = Min A . SET B = Min B. 

SET A DeltaStep = abs((A Max - A Min) / Step) 

SET B DeltaStep = abs((B Max - B Min) / Step) 

FOR 1 = 1 to 2 do { case of A and B parameter] 

FOR j = 1 to Step do {Step = quantity of losi] 

IF Step A then A = A Min + A DeltaStep * (j - 1) 

ELSE B = B Min + B DeltaStep * (j - 1) 

WHILE point no <= (points - 1) do 

IF Step A then B = NextGain2 
ELSE A = NextGain2 
FOR term = Initial Degree downto 0 do 
CONVERT ai from Infix to Polish 
SUBSTITUTE values for A and B 
COMPUTE ai 
END {FOR} 

CALL RootFinder 

CALL Results {make plot array and numerical data) 
END{WHILE} 

CALL PlotRootLOcus2 
END{FOR) 

Step A = False {Next case) 

END{FOR} 



Figure 3.1 Two Parameter Root Locus Algorithm 



'infix' notation. The available operators include (+) addition, (— ) subtraction, (*) 
multiplication, (/) division, and (') exponentiation. These operators follow a 
hierarchical precedence with exponentiation operations being done first, followed by 
multiplication and division, and finally addition and subtraction. Operations like 
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multiplication and division which have the same precedence are performed from left 
to right when conflicts arise. To change the order of precedence, parentheses may 
be used around any set of operations. These parenthetical expressions have the 
highest priority and, when nested, the innermost operations within parentheses are 
done first. This scheme follows closely the protocol used in most calculators and 
high level programming languages. 

Infix notation, while convenient for the program user, does not lend itself 
well to computer manipulation. A better way to represent equations for the 
computer is the so called 'reverse Polish notation.' In reverse Polish notation, the 
operands of an equation are entered first, followed by the operator. For example, 
the infix expression 

3*4 + 5 



would be represented as 

3 4*5 + 

in reverse Polish notation. The numbers 3 and 4 are entered and multiplied, then 5 
is entered and added to the previous result. Using the concept of a 'stack' the 
reverse Polish expression is easy to evaluate. 

Recall that a stack is a last— in— first— out queue whose operation is 
analogous to a stack of trays. To operate the stack, the program calls a 'push' 
procedure to place an item on the stack, and a 'pop' procedure to remove the top 
item. Now, using the example given above, an arithmetic evaluation procedure can 
be illustrated. Figure 3.2 demonstrates such an implementation. 

The basic equation evaluation algorithm can be outlined in three steps: 

(1) Scan the reverse Polish equation term-by-term. 

(2) If the term is a constant then push it onto the stack. 
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3 




4 


(4, 3 *)= 12 


5 


(5, 12 +)= 17 


1 7 




<-stack-> 


3 


|pop4} 

{pop3| 

(multiply) 


1 2 


(pop5) 

(pop12) 

(add) 















{push 3} {push 12} (push 17} 

(push 5) 



Figure 3.2 Example of the Stack Operation 

(3) If the term is an operator then pop the first two items off the 

stack, apply the operator, and push the result back onto the 
top of the stack. 

When the algorithm is completed, the answer to the expression will be on the top of 
the stack. 

To get the infix equation into reverse Polish form is a bit more difficult 
than simply evaluating the Polish expression. Of special consideration when 
building the Polish form of the equation are the operator priorities and the use of 
parenthesis to change those priorities. A set of rules can be written which outline 
the conversion procedure [Ref. 2]. These rules are discussed below with an 
illustrative eaxmple. The infix expression 

8 + (7-6/3) *2 

will be parsed using the following operator priority table. 

Operator Priority 

4 

*7 3 

2 

operand 1 

(,),space 0 
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Using a result string called RPN, and an operator stack the rules for infix to reverse 
Polish conversion follow: 

(1) If an operand is encountered, move it to RPN. 

(2) If an operator is encountered, move all higher priority 
operators on the stack to RPN and push the new operator 
onto the stack. 

(3) If a left parenthesis is encountered push it onto the stack. 

(4) If a right parenthesis is encountered, pop all operators off the 
stack and append them to RPN until a left parenthesis is 
encountered. Discard both parenthesis. 

(5) When finished with the infix expression, pop all remaining 
operators from the stack and append them onto RPN. 

Applying these rules to the example problem above is shown in Figure 3.3. 



PRN 


STACK 


INFIX 


RULE 






8 + { 7 - 6 / 3 ) * 2 
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8 7 6 3 / - 2 * + 









Figure 3.3 Conversion from Infix to Reverse Polish 
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2. NextGainl Function 

This function makes two kind of intervals, linear interval and 
logarithmic interval. This function is called by the 'GetRootl' procedure. 

3. NextGain2 Function 

This is like the 'NextGainl' function except it has two parameters. It 
supports two kind of intervals for each parameter. This function is called by the 
'GetRoot2' procedure. 

4. Results Procedure 

This procedure outputs the calculated roots to the device 'Out File'. 
Then it makes it possible for the user to access the numerical data from the hard 
disk. Also, it supports the plot array named 'GraphArray' for the plotting data. 

5. Plot Root Locus 1 and PIotRootIjOCus2 Procedure 

The structure of these procedures is exactly the same except for the 
auto— scale function of the coordinate to be supported by the 'PlotRootLocusl.' 
These procedures draw the roots and some information in the desired coordinate. 
First, the 'SelectWind' is called with the integer number and a boolean expression in 
order to select a window as visible. The 'Define Header' is called to draw the title. 

Next the 'OpenPic' opens a picture for a specific window and only shows 
the drawing if the boolean is set True. In order to defind the coordinate it calls 
'FindWorld' for the auto— scale or 'FindWorldl' for the manual— scale. Then 
'DrawAxis' draws an axis with Footers and optional arrows on the axis. Finally, 
'DrawPolygon' draws a polygon defined in the plot array with the predefined shape. 

D. ROOTFINDER UNIT 

The unit which solves for the roots of a polynomial is called 'RootFinder.' 
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This unit uses Laguerre's method with linear deflation. Since this method is a 
commercially available package of subroutines in the 'Turbo Pascal ToolBox 
(Numerical Methods)', a brief explanation is offered here. 

1. Laguerre's Method 

Laguerre's method attempts to approximate all the real and complex 
roots of a real or complex polynomial. Laguerre's method is very reliable and quick, 
even when converging to a multiple root. 

To motivate (although not rigorously derive) the Laguerre formulas we 
can note the following relations between the polynomial and its roots and 
derivatives. 



Pn(x) = (x - Xi)(x - X 2 ) . . . (x - Xn) 



(3.1) 



ln|Pn(x)| = ln|x - xi| + ln|x - X 2 I + . . . 
+ ln|x -Xn| 



(3.2) 



dln|Pn(x)| L 



"dF 



+ 



X - Xi X - X2 

I ^ _ Pn - p 

+ , . _ — p = Li 

X Xn r n 



+ 



(3.3) 



d^ln|Pn(x)| _ 1 I 1 

dx^ (x - xi)- (x - X 2 )' 

1 



+ 



(x - X„)2 



Pn 



A n 



(3.4) 



Starting from these relations, the Laguerre formulas make what Acton nicely calls 
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"a rather drastic set of assumptions", the root xl that we seek is assumed to be 
located some distance a from our current guess x, while all other roots are assumed 
to be located at a distance b. 



a = x — Xl b = x— xj i = 2, 3,. . n 



(3.5) 



Then we can express Eqs (3.3), and (3.4) as 



1 I n 1 p 

a + 5 



(3.6) 



— +— —= H 



a^ 



b2 



(3.7) 



which yields as the solution for a 

a = ^ (3.8) 

G ± '^(n - l)(nH - G^) 

where the sign should be taken to yield the largest magnitude for the denominator. 
Since the factor inside the square root can be negative, a can be complex. (A more 
rigorous justification of Eq. (3.8) is in [Ref. 3].) 

The method operates iteratively; for a trial value x, a is calculated by 
Eq. (3.8). Then x — a becomes the next trial value. This continues until a is 
sufficiently small. In the next section, a major procedure which handles this method 
will be explained. 
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2. InitAndTest Procedure 

This procedure sets the initial value of the input and output variables. 
This procedure also tests the tolerance (Tol), maximum number of iterations 
(Maxiter), and code. Finally, it examines the coefficients of Poly. If the constant 
term is zero, then zero is one of the roots and the polynomial is deflated accordingly. 
Also if the leading coefficient is zero, the degree is reduced until the leading 
coefficient is non— zreo. 

3. FindOneRoot Procedure 

This procedure approximates a single root of the polynomial Poly. The 
root must be approximated within Maxiter iterations to a tolerance of Tol. The 
root, a value of the polynomial at the root, and the number of iterations(Iter) are 
returned. If no root is found, the appropriate error code (Error) is returned. 

4. EvaluatePoIy Procedure 

This procedure applies the technique of synthetic division to determine 
the value (yValue), first derivative (yPrime) and second derivative (yDoublePrime) 
of the polynomial. Poly, at X. The 0th element of the first synthetic division is the 
value of Poly at X, the 1st element of the second synthetic division is the first 
derivative of Poly at X, and twice the 2nd element of the third synthetic division is 
the second derivative of Poly at X. 

5. ConstructDifference Procedure 

This procedure computes the difference between approximations; given 
information about the function and its first two derivatives. 

6. TestForRoot Function 

These are the stopping criteria. Four different ones are provided. If you 
wish to change the active criteria, simply comment off the current criteria 
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(including the appropriate OR) and remove the comment brackets from the criteria 
(including the appropriate OR) you wish to be active. 

7. ReducePoly Procedure 

This procedure deflates the polynomial Poly by factoring out the Root. 
Degree is reduced by one. 

E. MESSAGE UNIT 

This unit provides the several messages which inform the user with some 
warnings and help informations. There are 11 procedures, which have the same 
structure, to provide the content of a message and one procedure, 'SetupWindow,' 
to define the window. Here, the 'SetupWindow' procedure and one sample 
procedure to make message will be explained. 

1. SetupWindow Procedure 

This procedure defines all of the windows in the program. It calls the 
'DefineWindow' procedure from 'TurboGraph' Unit. There are five other standard 
window types: documentPrc, DocProc, dBoxProc, PlainDBox, and noGrowDocProc 
[Ref 4: P. 463]. We can insert window type, window ID, and size of window in this 
procedure. 

2. MakeInfoScreen Procedure 

This is one of nine procedures for messages. It is a sample to explain 
this kinds of procedure. 'MakeInfoScreen' procedure brings up a window with a 
description of what this program is. First, a 'GrafPort' is called. It is simply a 
Pascal record type with fields that control QuickDraw's behavior. [Ref 4; p. 405 — 
422] 'GrafPort' allows an application with windows to invoke drawing operations 
that are appropriate for each window. The 'SelectWind' chooses the window to be 
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defined in 'SetupWindow' using just window ID for this procedure. The selected 
window moves to the center of the screen using the 'Movewindow' and the 
'SetVisibility' sets the visibility of a window. The 'TextFont,' 'TextSize,' and 
'TextStyle' are used to define the letter to be drawn in the window. Finally the 
'MoveTo' assigns the location to be drawn. 'DrawString' draws the messages in the 
window. 

F. MYDIALOG UNIT 

This unit supports 4 dialog boxes for input data. Each dialog procedure has 
the same skeleton to make a program; the only major procedure to construct the 
program will be explained in this section. 

A dialog box communicates with the background text, controls, icons, and 
'editText' items in which the user can enter and edit text; the user talks back with 
the mouse and keyboard. 

The Dialog Manager [Ref. 5: p. 53] handles this communication. The Dialog 
manager leans heavily on the resource mechanism [Ref. 5: pp. 137 — 154] to provide 
it with data structures. Basically, each dialog is represented on disk as a resource of 
'resType' DLOG— a template describing the dialog's size, window type, and 
title— and a resource of type DITL (Dialog Item List), which lists the content of the 
dialog (controls background text, and so on). 

A program that intends to use dialogs sees to it that appropriate resources of 
type DLOG are available at runtime. Then, the program loads them into memory 
and draws them on the screen with 'GetNewDialog.' 

After 'GetNewDialog' has read a dialog template and its item list into memory 
and drawn it on the screen, a program enters a loop in which it repeatedly 
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calls the lynchpin of dialog processing, 'ModalDialog,' and acts on the integer value 
it returns. 'ModalDialog' allows the user to customize its operation with a routine 
specified by the user. The user informs 'ModalDialog' of this by passing a specific 
parameter that points to the routine, which then gets control every time an event is 
generated when the dialog is on the screen. It is up to the filter routine to decide 
what to do about it. 

To set the value of a control, you need to call 'SetCtlValue'. Since the 
'SetCtlValue' expects a handle to the control record, you will have to get a handle 
to the control that needs setting. This is a task for the Dialog Manager's 
'GetDItem' routine. 'GetDItem' takes in two bits of information and returns three. 
Given the indicated dialog and item number, it returns the information about that 
particular item: its type (such as 'radioButton' or 'staticText', encoded as an 
integer), a handle to its underlying data structure (for controls, this is a 
'ControlHandle'), and finally, its bounding box. The key working with 'editText' 
items is to use 'GetIText' procedures. An 'editText' item begins life with the 
starting content assigned by its definition in an RMaker file [Ref 5: p. 8]. After the 
user has played with it (as determined by a suitable 'itemHit' value returned by 
'ModalDialog'), 'GetIText' is used to see what the value is now. 

G. TURBOGRAPH UNIT 

This unit supports many procedures for graphics under the Turbo Pascal 
environment. It has commercially available packages in the Turbo Pascal Toolbox 
(Numerical Method). 

Many procedures and functions in this unit are called by several units. Since 
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these routines are already explained, they will not be explained again. The 
comments in source code of this unit give you enough information. 
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IV. EXAMPLE OF DESIGN 



A. OVERVIEW 

The root locus method is very valuable in the analysis of dynamic systems, 
and is also used for design. By inspection of the curve, we determine: 

(1) Whether any loci cross the ]u) axis into the right half s— plane ( If such a 
crossing exists, it defines a stability limit for the system.), 

(2) The frequency (value of uJ) at such a stability limit, 

(3) Whether the loci go through any area on the s— plane where we want 
dominant complex roots for our system (i.e., can we get what we want). 

We can also determine the value of the parameter at the stability limit, the 
value of the roots for any specified value of the parameter, and the value of the 
parameter required to place roots at any selected point on the loci. These latter 
items, however, are not done by inspection of the curves, but by inspection of the 
computer printout or by separate calculations. On the MacRootLocus, both plot 
and calculation data are supported. Since two parameters of the system are 
adjustable, it is possible to calculate and plot a family of root loci for the pair of 
parameters The technique for synthesizing a system utilizing the rootlocus method 
is demonstrated in this chapter. 

B. GRAPHICAL SOLUTION 

In this section, four different kinds of control problems will be demonstrated. 
Such problems include simple cascade compensators for systems, subject only to 
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step inputs and/or load disturbances, and feedback compensators with no more than 
two adjustable variables. 

1. Example 1 (cascade lead compensation) 

The block diagram of Figure 4.1 shows the general case of cascade 
compensation. 




Figure 4.1 Simple Casade Compensation Block Diagram 



Consider a plant with 



Gp(S) 



400(S + 1 ) 

S2(S + 1)(S + 100) 



(4.1) 



The uncompensated system is unstable. It is desired to stabilize the system with a 
low— pass filter to reduce bandwidth. Specifications are now in the form of a desired 
location for the roots of the characteristic equation. A cascade lead compensator is 
to be used, and we choose the simplest possible i.e.. 
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Gc(S) 



P (S + Z) 
Z (S + P) 



Also we define A and B as 

A = — , B = P. 

Z 



The characteristic equation for this system is 



+ (101 +B)S^ + (100 + 101B)S^ + (lOOB + 400A)S2 
+ 400A(1 + B/A)S + 400B = 0. 



The two parameter root locus family for this system is given on Fig 4.2 
desired roots are 

S = 1.684 ± j 3.579, 

an appropriate choice from this plot is A = 6.0, B = 4.8, from which Z 
P = 5.0. 



(4.2) 



(4.3) 



(4.4) 

Since the 



0.96, and 
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Figure 4.2 Plot of Example 1 
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2. Example 2 (cascade lag compensation) 

If we consider a type— one third-order plant for which 



Gp(s) = 1^5 


(4.5) 


S(S + 1)(S + 10) 


Gc(S) = P 


(4.6) 


Z (S + P) 



Then, letting 

A=-^, B = P, (4.7) 

Z 

the characteristic equation for the cascade lag compensated system is 

+ (11 + B)S3 + (10 + 11B)S2 + (lOB + 150A)S 

+ 150B = 0. (4.8) 



The two parameter root Locus for the lag relocation zone is shown on 
Figure 4.3. Choosing A = 0.1 and B = 0.1 gives P = 0.01, Z = 0.1 and provides 
complex roots at 

S = 0.0372 ±j 1.12, 
with real roots at S = —1.014, —10.15. 

Stabilization can also be achieved with a lead filter, but the resulting 
system will have a wider bandwidth and faster response. 
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Figure 4.3 Plot of Example 2 
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3. Example 3 (velocity feedback) 

The basic problem is given by the block diagram of Figure 4.4. G(s) 
may be any order. 




Figure 4.4 Velocity Feedback Compensation Block Diagram 



For a second— order system, let 



G(s) = 100 / S(S+2). 



We want the root loci for 



100 AS _ ^ 
S(S + 2) + 100 



(4.9) 



( 4 . 10 ) 



SO the characteristic equation is 
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+ (2 + 100A)S + 100 = 0. 



(4.11) 



The resulting loci are shown on Fig 4.5. Inspection of Fig 4.5 shows that increasing 
the feedback gain increases damping. Any desired value of C is available for the 
complex roots. With high feedback gain, overdamping (all real roots) is available, 
and no positive value of gain can make the system unstable. To design the system, 
pick a root location on the locus, use the magnitude rule to find the gain as in 
Figure 4.5, or, since a computer program was used to generate the locus, the 
tabulated data are supplied to get the desired information. Part of the tabulated 
data are shown below. 



A = 0.07971 

-4.98533003622524e+0 
-4.98533003622524e+0 

A = 0.08633 

-5.31655371460862e+0 
-5.31655371460862e+0 

A = 0.09351 

-5.67530563384649e+0 
-5.67530563384649e+0 

A = 0.10128 

-6.06387368606130e+0 
-6.06387368606130e+0 

A = 0.10969 

-6.48473591175408e+0 
-6. 484735911754086+0 

A = 0.11881 

-6.940576303174546+0 
-6.940576303174546+0 

A = 0.12869 

-7.434301921120686+0 
-7.434301921120686+0 



+ 8.668707194842296+0 j 
+-8.668707194842296+0 j 



+ 8.469607818528636+0 j 
+-8.469607818528636+0 j 



+ 8.233523301869646+0 j 
+-8.233523301869646+0 j 



+ 7.951693902527526+0 j 
+-7.951693902527526+0 j 



+ 7.612371519756976+0 j 
+-7.612371519756976+0 j 



+ 7.199194439644766+0 j 
+-7.199194439644766+0 j 



+ 6.688135386310706+0 j 
+-6.688135386310706+0 j 
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Figure 4.5 Plot of Example 3 
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4. Example 4 (velocity and acceleration feedback) 

The block diagram of Figure 4.6 feedback shows the general case of 
velocity and acceleration feedback. 




Figure 4.6 Velocity and Acceleration Feedback Compensation Block Diagram 



If we consider a third order plant for which 



Gp(s) 



8 

(S + 1)3 ’ 



(4.12) 



then the root locus form is 



8(AS + BS3) _ ^ 
(S + 1)3 + 8 



and the characteristic equation is 



49 



S3 + (3 + 8A)S3 + (3 + 8B)S + 9 = 0. 



(4.14) 



The two— parameter root locus family is given on Figure 4.7. We may select 
desired location for the complex roots such as s = —0.73^ j 1.07, for which Ai 
0.48, B 2 =0.8. The real root is at s = — 5.38 so the complex roots are dominant. 
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Figure 4.7 Plot of Example 4 
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V. CONCLUSION AND RECOMMENDATION 



MacRootLocus is a useful and powerful tool for the designer and a simple and 
easy to use package for the user. The program supports the two parameter root 
locus method intensively as well as the one parameter root locus method. Also, it 
offers high resolution plots, the tabulated data and a small text editor for recording 
some information of the design. 

The program does, however, have some areas where improvements are 
possible. These are discussed below. 

(1) According to the engineer's needs, the root finder can be replaced 
easily with a more reliable and faster algorithm, and 

(2) the capability to move the information window can be added. 

Since this program was developed under the standard interface philosophy the 
Macintosh was designed for, it is a simple matter to change and to append some 
subroutines. 

MacRootLocus is the user friendly CAD program available that is both simple 
enough for the beginning student to easily use as well as powerful and flexible 
enough to benefit the experienced system designer. 
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APPENDIX 



SOURCE CODE 

This appendix contains the source code of all modules in MacRootLocus 
except Standard Apple Macintosh libraries. Some subroutines in Numerical Method 
(Turbo Pascal Tool Box) are also included since they are slightly modified. A disk 
is available from Dr. Thaler that contains the MacRootLocus source code and the 
MacRootLocus resource file. 
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{= =} 

{= MacRootLocus = } 

{= version 1.0 =} 

{= =} 

{= 04 DEC 1989 =} 

{= =} 

{= author = } 

{= KO, SUNG HOON = } 

{= =} 

{= =} 

{= SYSTEM : Apple Macintosh SE = } 

{= LANGUAGE : Turtxj Pascal for Macintosh versioni .0 =} 

{= LIBRARIES : Numerical Method (Turbo Pascal Tool Box) = } 

{= RESOURCE : MacRootLocus.rsrc =} 

{= =} 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 

program MacRootLocus; 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 

{= This main program handles the system and integrate theresource =} 

{= and 6 units. =} 

1 = = = = = ^ = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 

{$B+} { Set the bundle bit } 

{$R MacRootLocus. Rsrc} { Identify resource file for menu and icon info } 

{$T APPLFFTD} { Set the application type and creator } 

{$S+} { Generate segmented code} 

{$!-} { Turn off I/O error checking } 

{$U SpecVar) 

{$U RootsFinder} 

{$U MakeRoot) 

{$U TurboGraph} 

{$U Message) 

{$U MyDialog) 

uses 

MemTypes, QuickDraw, OSIntf, ToolIntf.PackIntf, 

PasPrinter.SANE.MacPrint, RootsFinder, SpecVar, 

{$S Second Segment) 

TurboGraph, 

Message, 

{$S Third Segment) 

MakeRoot, 

MyDialog; 
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function UpCase(Ch : char) : char; 



{=========== ===} 

{= Returns the upper case of Ch =} 

inline 



$301 F. { UpCase MOVE.W (SP)+,D0 ; GetCh } 
S0C40, 

$0061 , { CMP.W #’a’,D0 ; skip if not lower case } 
$6D0A, { BLT.S @1 } 

$0C40, 

$007A, { CMP.W #’z’,D0 } 

$6E04, { BGT.S @1 } 

$0440, 

$0020, { SUB.W #$20,D0 } 

$3E80; { @1 MOVE.W D0,(SP) } 



procedure HideAllWindows; 

{= = = = = = = = = = = = = = = = =} 
{= Hide the three main windows =} 

{= = = = = = = = = = = = = = = = =) 
var 

WindNum : integer; 
begin 



for WindNum := RegendBox downto InfoScrnWind do 
SetVisibility(WindNum, FALSE); 
end; { HideAllWindows } 

procedure ShowAIIWindows; 



{= = = = = = = = = = = = = = = = = = = = = =} 
{= Make the three main windows visible =} 

= = = = ^ = = = = = = = = = = = = = = = =} 



var 

WindNum : integer; 
begin 

for WindNum := RootLocusWind to RegendBox do 
SelectWind(WindNum, TRUE); 
end; { ShowAIIWindows } 

procedure MakeWatchCursor; 
var 

CursorHandle : CursHandle; 
begin 

CursorHandle := GetCursor(WatchCursor); 
SetCursor(CursorHandle'^''); 
end; { MakeWatchCursor } 

procedure MakeArrowCursor; 
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begin 

InitCursor; 

end; { MakeArrowCursor } 
procedure PrintScreen; 



= = = = = = = = = = = = = = = = = =} 

{= Dump the entire screen to the printer. =} 

{= = = = = = = = = = = = = = = = = } 



begin 

HideCursor; 
HardCopy(FALSE); 
ShowCursor; 
end; { PrintScreen } 



procedure DoDeskAcc(ltem : Integer); 

{= = = = = = = = = = = = = = = = = = = = = = = =} 

{= start up desk accessory from Apple menu =} 

{= = = = = = = = = = = = = = = = = = = = = = = ^ 
var 

SavePort : GrafPtr; 

RefNum : integer; 

□Name : string; 
begin 

GelPort(SavePort); { save port before starting it } 

Getltem(MenuUst[AM], Item, DName); { get name of desk accessory } 
RefNum := OpenDeskAcc(DName); { and start that sucker up! } 
SetPort(SavePort): { restore grafport and continue } 

end; { DoDeskAcc } 

procedure SetltemState(Mndx,lndx : Integer; Flag : Boolean); 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 

{= if Flag is true, enables item Indx of menu Mndx; else disables =} 

= = = = = = _ = = = = = = = = = = = = = = = = = = = = = = = = = ^ 

begin 

if Flag then 

Enableltem (MenuList[Mndx],lndx) 
else 

Disableltem(MenuList[Mndx],lndx); 
end; { SetItemState } 

procedure HandleMenu(Menulnfo : Longint); 

{= Decode Menuinfo and carry out command =} 
var 

Menu : Integer; { menu number that was selected } 
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Item : Integer; { item in menu that was selected } 

WindNum : Integer; 
begin 

if Menuinfo <> 0 then 
begin 

PenNormal; { set the pen back to normal } 

Menu := HiWord(Menulnfo); { find which menu the command is in } 
Item := LoWord(Menulnfo); { get the command number } 

case Menu of { and carry it out } 

AppIMenu : if Item = 1 then 

DoAbout { bring up "About..." window} 

else 

DoDeskAcc(ltem); { start desk accessory } 

FileMenu : case Item of { File menu } 

1 ; GetEQParameter; 

2 : GetCoeff; 

3 ; PrintScreen; {Print screen} 

4 : begin 

HideCursor; 

HardCopy(TRUE): { Print top window } 

ShowCursor; 

end; 

5 : Finished := TRUE; { Quit command } 
end; 



EditMenu : case Item of { Edit Menu } 

1..5 : if not SystemEdit(ltem-l) then 



{ Do nothing } 

end: 

PlotMenu : case Item of 

1 : PlotOneParameter; 

2 : PlotTwoParameter; 
end; 



{ bring up the dialog } 

{ for one parameter plot data} 

{ for two parameter plot data} 



HelpMenu : 

case Item of { Help menu } 

1 : InfoGetEQParameter; 

2 : InfoGetCoeff; 

3 : InfoPlotOneParameter; 

4 : InfoPlotTwoParameter; 

5 : InfoPrint 



end;{ case } 

HiliteMenu(O); { reset menu bar } 
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end; 

end; { HandleMenu } 

procedure SelectOurWindow(WPtr : WindowPtr); 

{= Select the window pointed to by WPtr if it’s one of our’s =} 
begin 

if (WPtr = Window[RootLocusWind].P) then 
SelectWind(RootLocusWind, TRUE); 
end; { SelectOurWindow } 

procedure HandleClick(WPtr : WindowPtr; MLoc : Point); 

{= Handle a mouse click within a window =} 
begin 

if WPtr <> FrontWindow then { if it’s not in front... } 
SelectOurWindow(WPtr); 
end; { HandleClick } 

procedure HandleGoAway(WPtr : WindowPtr; MLoc : Point); 

{= Handle a mouse click in the go-away box of a window =} 
var 

WPeek : WindowPeek; { for looking at windows } 
begin 

if WPtr = FrontWindow then { if it’s the active window } 
begin 

WPeek := WindowPeek(WPtr); { peek at the window } 
if TrackGo Away( WPtr, MLoc) then { and the box is clicked } 
begin 

if WPeek''. WindowKind = userKind then { if it’s our window } 
HideAllWindows { time to stop } 

else 

C!oseDeskAcc(WPeek''.WindowKind) { close DeskAcc } 
end 
end 
else 

SelectOurWindow(WPtr); 
end; { HandleGoAway } 

procedure HandleZoom(WPtr : WindowPtr; MLoc : Point; WLoc : integer); 
{= Handle a mouse click in zoom box of a window =} 
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var 

WPeek : WindowPeek; { for looking at windows } 
begin 

if WPtr = FrontWindow then { if it’s the active window } 
begin 

WPeek := WindowPeek(WPtr); { peek at the window } 
if TrackBox(WPtr, MLoc, WI_oc) then { and the box is clicked } 
begin 

if WPeek'^-WindowKind = userKind then { if it’s our window } 
ZoomWindow(WPtr, WLoc, FALSE); 
end 
end 
else 

SelectOurWindow(WPtr): 
end; { HandleZoom } 

procedure HandleGrow(WPtr : WindowPtr; MLoc : Point); 



{= = = } 

{= Handle mouse click in the grow box of a window =} 



type 

GrowRec = record 

case integer of 

0 : (Result : Longint); 

1 : (Height, Width : integer); 

end; 

var 

Growinfo : GrowRec; 

WindNum : integer; 
begin 

if WPtr = FrontWindow then { if it’s the active window } 
with Growinfo do 
begin 

Result := GrowWindow(WPtr, MLoc, GrowArea); { get amount of growth } 
SizeWindow(WPtr, Width, Height, TRUE); { resize window } 
lnvalRect(WPtr''.portRect); { set up for update } 

WindNum := 2 ; 

if Window[WindNum].P <> FrontWindow then 

DrawGrowlcon(Window[WindNum].P); { draw grow icons } 
end 
else 

SelectOurWindow(WPtr) ; 
end; { HandleGrow } 

procedure HandleDrag(WPtr : WindowPtr; MLoc : Point); 

{= Handle the dragging of a window =} 
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var 

WindNum : integer; 
begin 

if WPtr = FrontWindow then 
begin 

DragWindow(WPtr, MLoc, DragArea); { in the drag bar } 
WindNum := 2 ; 

DrawGrowlcon(Window[WindNum].P); { draw grow icons } 
end 
else 

SelectOurWindow(WPtr); 
end; { HandleDrag } 

procedure DoMouseDown(TheEvent : EventRecord); 



{= = = = = = } 

{= identify where the mouse was clicked and handle it =} 

{= = = = = = = = = = = = = = = = = = = = = = = = =} 



var 

theWindow : 

WindowPtr; 

MLoc : 

Point; 

WLoc ; 
integer; 
begin 

MLoc ;= TheEvent.Where; { get mouse position } 

WLoc := FindWindow(MLoc, theWindow) ;{ get the window and the location } 

case WLoc of 

InMenuBar : HandleMenu(MenuSelect(MLoc)); { in the menu } 
InContent : HandleClick(theWindow, MLoc); { inside the window } 
InZoomln : HandleZoom(theWindow, MLoc, WLoc);{ in the zoom box } 
InZoomOut : HandleZoom(theWindow, MLoc, WLoc);{ in the zoom box } 
InGoAway : HandleGoAway(theWindow, MLoc); { in the go away box } 
InGrow : HandleGrow(theWindow, MLoc); { in the grow box } 
InDrag 

HandleDrag(theWindow, MLoc); { in the drag bar } 

InSysWindow ; SystemClick(TheEvent, theWindow); { in a DA window } 
end; 
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end; { DoMouseDown } 

procedure DoUpdate(TheEvent : EventRecord); 



{= = = = = ^ = = = = =) 

{= handles window update event =} 

{= = = = = = = = = = = = =} 

var 

SavePort, 

the Window : WindowPtr; 
begin 



theWindow := WindowPtr(TheEvent. Message): { find which window } 

if (theWindow = Window[RootLocusWind].P) or 
(theWindow = Window[RegendBox].P) 
then { only update our windows } 

begin 

MakeWatchCursor; 

GetPort(SavePort): { save current grafport } 

SetPort(the Window); { set as current port } 

BeginUpdate(theWindow); { signal start of update } 

{ and here’s the update stuff! } 

if theWindow = Window[RootLocusWind].P then 
begin 

C le arWi ndo w( Root Locu s Wind) ; 

DrawGrowlcon(theWindow); 

DrawPic(RootLocusWind); 

end; 

{ now, back to our progrann...} 

EndUpdate(theWindow): { signal end of update } 

SetPort(SavePort): { restore grafport } 

MakeArrowCursor; { restore cursor } 

end 

end; { DoUpdate } 

procedure DoActivate(TheEvent : EventRecord): 

{= Handles window activation event =} 
var 

AFlag : boolean; 
theWindow : WindowPtr; 

WindNum : integer; 
begin 

with TheEvent do 
begin 
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theWindow := WindowPtr(Message); { get the window } 

AFlag := Odd(Modifiers); { get activate/deactive } 

if AFlag then 

begin { if it's activated... } 

SetPort(the Window); { make it the port } 

WindNum := RootLocusWind ; 

DrawGrowlcon(Window[WindNuml.P); { draw grow icons } 
end; 
end 

end; { DoActivate } 

procedure DoKeypress(theEvent : EventRecord); 



{= handles keypress (keyDown, autoKey) event =} 

= _ = = = = = = = = = 

var 



KeyCh : char; 
begin 

KeyCh := Chr(theEvent. Message and charCodeMask); { decode character } 
if (theEvent.modifiers and cmdKey) <> 0 then 
begin { menu key command } 

HandleMenu(MenuKey(KeyCh)) { get menu and item} 

end 
else 

if TextInputEnabled then 
begin 

TEKey(KeyCh, textH); 

TEUpdate(thePort''.portRect, textH); 
end 
else 

SysBeep(t); { do ’something* } 

end; { DoKeypress } 

procedure Initialize; 



= = = = = = = = = = = = = = ^ = = = = = = = = 

{= Initialize everything for the program =} 
= = = = = = = = = = = = = = = = = = = = = 
var 



Indx 

integer; 

Result : real; 
FileErr : byte; 
begin 

lnitGraf(@thePort); 

InitFonts; 

I nit Windows; 



{ initialize all managers used } 

{ create a grafport for the screen } 
{ start up the font manager } 

{ start up the window manager } 
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InitMenus; { start up the menu manager } 

TEInit; 

InitGraphic; 

SetUpWindows; 

TextFont(SystemFont): { initialize the font } 

{ set up menus } 

MenuList[AM] := GetMenu(ApplMenu);{ read menus in from resource fork } 
MenuList[FM] := GetMenu(FileMenu); 

MenuListfEM] := GetMenu(EditMenu); 

MenuList[PM] := GetMenu(PlotMenu); 

MenuListfHM] := GetMenu(HelpMenu); 

AddResMenu(MenuList{AM],'DRVR'): { pull in all desk accessories } 
for Indx := 1 to 5 do { place menus in menu bar } 
lnsertMenu(MenuList(lndx],0); 
for Indx := 1 to 6 do 

SetltemState(EM, Indx, FALSE): { deactivate items in Edit menu } 

DrawMenuBar; { draw updated menu bar to screen } 

Finished := False; { set program terminator to false } 

{ set drag region } 

SetRect(DragArea, XMinGlb+1, YMinGlb+38, XMaxGlb-1, YMaxGlb-1); 

{ set grow region } 

SetRect(GrowArea, XMinGlb+1, YMinGlb+38, XMaxGlb-1, YMaxGlb-1): 
InitGuess.Re := 1.0; 

InitGuess.lm := 0.0; 

Tolerance := 1e-6; 

Maxiter := 100; 

InitDegree ;= 0; 

AMinGain := 0.1 ; 

AMaxGain ;= 10000; 

BMinGain := 0.1 ; 

BMaxGain := 10000; 

XMn :=-10; 

XMx ;= 5; 

YMn := -10; 

YMx := 10; 

Step := 4; 

Points := 50; 

FillChar(lnitPoly, SizeOf(lnitPoly), 0); 

FillChar(xAnswer, SizeOf(xAnswer), 0); 

FillChar(xlnitPoly, SizeOf(xlnitPoly), 0); 

FillChar(lnfixArray, SizeOf(lnfixArray), 0); 

MakeWatch Cursor; 

MakeInfoScreen; 
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MakeArrowCursor; 

TextInputEnabled := False; 

InitDegreeStatus := False; 

GetCoeffStatus := False; 

end; { Initialize } 

procedure HandleEvent(TheEvent : EventRecord); 

{= Decodes an event and handles it =} 
begin 

case The Event. What of 
mouse Down 

DoMouseDown(TheEvent); { mouse button pushed } 

key Down ; DoKeyPress(TheEvent); { key pressed down } 

autoKey : DoKeyPress(TheEvent); {key held down } 
updateEvt : DoUf^ate(TheEvent); { window needs updating } 
(activateEvt : DoActivate(TheEvent);} { window made act/inact } 
end 

end; { HandleEvent } 
procedure Cleanup; 



{= = = = = = = } 

{= Do any last minute clean up work =} 
= ^ = ^ = = = = = = = = = = = 



begin 

end; { Cleanup } 



begin { MacRootLocus } 

Initialize; { set everything up } 

repeat { keep doing the following } 

SystemTask; { update desk accessories } 

if GetNextEvent(everyEvent.theEvent) then { if there’s an event... } 
HandleEvent(theEvent); { ...then handle it } 

until Finished; { until user is done } 

Cleanup; 

end. { MacRootLocus } 
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unit GlobalVar (5000); 

{= This unit declares the whole variables tobe used =} 
{= in MacRootLocus except a part of local variables. =} 



{$U RootsFinder} 

interface 

uses 

MemTypes, QuickDraw, OSIntf, ToolIntf.PackIntf.PasPrinter, 
SANE.MacPrint, RootsFinder; 



const 

InfoScrnWind = 1; 
RootLocusWind = 2; 
RegendBox = 3; 
AboutBoxWind = 4; 
AlertBox = 5; 
HeIpWind = 6; 



{ Window number of info window } 

{ Window number of root locus window } 

{ Window number of informatiom box window } 
{ Window number of about box window } 

{ Window number of alert box window } 

{ Window number of help window) 



MenuCnt = 5; { total # of menus } 

AppIMenu = 1000; { resource ID of Apple Menu } 

FileMenu = 1001; { resource ID of File Menu } 

EditMenu = 1002; { resource ID of Edit Menu } 

PlotMenu = 1003; { resource ID of Plot Menu } 

HelpMenu = 1004; { resource ID of Help Menu } 



AM 

FM 

EM 

PM 

HM 



1 ; { index into MenuList for Apple Menu } 



2; { ditto for File Menu } 

3; { ditto for Edit Menu } 

4; { ditto for Plot Menu } 

5; { ditto for Help Menu } 



MaxPlotGIb 



= 1024;{ The maximum number of points in a plot array } 



type 

str255 = string [255] ; 

StringArray = array[1..20] of str255; { Storage for plot data } 
NodePtr= ''Node; 

Node = 



record 

Data : Extended; 

Next : NodePtr; 
end; 

PlotArray = arrayfl ..MaxPlotGIb, 1..2] of Extended; 
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var 

Answer, polish : str255; 

InfixArray : StringArray; 

A,B,ADeltaStep, AIncrem .BDeltaStep, BIncrem: Extended; 

InitGuess ; TNComplex; { Initial approximation } 

Tolerance : Extended; { Tolerance of approximation } 

AMaxGain, BMaxGain : Extended; { Range of unknown parameter values } 
AMinGain, BMinGain : Extended; 

Root, I mag. Value, Deriv : TNvector; { Resulting roots and other info } 

Iter : TNIntVector; { Iterations to find each root } 

Maxiter ; integer; { Maximum number of iterations } 

InitDegree, Degree : integer; { Initial and final degree } 
xInitPoly : TNvector; { of polynomial } 

InitPoly, Poly : TNCompVector; { Initial and final coefficients } 

{ of the polynomial } 
xAnswer : TNCompVector; 
yAnswer : TNCompVector; 

NumRoots, Arrayindex, pointno : integer; { Number of roots } 

Error : byte; { Error flag } 

StepA, StepB, linear : Boolean; 

Step,Points : Longint; 



OutFile : text; 
OutName : string; 
lOerr : boolean; 



{ The output file used by an application} 
{ The out file name} 

{ Flags I/O errors } 



Finished : boolean; { used to terminate the program } 

theEvent : EventRecord; { event passed from operating system } 

MenuList ; array[1..MenuCnt] of MenuHandle; { holds menu info } 
DragArea : Rect ; { Area in which window can be dragged } 

GrowArea : Rect; { Area in which a window’s size can change } 



DataPicture : PicHandle; 
the Dialog ; DialogPtr; 



item Hit 


: Integer; 


theType 


: Integer; 


r ; 


Rect; 


done, AutoScale 


n 


: Integer; 


h,h3,h4,h5,h6 


s 


: Str255; 


XMn 


: Extended; 


XMx 


: Extended; 


YMn 


: Extended; 


YMx 


: Extended; 



Boolean; 

: Handle; 



{ Holds a picture of the plotted data } 



GraphArray : ''PlotArray; 
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Ds : DecStr ; 

ARJustification, AMarkStatus.BRJustification, BMarkStatus : Boolean; 
Da, Db, Dan, Dax, Dbn, Dbx, Dss : DecStr; 
txRect : Rect; 
textH :TEHandle; 

TextInputEnabled : Boolean; 

InitDegreeStatus : Boolean; 

GetCoeffStatus : Boolean; 

implementation 

begin 

end.{GlobalVar} 
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unit MakeRoot(3000); 



{= This unit provides several procedures and functions to = } 

{= make the roots and the array vectors for plot. = } 

{= = = = = = = = = = = = = = = = = = == = } 



{$!-} { Disable I/O error trapping } 

{$S+} { Enable segmentation of code } 

{$U SpecVar} 

{$U Message} 

{$U RootsFinder) 

{$U TurboGraph} 

interface 

uses 

MemTypes, QuickDraw, OSIntf, Toolintf, PackIntf.PasPrinter, 
SANE, MacPrint, RootsFinder, SpecVar, 

{$S SecondSegment} 

TurboGraph, 

Message: 

var 

Error : byte; 

procedure GetRootl ; 
procedure GetRoot2; 

implementation 



procedure lnfixToPolish( Answer:str255;var RPN;str255); 

{= = = = = = = = = = = = = = = = = = = = = = = = = =} 

{= This procedure converts an algebraic expression(infix) = } 
{= to reverse Polish notation. = } 

= = = ^ = = ^ = = = = = = = = = = = = == = = = = =} 



var 

Stack : Array[1 ..50] of char; 
Top,l,N,p : Integer: 

Ch ; char; 

Ch1 ; string; 

Firstchar,PrevDigit : Boolean; 

Function Priority(Ch:char): Integer; 
Begin 
Case Ch of 
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'''■ : Priority:= 4; 

: Priority := 3; 

'/' : Priority:= 3; 

: Priority:= 2; 

V : Priority:= 2; 

'a’..'b' : Priority:= 1; 

: Priority:= 1 ; 

: Priority:= 1 ; 

: Priority:= 1 ; 

: Priority:= 0; 

’)' : Priority:= 0; 

' ' : Priority:= 0; 

end; 
end; 

begin 

RPN := 

Top :=0; 

FirstChar :=True; 

PrevDigit := False; 

for i:= 1 to Length(Answer) do 
begin 

ch1 := copy(Answer,i,1); 
ch:= ch1(1]; 

P;= Priority (ch); 

if FirstChar and (Ch = then P:= 1 ; 

if P = 1 then 
begin 

if PrevDigit then 
RPN := concat(RPN.ch) 
else 

RPN := concat(RPN,’ *,ch); 
firstchar ;= False; PrevDigit := True; 
end; 

if P>1 then 
begin 

while (Top>0) and (Priority(stack[top]) >= P) do 
begin 

RPN := concat(RPN,' ',Stack(Top]); 

Top := Top - 1; 
end; 

Top := Top + 1; 
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StackfTop] := ch; 

Firstchar := True; PrevDigit := False; 
end; 

if ch = ’(’ then 
begin 

Top := Top +1 ; 

Stack[Top] := 

FirstChar := True; PrevDigit :=False; 
end; 

if ch = ’)’ then 
begin 

while Stack[Top]<> ’(’ do 
begin 

RPN := concat(RPN,’ ’.StackfTop]); 

Top := Top -1; 
end; 

Top := Top -1; 

Firstchar := False; PrevDigit ;= False; 
end; 
end; 

while Top > 0 do 
begin 

RPN := concat(RPN,’ '.StackfTop]); 

Top := Top -1; 
end; 

end; 

procedure ComputePolish(polish:str255;a.b:Extended; var EvalArray:Extended); 



{= = = = = = = = = = = = = = =} 

{= This procedure uses the string generated in Procedure = } 

{= InfixToPolish to evaluate the numeric expression =} 

{= = = = = = = = = = _ = = = = = = = = = = = = = = = = = = = = = =} 



var 

i.N : integer; 

ch.ch2 : char; 

ch1 : string; 

temp : string[255]; 

Value1.Value2.Value3 : Extended; 

StackPtr ; NodePtr; 

StackEmpty : Boolean; 

procedure CreateStack; 



{= = = = = = = = = = = = =} 
{= Initialize stack =} 

{= = = = = = =} 



begin 
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New(StackPtr); 
with Stack Ptr^ cJo 
begin 

Next ;= nil; 

Data := 0.0; 
end; 

StackEmpty := True; 
end; 

procedure Pop(var Val : Extended);{Push a number onto numeric stack} 



{= = = = = = = = = = = = = = = = = = = = = =} 
{= pop a number off numeric stack =} 
= = ^ = = = = = = = = = = 



var 

NPtr :NodePtr; 
begin 

if not StackEmpty then 
begin 

NPtr := StackPtr'^.Next; 

StackPtr'^.Next := NPtr'^.Next; 

Val := NPtr^.Data; 

Dispose(NPtr): 

StackEmpty := (StackPtr'^.Next = nil); 
end; 
end; 

procedure Push(Val ; Extended); 



{ = = = = = = = = = = = = } 

{= pop a number off numeric stack =} 

{= = = = = = = = = = = = = =} 



var 

NPtr ; NodePtr; 
begin 

StackEmpty := False; 

New(NPtr); 

NPtr^.Data := Val; 

NPtr'^.Next := StackPtr''.Next; 
StackPtr'^.Next := NPtr; 

end; 

procedure DeleteStack;{Delete stack} 



{= = = = = = = = = = =} 
{= Delete stack =} 

{= = = = =} 



var 

Temp : Extended; 
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begin 

while not StackEmpty do 
Pop(Temp): 
Dispose(StackPtr): 
end; 



function Expon(y,x:Extended) Extended; 



{= = = = = = } 

{= Computes Y raised to X power =} 

{= = = = = = = = = = = = = = = = = = = = =} 



begin 

Expon := exp( x * (ln(y))); 
end; 



begin 

CreateStack;{lnitialize} 
temp :=’ 

for i:= 1 to Length(polish) do{do one char at a time} 
begin 

ch1 := copy(polish,i,1): 
ch:= ch1[1]:{Get a char} 

case ch of {and evaluate it} 

*0’..’9’ ; temp :=concat(temp,ch);(Real constant} 

: temp := concat(temp,ch); 

: begin 

ch1 := copy(polish,i+1,1): 
ch2;= ch1[1}:{Get a char} 

if (ch2 <> ' ’) and (i < length(polish)) then{Unary minus} 
temp := concat(temp.ch) 
else 

begin {minus operator} 

POP( Valuel ) ;POP{Value2) ; 

Values := Value2 - Valuel ; 

PUSH(Value3): 

end; 

end; 



•a’ 


: PUSH(a): 


’A’ 


: PUSH(A): 


’b* 


: PUSH(b): 


’B' 


: PUSH(B): 
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V :begin 

POP(Valuel): POP(Value2); 
Values := Value2 + Valuel; 

PUSH( Values): 
end; 

:begin 

POP(Valuel): POP(Value2); 
Values := Values * Valuel ; 
PUSH( Values): 

end: 



V :begin 

POP(Valuel): POP( Values): 
Values := Values / Valuel : 
PUSH( Values): 
end: 



:begin 

POP(Valuel): POP(ValueS): 
Values :=Expon( Values , Valuel): 
PUSH(ValueS): 
end: 

:begin 

if temp <> ' 'then 
begin 

Valuel ;= StrSNum(temp): 
PUSH(Valuel): 
temp := ’ ’: 
end: 
end: 



end: 

end: 

if temp <> ' ' then 
begin 

Values := StrSNum(temp): 
PUSH(ValueS): 
end: 

POP(Valuel): 
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EvalArray := Valuel ; 

DeleteStack; 

end; 

function Ten2 (power: Extended) ; Extended; 
begin 

Ten2 := Exp(power * Ln(10)); 
end; 



function NextGainI : Extended; 

{= This function makes two kind of intervals for one = } 

{= parameter root locus method. = } 

{= = = = = = = = = = = = = = = = = = = = =} 

var 

GainOut, logmin, logmax : Extended; 



begin 

if linear then 
begin 

AIncrem := abs((AMaxGain - AMinGain)/(points - 1)); 
GainOut := AMinGain + AIncrem * pointno; 
end 
else 
begin 

if (AMinGain < IE-5) and (AMaxGain > IE-1) then 
logmin := -5 
else 

logmin := Ln(AMinGain)/Ln(10); 
logmax := ln(AMaxGain)/Ln(10); 

AIncrem ;= (logmax - logmin)/(Points -1); 

GainOut := Ten2(logmin + pointno * AIncrem); 
end; 

NextGainI := GainOut; 
end; 



function NextGain2 : Extended; 

= = = = = = = = = = = =} 

{= This function makes two kind of intervals for each = } 

{= parameter of two parameter root locus method. =} 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 

var 

GainOut, logmin, logmax : Extended; 



begin 

if StepA then 
begin 
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if linear then 
begin 

BIncrem := abs((BMaxGain - BMinGain)/(points - 1)); 
GainOut := BMinGain + BIncrem * pointno; 
end 
else 
begin 

if (BMinGain < 1E-5) and (BMaxGain > 1E-1) then 
logmin := -5 
else 

logmin := Ln{BMinGain)/Ln(10); 
logmax := ln(BMaxGain)/Ln(10); 

BIncrem := Oogmax - logmin)/(Points -1); 

GainOut := Ten2(logmin + pointno * BIncrem); 
end; 

NextGain2 := GainOut; 
end 
else 
begin 

if linear then 
begin 

AIncrem := abs((AMaxGain - AMinGain)/(points - 1)); 
GainOut := AMinGain + AIncrem * pointno; 
end 
else 
begin 

if (AMinGain < IE-5) and (AMaxGain > IE-1) then 
logmin := -5 
else 

logmin := Ln(AMinGain)/Ln(10); 
logmax := ln(AMaxGain)/Ln(10); 

AIncrem := (logmax - logmin)/(Points -1); 

GainOut := Ten2(logmin + pointno * AIncrem); 
end; 

NextGain2 := GainOut; 
end; 
end; 

procedure Results(A,B : Extended; 

xAnswer ; TNCompVector; 

Error : byte); 



(=========================== 

{= This procedure outputs the results to the device OutFile =} 
{= and make the array for plot data. =} 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 
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var 

Term : integer; 
begin 

if Degree > 0 then 
begin 

Writeln(OutFile,'The deflated polynomial:'); 
for Term := Degree downto 0 do 
Writeln(OutFile, 'Poly[',Term,']:', 

Poly[Term].Re:23, ’ +’, Poly[Term].lm:23,'i’); 

end; 

if Error <= 1 then 
begin 

Writeln(OutFile); 

Writeln(OutFile, ’ A = ’ . A :10:5,' B = '. B :10:5 ); 
for Term := 1 to NumRoots do 
begin 

Writeln(OutFile, xAnswer[Term].Re:23, ’ +’, xAnswer[Term].lm:23, 'i'); 
GraphArray''[Arraylndex,1] := xAnswer[Term].Re; 
GraphArray^[Arraylndex,2] := xAnswer[Term].lm; 

Arrayindex := Arrayindex + 1 ; 
end; 
end; 

end; { procedure Results } 
procedure PlotRootLocusI ; 



{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 

{= This procedure draw a plot of one parameter root locus method =} 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 



begin 

SelectWind(2, TRUE); 

DefineHeader(2, ’One Parameter RootLocus’); 
DrawGrowlcon(Window[2].P); 

OpenPic(2, TRUE); { Open up a picture for window #1 } 

{ Draw the axis and data points } 
if AuloScaie then 

FindV/orld(2, GraphArray'^, Arrayindex) 

else 

FindWorldt (2,XMn,YMn,XMx,YMx); 

DrawAxis( Da , ", True): 

DrawPolygon(GraphArray^, 1, -Arrayindex, -3, 1, 0, TRUE); 
ClosePicture; { Close the picture for window #2} 
ValidRect(Window[2].P''.portRect); 
end; { PlotPowerExpLog } 
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procedure PlotRootLocus2; 

{= This procedure draw a plot of two parameter root locus method =} 
begin 

SelectWind{2, TRUE); 

DefineHeader(2, Two Parameter RootLocus’); 

DrawGrowlcon(Window[2].P); 

OpenPic(2, TRUE); { Open up a picture for window #2 } 

{ Draw the axis and data points } 

FindWorld1(2,XMn.YMn,XMx,YMx); 

DrawAxis(Da, Db, True); 

DrawPolygon(GraphArray'', 1, -Arrayindex, -3, 1, 0, TRUE); 

ClosePicture; { Close the picture for window #2 } 
ValidRect(Window[2].P'^.portRect); 
end; { } 



procedure GetRoott ; 

{= = = = = = = = = = = = = = = = = = = = = } 

{= This procedure gets the root of the characteristic equation for =} 

{= one parameter root locus method. =} 

(= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 



var 

term : integer; 
f : DecForm; 
begin 

Num2Str(f,AMinGain,Dan); 

Num2Str(f,AMaxGain,Dax); 

Da := ’A : ’+ Dan + ’ To '+Dax; 

Ds := ’ ’; 

Arrayindex := 1; 

New (Graph Array); 
pointno := 0; 

while pointno <= ( Points - 1) do 
begin 

A := NextGaint; 

B := 0; 

for Term := InitDegree downto 0 do 
begin 

infixToPolish(lnfixArray[Term], polish); 
ComputePolish(polish,a,b,xlnitPoly[Term]); 
lnitPoly[Term].Re := xlnitPoly[TermJ; 
lnitPoly[Term].lm := 0.0; 
end; 

Degree := InitDegree; 

Poly := InitPoly; 
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Laguerre(Degree, Poly, InitGuess, Tolerance, Maxiter, 
NumRoots, xAnswer, yAnswer, Iter, Error); 
Results(A, B, xAnswer, Error); 
pointno := pointno + 1 ; 
end; 

Arrayindex := Arrayindex - 1; 

PlotRootLocusI ; 
end; 

procedure GetRoot2{(var StepA : boolean; A,B : Extended)}; 



{= This procedure gets the root of the characteristic equation for =} 
{= two parameter root locus method. =} 

!= = = = _ = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 



var 

i, j, k, term : integer; 
f : DecForm; 
begin 

ADeltaStep := abs((AMaxGain - AMinGain)/(Step)); 
BDeltaStep := abs({BMaxGain - BMinGain)/(Step)); 
Num2Str(f,AMinGain,Dan); 
Num2Str(f,AMaxGain,Dax); 

Da := ’A : ’+ Dan + ’ To ’+Dax; 
Num2Str(f,BMinGain,Dbn) ; 
Num2Str(f,BMaxGain,Dbx) ; 

Db := ’B ; ’+ Dbn + ’ To ’+Dbx; 

New (Graph Array); 

StepA := True; 
for i := 1 to 2 do 
begin 

for j := 1 to Step do 
begin 

Arrayindex := 1; 

if StepA then A := AMinGain + ADeltaStep*(j-1) 
else B := BMinGain + BDeltaStep*(j-1); 
if StepA then 
begin 

Mum2Str(f,A,Dss); 

Ds := ’A=’+ Dss; 
end 
else 
begin 

Num2Str(f,B,Dss); 

Ds := ’B=’+ Dss; 
end; 

pointno := 0; 

while pointno <= ( Points - 1) do 
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begin 

if StepA then B := NextGain2 
else A := NextGain2; 
for Term := InitDegree downto 0 do 
begin 

infixToPolish(lnfixArray[Term], polish): 
ComputePolish{polish,a,b,xlnitPoly[Term]); 
lnitPoly[Term].Re := xlnitPoly[Term]: 
InitPolyfTermJ.lm := 0.0; 
end; 

Degree := InitDegree; 

Poly := InitPoly: 

Laguerre(Degree, Poly, InitGuess, Tolerance, Maxiter, 
NumRoots, xAnswer, yAnswer, Iter, Error); 
Results(A, B, xAnswer, Error); 
pointno := pointno + 1 ; 
end; 

Arrayindex := Arrayindex - 1; 

PlotRootLocus2; 

end; 

StepA ;= False; 
end; 
end; 
begin 

end. { Unit MakeRoot } 
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unit RootsFinder(2000); 



{= Turbo Pascal Numerical Methods Toolbox = } 

{= Copyright (C) 1987 Borland International =} 

{= =} 
{= This unit provides procedures for finding the roots = } 

{= of a single equation in one real variable. = } 

(= = = = = = = = = = = = = = = = = = = = = = = = = = = = = _ = = } 



{$R+} { Enable range checking } 

interface 

uses 

MemTypes; 

const 

TNNearlyZero = IE-015; { Close to zero } 

TNArraySize =10; { Maximum size of vectors } 

type 

TNvector = array[O..TNArraySize] of Extended; 

TNIntVector = array[0. .TNArraySize] of integer; 

TNcomplex = record 

Re, Im : Extended; 
end; 

TNCompVector = array[0.. TNArraySize] of TNcomplex; 



procedure Laguerre(var Degree : integer; 
var Poly : TNCompVector; 
InitGuess : TNcomplex; 

Tol : Extended; 

Maxiter : integer; 
var NumRoots : integer; 
var Roots : TNCompVector; 
var yRoots : TNCompVector; 
var Iter : TNIntVector; 



var Error : byte); 

{= = = = _ = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = _=} 
{= =} 
{= Input: Degree, Poly, InitGuess, Tol, Maxiter = } 
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{= Output: Degree, Poly, NumRoots, Roots, yRoots, Iter, Error = } 

{= =} 

{= Purpose: This unit provides a procedure for finding all the = } 

{= roots (real and complex) to a polynomial. = } 

{= Laguerre's method with deflation is used. = } 

{= The user must input the coefficients of the quadratic = } 

{= and the tolerance in the answers generated. = } 

{= =} 

{= Pre-defined Types: TNcomplex = record =} 

{= Re, Im : Extended; = } 

{= end; =} 

{= =} 

{= TNIntVector = array [O..TNArraySize] of integer; = } 

{= TNCompVector = array[O..TNArraySize] of TNcomplex; = } 

{= =} 

{= Variables: Degree : integer; degree of deflated = } 

{= polynomial = } 

{= Poly : TNCompVector; coefficients of deflated = } 

{= polynomial where Poly[n] is = } 

{= the coefficient of X^'n = } 

{= InitGuess : TNcomplex; initial guess to a root = } 

{= (may be very crude) = } 

{= Tol : Extended; tolerance in the answer = } 

{= Maxiter : integer; number of iterations = } 

{= NumRoots : integer; number of roots calculated = } 

{= Roots : TNCompVector; the roots calculated = } 

{= yRoots : TNCompVector; the value of the function = } 

{= at the calculated roots = } 

{= Iter : TNIntVector; number iteration it took to = } 

{= find each root = } 

{= Error : byte; flags an error = } 

{= =} 

{= Errors: 0: No error = } 

{= 2: Degree <= 0 = } 

{= 3:Tol<=0 =} 

{= 4: Maxiter <0 = } 

{= =} 

= = = = 

implementation 

{====== ===============================} 

{= The following inline procedure and function are used to call the user = } 

{= defined procedures and functions pointed to by the ProcAddr parameter. = } 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = ==} 



81 



function UserFunction(X : Extended; ProcAddr : ProcPtr) : Extended; 
inline 

$205F, { MOVE.L (A7)+, AO } 

$4E90; { JSR (AO) } 

procedure UserProcedure(X : TNcomplex; var Y : TNcomplex; ProcAddr : ProcPtr); 
inline 

$205F, { MOVE.L (A7)+, AO } 

$4E90; { JSR (AO) } 



procedure Laguerre{(var Degree : integer; 
var Poly : TNCompVector; 
InitGuess : TNcomplex; 

Tol : Extended; 

Maxiter : integer; 
var NumRoots : integer; 
var Roots : TNCompVector; 
var y Roots : TNCompVector; 
var Iter : TNIntVector; 
var Error : byte)); 



type 

TNquadratic = record 

A, B, C : Extended; 
end; 



var 

Additer : integer; 

InitDegree : integer; 

InitPoly : TNCompVector; 

GuessRoot : TNcomplex; 

{ Here are a few complex operations } 

procedure Conjugate(var Cl, C2 : TNcomplex); 
begin 

C2.Re := C1.Re; 

C2.lm ;= -1.0 * Cl.lm; 
end; { procedure Conjugate } 

function Modulus(var Cl : TNcomplex) : Extended; 
begin 

Modulus := Sqrt(Sqr(C1.Re) + Sqr(Cl.lm)); 
end; { function Modulus } 
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procedure Add(var C1, C2, C3 : TNcompIex); 
begin 

03. Re 1 = 01 .Re + 02.Rej 
03.lm := 01 .Im + 02.1m; 
end; { procedure Add } 

procedure Sub(var 01, 02, 03 : TNcompIex): 
begin 

03.Re :=01.Re-02.Re; 

03.1m := 01.1m - 02.1m; 
end; { procedure Sub } 

procedure Mult(var 01, 02, 03 : TNcompIex); 
begin 

03.Re := 01 .Re * 02.Re - 01.1m * 02.1m; 

03.1m := 01.1m * 02.Re + 01. Re * 02.1m; 
end; { procedure Mult } 

procedure Divide(var 01, 02, 03 : TNcompIex); 
var 

Dumi, Dum2 : TNcompIex; 

E : Extended: 
begin 

0onjugate(02, Dumi); 

Mult(Ol, Dumi, Dum2): 

E := Sqr(Modulus(02)); 

03.Re := Dum2.Re / E; 

03.1m := Dum2.lm / E; 
end; { procedure Divide } 

procedure SquareRoot(var 01, 02 : TNcompIex); 
var 

R, Theta : Extended: 
begin 

R := Sqrt(Sqr(01.Re) + Sqr(Ol.lm)); 
if ABS(OI.Re) < TNNearlyZero then 
begin 

if 01.1m < 0 then 
Theta := Pi / 2 
else 

Theta := (-1.0 * Pi) / 2; 
end 
else 

if 01 .Re < 0 then 

Theta ;= ArcTan(Ol.lm / 01 .Re) + Pi 
else 

Theta := ArcTan(01.lm / 01 .Re); 



83 



C2.Re := Sqrt(R) * Cos(Theta / 2); 

C2.lm := Sqrt(R) * Sin(Theta / 2); 
end; { procedure SquareRoot } 

procedure lnitAndTest(var Degree : integer; 

var Poly : TNCompVector; 

Tol : Extended; 

Maxiter : integer; 

InitGuess rTNcomplex; 
var NumRoots : integer; 
var Roots ; TNCompVector; 
var yRoots : TNCompVector; 
var Iter : TNIntVector; 

var GuessRoot : TNcomplex; 
var InitDegree : integer; 
var InitPoly : TNCompVector; 
var Error : byte); 



{========================================} 

{= Input: Degree, Poly, Tol, Maxiter, InitGuess = } 

{= Output: InitDegree, InitPoly, Degree, Poly, NumRoots, = } 

{= Roots, yRoots, Iter, GuessRoot, Error = } 

{= =} 
{= This procedure sets the initial value of the above = } 

{= variables. This procedure also tests the tolerance = } 

{= (Tol), maximum number of iterations (Maxiter), and =} 

{= code. Finally, it examines the coefficients of Poly. =} 

{= If the constant term is zero, then zero is one of the = } 

{= roots and the polynomial is deflated accordingly. Also = } 

{= if the leading coefficient is zero, then Degree is = } 

{= reduced until the leading coefficient is non-zero. = } 

{========================================] 



var 

Term : integer; 
begin 

Error := 0; 
if Degree <= 0 then 

Error := 2; { degree is less than 2 } 

if Tol <= 0 then 
Error := 3; 
if Maxiter < 0 then 
Error := 4; 

if Error = 0 then 
begin 
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NumRoots := 0; 

GuessRoot := InitGuess; 

InitDegree := Degree; 

InitPoly := Poly; 

{ Reduce degree until leading coefficient <> zero } 
while (Degree > 0) and {Modulus(Poly[Degree)) < TNNearlyZero) do 
Degree := Pred(Degree); 

{ Deflate polynomial until the constant term <> zero } 
while (Modulus(Poly[0)) = 0) and (Degree > 0) do 
begin 

{ Zero is a root } 

NumRoots := Succ(Num Roots); 

Roots[NumRoots].Re := 0; 

Roots[NumRoots).lm := 0; 
yRoots[NumRoots].Re := 0; 
yRoots[NumRoots).lm := 0; 
lter[NumRoots) := 0; 

Degree := Pred(Degree); 
for Term := 0 to Degree do 

Poly[Term) := Poly[Term + 1]; 

end; 

end; 

end: { procedure InitAndTest } 

procedure FindOneRoot(Degree ; integer; 

Poly : TNCompVector; 

GuessRoot : TNcomplex; 

Tol : Extended: 

Maxiter : integer; 
var Root : TNcomplex; 
var yValue : TNcomplex; 
var Iter : integer; 
var Error : byte); 



{==== =============================} 

{= Input: Degree, Poly, GuessRoot, Tol, Maxiter =} 

{= Output: Root, y Value, Iter, Error =} 

{= =} 
{= This procedure approximates a single root of the polynomial = } 

{= Poly. The root must be approximated within Maxiter =} 

{= iterations to a tolerance of Tol. The root, value of the = } 

{= polynomial at the root (yValue), and the number of iterations = } 

{= (Iter) are returned. If no root is found, the appropriate error =} 

{= code (Error) is returned. = } 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 



var 
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Found : boolean; 

Dif : TNcomplex; 

yPrime, yDoublePrime : TNcomplex; 

procedure EvaluatePoly(Degree : integer; 

Poly : TNCompVector; 

X : TNcomplex; 

var yValue : TNcomplex; 
var yPrime : TNcomplex; 
var yDoublePrime : TNcomplex); 



{= Input: Degree, Poly, X = } 

{= Output: yValue, yPrime, yDoublePrime = } 

{= This procedure applies the technique of synthetic division to = } 

{= determine value (yValue), first derivative (yPrime) and second =} 

{= derivative (yDoublePrime) of the polynomial. Poly, at X. = } 

{= The 0th element of the first synthetic division is the = } 

{= value of Poly at X, the 1st element of the second synthetic = } 

{= division is the first derivative of Poly at X, and twice the = } 

{= 2nd element of the third synthetic division is the second = } 

{= derivative of Poly at X. = } 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = _ = = = = = = = =} 



var 

Loop : integer; 

Dummy, yDPdummy : TNcomplex; 

Deriv, Deriv2 : TNCompVector; 

begin 

Deriv[Degree] := PolyfDegree]; 
for Loop := Degree - 1 downto 0 do 
begin 

Mult(Deriv(Loop + 1], X, Dummy); 
Add(Dummy, Poly[Loop], Deriv[Loop]); 
end; 

yValue := Deriv[0]; { Value of Poly at X } 

Deriv2[Degree] := DerivfDegree]; 
for Loop := Degree - 1 downto 1 do 
begin 

Mult(Deriv2[Loop + 1], X, Dummy); 
Add(Dummy, Deriv[Loop], Deriv2[Loop]); 
end; 

yPrime ;= Deriv2[1]; { 1st deriv. of Poly at X } 
yDPdummy := Deriv2[Degree]; 
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for Loop .= Degree - 1 downto 2 do 
begin 

Mult(yDPdummy, X, Dummy); 

Add(Dummy, Deriv2[Loop], yDPdummy); 
end; 

yDoublePrime.Re := 2 * yDPdummy.Re; { 2nd derivative of Poly at X } 
yDoublePrime.lm := 2 * yDPdummy. Im; 
end; { procedure EvaluatePoly } 

procedure ConstructDifference(Degree : integer; 
y Value : TNcomplex; 

yPrime : TNcomplex; 

yDoublePrime : TNcomplex; 
var Dif : TNcomplex); 



{==== ===============================} 

{= Input: Degree, yValue, yPrime, yDoublePrime = } 

{= Output: Dif = } 

{= =} 

{= This procedure computes the difference between approximations: = } 

{= given information about the function and its first two = } 

{= derivatives. = } 



var 

yPrimeSQR, yTimesyDPrime, Sum, SRoot, 

Numert, Numer2, Numer, Denom : TNcomplex; 

begin 

Mult(yPrime, yPrime, yPrimeSQR); 
yPrimeSQR.Re := Sqr(Degree - 1) * yPrimeSQR.Re; 
yPrimeSQR. Im := Sqr(Degree - 1) * yPrimeSQR. Im; 

Mult(yValue, yDoublePrime, yTimesyDPrime); 
yTimesyDPrime. Re := (Degree - 1) * Degree * yTimesyDPrime. Re; 
yTimesyDPrime. Im := (Degree - 1) * Degree * yTimesyDPrime. Im; 
Sub(yPrimeSQR, yTimesyDPrime, Sum); 

SquareRoot(Sum, SRoot); 

Add(yPrime, SRoot, Numert); 

Sub(yPrime, SRoot, Numer2); 
if Modulus(Numerl) > Modulus(Numer2) then 
Numer ;= Numert 
else 

Numer := Numer2; 

Denom.Re := Degree * yValue.Re; 

Denom.lm := Degree * yValue.lm; 
if Modulus(Numer) < TNNearlyZero then 
begin 
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Dif.Re := 0; 

Dif.lm := 0; 
end 
else 

Divide(Denom, Numer, Di‘0: { The difference is the } 
{ inverse of the fractton } 
end; { procedure ConstructDifference } 

function TestForRoot(X, Dif, Y, Tol : Extended) : boolean; 



= = = = = = = = = = = = = = = = = = = = } 

{= These are the stopping criteria. Four different ones are = } 

{= provided. If you wish to change the active criteria, simply = } 

{= comment off the current criteria (including the appropriate OR) = } 

{= and remove the comment brackets from the criteria (including = } 

{= the appropriate OR) you wish to be active. = } 

{= = = = = = = = = = = = = = = = = = = = ^ = = = = = = = = = = =} 



begin 

TestForRoot := 

(ABS(Y) <= TNNearlyZero) 

or 

(ABS(Dif) < ABS(X * Tol)) 



(‘ or *) 

(* *) 

(* (ABS(DiO < Tol) *) 

(* *) 

(* or *) 

(* *) 
(* (ABS(Y) <= Tol) *) 



{ } 

{- Y=0 -} 

{- -} 

{- -} 

{- -} 

{- Relative change in X -} 

{- -} 

{- -} 

{- -} 

{■ -} 

{- Absolute change in X -} 

{- -} 

{- -} 

{- -} 

{- Absolute change in Y -} 

{ } 



{===============================.= 

{= The first criteria simply checks to see if the value of the 
{= function is zero. You should probably always keep this criteria 
{= active. 

{= 

{= The second criteria checks the relative error in X. This criteria 
{= evaluates the fractional change in X between interations. Note 
{= that X has been multiplied throught the inequality to avoid divide 
{= by zero errors. 

{= 



= } 
= } 
= } 
= } 
= } 
= } 
= } 
= } 
= } 
= } 
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{= The third criteria checks the absolute difference in X between = } 

{= iterations. = } 

{= =} 

{= The fourth criteria checks the absolute difference between = } 

{= the value of the function and zero. = } 



end; { procedure TestForRoot } 

begin { procedure FindOneRoot } 

Root := GuessRoot; 

Found ;= false; 

Iter := 0; 

EvaluatePoly(Degree, Poly, Root, y Value, yPrime, yDoublePrime): 

while (Iter < Maxiter) and not(Found) do 

begin 

Iter := Succ(lter); 

ConstructDifference(Degree, y Value, yPrime, yDoublePrime, Dif); 

Sub(Root, Dif, Root): 

EvaluatePoly(Degree, Poly, Root, yValue, yPrime, yDoublePrime); 

Found := TestForRoot(Modulus(Root), Modulus(Dif), Modulus(yValue), Tol); 
end; 

if not(Found) then Error := 1 ; { Iterations execeeded Maxiter } 

end; { procedure FindOneRoot } 

procedure ReducePoly(var Degree : integer; 
var Poly : TNCompVector; 

Root : TNcomplex); 



{========================================} 

{= Input: Degree, Poly, Root = } 

{= Output: Degree, Poly = } 

{= =} 

{= This procedure deflates the polynomial Poly by = } 

{= factoring out the Root. Degree is reduced by one. = } 



var 

Term : integer; 

NewPoly : TNCompVector; 

Dummy : TNcomplex; 

begin 

NewPoly[Degree - 1] := Poly[Degree]; 
for Term := Degree - 1 downto 1 do 
begin 

Mult(NewPoly[Term], Root, Dummy); 
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Add(Dummy, Poly[Term], NewPoly[Term - 1]); 
end; 

Degree := Pred(Degree); 

Poly := NewPoly: 
end; { procedure ReducePoly } 

begin { procedure Laguerre } 

lnitAndTest(Degree, Poly, Tol, Maxiter, InitGuess, NumRoots, Roots, 
yRools, Iter, GuessRool, InitDegree, InitPoly, Error); 
while (Degree > 0) and (Error = 0) do 
begin 

FindOneRoot(Degree, Poly, GuessRoot, Tol, Maxiter, 
Roots[NumRoots + 1), yRoots[NumRoots + 1], 
lter[Num Roots + 1], Error); 
if Error = 0 then 
begin 



{= = = = = = = = = = = = = = = = = = = = = = = = = =} 

{= The next statement refines the approximate root by =} 

{= plugging it into the original polynomial. This =} 

{= eliminates a lot of the round-off error =} 

{= accumulated through many iterations =} 

{= = = = = = = = = = = = = = = = = = = = = = = = = = } 



FindOneRoot(lnitDegree, InitPoly, Roots[N urn Roots -t- 1], 

Tol, Maxiter, Roots[Num Roots -h 1], 
yRoots[NumRoots + 1], Additer, Error); 
lter[NumRoots + 1] := lter[NumRoots + 1] + Additer; 

NumRoots := Succ(N urn Roots); 

ReducePoly(Degree, Poly, Roots[NumRoots]); { Reduce polynomial } 
end; 

GuessRoot := Roots[NumRoots]; 
end; 

end; { procedure Laguerre } 
begin 

end. { RootsOfEquat } 
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unit Message (6000); 



{= This unit provides the several messages which inform the = } 

{= user with some warnings and the help informations. =} 



{$T APPLFFTD} { Set the application type and creator } 

{$S+} 

{$!-} { Turn off I/O error checking } 

{$U RootsFinder} 

{$U SpecVar} 

{$U TurboGraph} 

interface 

uses 

MemTypes, QuickDraw, OSIntf, ToolIntf.PackIntf, 

PasPrinter.SANE.MacPrint, RootsFinder, SpecVar, 

{$S Second Segment} 

TurboGraph; 

procedure SetUpWindows; 
procedure MakeInfoScreen; 
procedure DoAbout; 
procedure AlertBoxI ; 
procedure AlertBox2; 
procedure AlertBoxS; 
procedure InfoGetEQParameter; 
procedure InfoGetCoeff; 
procedure InfoPlotOneParameter; 
procedure InfoPlotTwoParameter; 
procedure InfoPrint; 

implementation 

procedure SetUpWindows; 

{= Define all of the windows used in the program =} 
begin 

{ The initial information screen } 

DefineWindow(lnfoScrnWind, 100, 100, 362, 200, dBoxProc); 

{ The real transform window } 

DefineWindow(RootLocusWind, XMinGlb+5, YMinGlb+43, XMaxGlb-5, YMaxGlb-5, 
documentProc); 

DefineHeader(RootLocusWind, 'Root Locus'); 
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DefineWindow(RegendBox,100, 100,362,150, dBoxProc); 

{ The about box window } 

DefineWindow(AboutBoxWind, 125, 60, 387, 272, dBoxProc); 
DefineWindow(AlertBox,50, 100,360,170, dBoxProc); 

Define Window(HelpWind, 20, 60, 480,330, dBoxProc); 

end; { SetUpWindows } 



procedure MakeInfoScreen; 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 

{= Bring up a window with a description of what this program does =} 

{= = = = = = = = = = = = = = = = = =} 



var 

SavePort : GrafPtr; 

InLeft : integer; { Window offset from left } 

InTop : integer; { Window offset from top } 

begin 

GetPort(SavePort); { save the current grafport } 

SelectWind(lnfoScrnWind, FALSE); 
with Window[!nfoScrnWind].P''.portRect do 
begin 

InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 

InTop ;= (YMaxGlb - (Bottom-Top)) div 2; 
end; 

MoveWindow(Window[lnfoScrnWind].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(lnfoScrnWind, TRUE); 

TextFont(SystemFont); 

TextSize(14); 

TextStyle([]); 

MoveTo(20, 20); 

DrawStringC Welcome To'); 

TextSize(16); 

MoveTO(70, 50); 

DrawStringCMAC RootLocus'); 

TextSize(12); 

TextStyle([l); 

MoveTo(90, 80); 

DrawStringCHave a fun I!’); 
repeat until Button; 

HideWindow(Window[lnfoScrnWind].P); 

GetPort(SavePort); { save the current grafport } 
end; { MakeInfoScreen } 
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procedure DoAbout; 

{= Bring up the about box =} 
var 

SavePort : GrafPtr; 

InLeft : integer; { Window offset from left } 

InTop : integer; { Window offset from top } 

begin 

GetPort(SavePort); { save the current grafport } 

SelectWind(AboutBoxWind, FALSE); 

with Window[AboutBoxWind].P''.portRect do 

begin 

InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 

InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 

M9veWindow(Window[AboutBoxWind].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(AboutBoxWind, TRUE); 

T extFont(SystemFont) ; 

TextSize(18); 

TextStyle((]); 

MoveTo(60, 40); 

DrawStringCMAC RootLocus'); 

TextStyle([]); 

TextSize(12); 

MoveTo(85, 60); 

DrawStringCversion 1.00’); 

MoveTo(35, 90); 

DrawString(’Computer Aided Desgin Tool’); 

MoveTo(35, 110); 

DrawString('For RootLocus Analysis'); 

MoveTo(97, 160); 

DrawString('written by’); 

TextSize(13); 

MoveTo(85, 185); 

DrawStringCKO SUNG HOON’); 

T extFont(SystemFont) ; 
repeat until Button; 

HideWindow(Window[AboutBoxWind].P); 

SetPort(SavePort); 
end; { DoAbout } 

procedure AlertBoxI; 



(= = = = = = = = = = ^ = = = = = = = = = = = = = = = = = = = = = =} 

{= Bring up a window with a description of warning message. =} 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 
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var 

SavePort : GrafPtr; 

InLeft : integer; { Window offset from left } 

InTop ; integer; { Window offset from top } 

begin 

GetPort(SavePort); { save the current grafport } 

SelectWind(3, FALSE); 

with Window[3).P''.portRect do 

begin 

InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 
InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 

MoveWindow(Window[3J.P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(3, TRUE); 

T extFont(SystemFont) ; 

TextSize(12); 

TextStyle([]); 

MoveTo(30, 30); 

Drawstring (There is no Initial Degree.'); 
repeat until Button; 

HideWindow(Window[3].P); 

end; 



procedure AlertBox2; 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = _ = = = = = } 
{= Bring up a window with a description of warning message. =} 

= = = = = = = = = = = = = = = = = = = = = = = = = = } 

var 

SavePort : GrafPtr; 

InLeft : integer; { Window offset from left } 



InTop : integer; { Window offset from top } 
begin 

GetPort(SavePort); { save the current grafport } 

SelectWind(5, FALSE); 

with Window[5].P''.portRect do 

begin 

InLeft ;= (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 
InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 

MoveWindow(Window[5].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(5, TRUE); 

TextFont(SystemFont); 

TextSize(12); 

TextStyle([J); 

MoveTo(30, 30); 
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DrawStringCThere is no Initial Degree or’); 

MoveTo(50,50); 

DrawStringCCharacteristic Equation Coefficients.'); 
repeat until Button; 

HideWindow(Window[5].P); 

end; 

procedure AlerlBox3: 

{= Bring up a window with a description of warning message. =} 
var 

SavePort : GrafPtr; 

InLeft : integer; { Window offset from left } 

InTop : integer; { Window offset from top } 

begin 

GetPort(SavePort); { save the current grafport } 

SelectWind(5, FALSE); 

with Window[5].P''.portRect do 

begin 

InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 
InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 

MoveWindow(Window(5].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(5, TRUE); 

T extFont(SystemFont) ; 

TextSize(12); 

TextStyledl); 

MoveTo(30, 30); 

DrawString(’You have a wrong Input.'); 

MoveTo(50,50); 

DrawStringCCheck the Help Menu and try again.'); 
repeat until Button; 

HideWindow(Window[5].P); 

end; 

procedure InfoGetEQParameter; 



{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 
{= Bring up the Help information of EQ parameter dialog. =} 

= = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 



var 

SavePort : GrafPtr; 

InLeft : integer; { Window offset from left } 

InTop : integer; { Window offset from top } 

begin 

GetPort(SavePort); { save the current grafport } 
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SelectWind(HelpWind, FALSE); 

with Window[HelpWind].P''.portRect do 

begin 

InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 

InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 

MoveWindow(Window[HelpWind].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(HelpWind, TRUE); 

T extFont (System Font) ; 

TextSize(l4); 

TextStyie(n); 

MoveTo(150, 15); 

DrawString(’EQ PARAMETER'); 

TextStyle((]); 

TextSize(12); 

MoveTo(20, 40); 

DrawString(The EQ Parameter stands for equation parameter. It will allow you’); 
MoveTo(10, 55); 

Drawstring ('to input the degree of polynomial and some equation parameters. '); 
MoveTo(10,70); 

DrawString(’The degree of the polynomial should be 1 up to 10. Then there are ’); 
MoveTo(10, 85): 

DrawStringCthe default values for the other parameters. These values avoid '); 
MoveTo(10, 100); 

DrawStringCthe convergence error for almost all polynomials. But if '); 
MoveTo(10, 115): 

DrawStringCconvergence error messages appears on the screen, you can change ’); 
MoveTo(10, 130): 

DrawStringCthese parameter values. These parameters must satisfy the '); 
MoveTo(10, 145): 

DrawStringCfollowing conditions:'); 

MoveTo(30, 160); 

DrawStringC(l) Initial guess >= 0 ’); 

MoveTo(30, 175); 

DrawString(’(2) Maximum Iteration >= 0 ’); 

MoveTo(30, 190); 

DrawString('(3) Tolerance > 0 '); 

MoveTo(40, 260): 

DrawStringC" CLICK THE MOUSE ONCE TO RETURN THE MAIN MENU "’); 
repeat until Button; 

HideWindow(Window[HelpWind].P): 

SetPort(SavePort): 

end; 

procedure InfoGetCoeff; 



{ =============================} 

{= Bring up the Help information of Get Coeff dialog. =} 
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var 

Save Port : GrafPtr; 

InLeft : integer; { Window offset from left } 

InTop : integer; { Window offset from top } 

begin 

GetPort(SavePort); { save the current grafport } 

SelectWind(HelpWind, FALSE); 

with Window[HelpWind].P''.portRect do 

begin 

InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 

InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 

MoveWindow(Window[HelpWind].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(HelpWind, TRUE); 

T extFont(System Font) ; 

TextSize(14); 

TextStyle([]); 

MoveTo(170, 15); 

DrawStringCGET COEFF'); 

TextStyle([]); 

TextSize(12); 

MoveTo(20, 40); 

DrawString(The Get Coeff stands for Get Coefficients. It will allow you to'); 
MoveTo(10, 55); 

DrawString(’input the algebraic expression for the coefficients of the '); 
MoveTo(10,70); 

DrawString(’characteristic equation. It may have up to two undetermined ’); 
MoveTo(10, 85); 

DrawString('parameters(A and B). In case of the one parameter root locus '); 
MoveTo(10, 100); 

DrawString(’method, you use only one undetermined paramerter(A). The routine'); 
MoveTo(10, 115); 

DrawStringCuses standard algebraic, or infix, notation with parenthesis'); 
MoveTo(10, 130); 

DrawStringCallowed. Operater can include +, -, *, /, and ^(expomentiation).'); 
MoveTo(10, 145); 

DrawString('The unary minus sign is allowed.'); 

MoveTo(10, 160); 

DrawString('lf you choose the Get Coeff command without the degree of the'); 
MoveTo(10, 175); 

DrawString('polynomial,the message appears on the screen. It tells you that'); 
MoveTo(10, 190); 

DrawString(’the degree of the polynomial has not been entered yet.'); 

MoveTo(40, 260); 

DrawStringC" CLICK THE MOUSE ONCE TO RETURN THE MAIN MENU '"); 
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repeat until Button; 

Hide Window(Window[HelpWind].P) ; 
SetPort(SavePort); 
end; 

procedure InfoPlotOneParameter; 



{= = = = = = = = = = = = = = = = = = = = = = = =} 

{= Bring up the Help information of One parameter dialog. =} 

{ ^ = = = =} 



var 

SavePort : GrafPtr; 

InLeft : integer; { Window offset from left } 

InTop : integer; { Window offset from top } 

begin 

GetPort(SavePort); { save the current grafport } 

SelectWind(HelpWind. FALSE); 
with Window[HelpWind].P''.portRect do 
begin 

InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 

InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 

MoveWindow(Window[HelpWind].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibi!ity(HelpWind, TRUE); 

TextFont(SystemFont); 

TextSize(14); 

TextStyle([]); 

MoveTo(145, 15); 

DrawStringCONE PARAMETER'); 

TextStyle([]); 

TextSize(12); 

MoveTo(20, 35); 

DrawStringCOne parameter command will allow you to input the plot data of); 
MoveTo(10, 50); 

DrawStringCthe one parameter root locus method. The default values are shown’); 
MoveTo(10,65); 

DrawStringCin the dialog box. They can be changed as desired.'); 

MoveTo(20, 80); 

DrawStringC-. Enter the minimum and maximum gain values into the Min Gain'); 
MoveTo(10, 95); 

DrawStringCand the Max Gain insertion box.’); 

MoveTo(20, 110); 

DrawStringC-. Choose the type of interval. There are two types of interval:'); 
MoveTo(10, 125); 

DrawStringCLinear and Logarithmic interval. When you click the radio button,'); 
MoveTo(10, 140); 

DrawStringCthe desired type of interval is choosen.’); 
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MoveTo(20, 155); 

DrawStringC*. Choose the type of scale for the axis. There are two types of ’); 
MoveTo(10, 170): 

DrawString(’scale: Auto and Manual scale. When you click the radio button,’); 
MoveTo(10, 185); 

DrawString(’the desired type of scale is choosen.'); 

MoveTo(20, 200); 

DrawStringC-. Enter some number into Points to Plot inserton box in order to '); 
MoveTo(10, 215); 

DrawStringCdecide the plot resolution.’); 

MoveTo(40. 260); 

DrawStringC" CLICK THE MOUSE ONCE TO RETURN THE MAIN MENU "’); 
repeat until Button; 

HideWindow(Window[HelpWind].P); 

SetPort(SavePort); 

end; 



procedure InfoPlotTwoParameter; 

{====================================} 

{= Bring up the Help information of the two parameter dialog. =} 

{====================================} 

var 

Save Port : GrafPtr; 

InLeft : integer; { Window offset from left } 

InTop : integer: { Window offset from top } 



begin 

GetPort(SavePort); { save the current grafport } 

SelectWind(HelpWind, FALSE): 

with Window[HelpWind].P''.portRect do 

begin 

InLeft ;= (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 

InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 

MoveWindow(Window[HelpWind].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(HelpWind, TRUE); 

TextFont(SystemFont): 

TextSize(14); 

TextStyle([]): 

MoveTo(145, 15); 

DrawString(’TWO PARAMETER’); 

TextStyle([]): 

TextSize(12); 

MoveTo(20, 35); 

DrawString('Two parameter command will allow you to input the plot data of); 
MoveTo(10, 50); 

DrawString(’the two parameter root locus method. The default values are shown'); 
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MoveTo(10,65): 

DrawStringCin the dialog box. They can be changed as desired.'); 

MoveTo(20, 80); 

DrawString(The items in this dialog are similar to those in the one parameter'); 
MoveTo(10, 95); 

DrawStringCpIot data dialog box, but several items are different. ’); 

MoveTo(20, 110); 

DrawString(There exists one more undetemined parameter B to be inserted.The'); 
MoveTo(10, 125); 

DrawStringCHow many loci item lets you decide how many loci are to be drawn '); 
MoveTo{10, 140); 

DrawStringCfor each parameter. It will be 1 up to 10 for each parameter.’); 
MoveTo(10, 155); 

Drawstring ('There is no auto scale for axes. Only the manual scale is available'); 
MoveTo(10, 170); 

DrawString(’The last item is the marking and justification in order to draw'); 
MoveTo(10, 185); 

DrawStringCthe selected A and B values on the plot. Two buttons are chosen'); 
MoveTo(10, 200); 

DrawStringCeach time for each parameter; one for position, the other justifi- '); 
MoveTo(10, 215); 

DrawStringCcation to draw.'); 

MoveTo(40, 260); 

DrawStringC" CLICK THE MOUSE ONCE TO RETURN THE MAIN MENU ’"); 
repeal until Button; 

HideWindow(Window[HelpWind].P); 

SetPort(SavePort); 

end; 

procedure InfoPrint; 



{ = = = = = = = } 

{= Bring up the Help information of print out. = } 

= = = ^ = = = = = = = = = = = = _ = = = = = = 



var 

SavePort : GrafPtr; 

InLeft : integer; { Window offset from left } 

InTop : integer; { Window offset from top } 

begin 

GetPort(SavePort); { save the current grafport } 

SelectWind(HelpWind, FALSE); 
with Window[HelpWind].P''.portRect do 
begin 

InLeft := (XMaxGlb - (Right-Left)) div 2; { Calculate window offsets } 

InTop := (YMaxGlb - (Bottom-Top)) div 2; 
end; 

MoveWindow(Window[HelpWind].P, InLeft, InTop, TRUE); { Center the window } 
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SetVisibility(HelpWind, TRUE); 

T extFont(System Font) ; 

TextSize(14); 

TextStyle([]); 

MoveTo(160, 15); 

DrawStringCPRINT OUT'); 

TextStyle([]): 

TextSize(12): 

MoveTo(20, 35); 

DrawString('ln order to print out your work, there are a few ways available.'); 
MoveTo(20, 50); 

DrawStringC-.Select the print command in the File menu. This allows the ’); 
MoveTo(10,65); 

DrawStringCuser to get a hard copy of any plot displayed.'); 

MoveTo(20, 80); 

DrawStringC-. Hold the command key and then type the number 4 to print '); 
MoveTo(10, 95); 

DrawStringCthe contents of the active window immediately.'); 

MoveTo(20, 110); 

DrawStringC-. Press the command and shift key and type the number 3 to '); 
MoveTo(10, 125); 

DrawStringCcreate a MacPaint document.'); 

MoveTo(40, 260); 

DrawStringC" CLICK THE MOUSE ONCE TO RETURN THE MAIN MENU '"); 
repeat until Button; 

HideWindow(Window[HelpWind].P); 

SetPort(SavePort); 

end; 

begin 

end.{unit message} 
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unit MyDialog (7000); 



(= = = = = = = = = = = = = = = = = _ = = = = = = _ = = = = = =) 
{= This unit supports four dialog boxes for input data. =} 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 



{$B+} { Set the bundle bit } 

{$R MacRootLocus.Rsrc} { Identify resource file for menu and icon info } 
{$T APPLFFTD} { Set the application type and creator } 

{$S+} { generate segmented code} 

{$!-} { Turn off I/O error checking } 

{$U SpecVar} 

{$U RootsFinder} 

{$U MakeRoot} 

{$U TurboGraph} 

{$U Message} 

interface 

uses 

MemTypes, QuickDraw, OSIntf, ToolIntf.PackIntf, 
PasPrinter,SANE,MacPrint, RootsFinder, SpecVar, 

{$S Second Segment} 

TurboGraph, 

Message, 

{$S Third Segment} 

MakeRoot: 

procedure MakeRegend; 
procedure GetEQParameter; 
procedure GetCoeff; 
procedure PlotOneParameter; 
procedure PLotTwoParameter; 

implementation 



procedure MakeRegend: 

{= = = = = = = = = = = = = = = = = = = = = =} 

{= This procedure provides the informaton box with =} 
{= the capabilety of the text editor =} 

{===============================} 



var 

SavePort : GrafPtr; 

In Left : integer; { Window offset from left } 

InTop : integer; { Window offset from top } 

begin 

GetPort(SavePort): { save the current grafport } 
SelectWind(3, FALSE); 
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Textsize(?); 

with Window[3].P''.portRect do 
begin 

InLeft ;= Round((XMaxGlb - (Right-Left)) / 1.38); { Calculate window offsets ) 
InTop := Round((YMaxGlb - (Bottom-Top)) /1.05); 
end; 

MoveWindow(Window[3].P, InLeft, InTop, TRUE); { Center the window } 
SetVisibility(3, TRUE); 
txRect ;= thePort''.portRect; 
textH := TENew(txRect, txRect); 

TEIdle(textH): 

TextInputEnabled := True; 
end; 



procedure GetEQParameter; 

(= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 
{= This procedure provides the dialog box for the = ) 

{= equation parameter of the characteristic equation. =} 

{= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 



const 

CancelBtn = 1 ; 

OKBtn = 2; 

InitDegreeText = 3; 

InitGuessText = 4; 

MaxIterText = 5; 

ToleranceText = 6; 
var 

tempdegree,templter : Longint; 
tempGuess : TNcomplex; 
tempTolerance : Extended; 

BoxPlot : Boolean; 
begin 

theDialog := GetNewDialog(256, nil, pointer(-l)); 
GetDltem(theDialog, 3, theType, h3, r); 
GetDltem(theDialog, 4, theType, h4, r); 
GetDltem(theDialog, 5, theType, h5, r); 
GetDltem(theDialog, 6, theType, h6, r); 

Done := False; 
tempGuess. Re := 1 .0; 
tempGuess.lm := 0.0; 
tempiter := 100; 
tempTolerance := IE-7; 

BoxPlot := False; 
repeat 

ModalDialog(nil, itemHit); 
case itemHit of 
CancelBtn : done := True; 
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OKBtn 



; begin 

InitDegree := lempdegree; 

InitGuess.Re := tempGuess.Re; 

InitGuess.lm := tempGuess.lm; 

Maxiter := tempiter; 

Tolerance := tempTolerance; 
done := True; 
end; 

InitDegreeText : begin 

GetlText(h3, s); 

StringToNum(s, tempdegree); 

InitDegreeStatus := True; 
if (tempdegree <1 ) or (tempdegree >10) then 
begin 

sysbeep(IO): 

end; 

end; 

InitGuessText : begin 

GetlText(h4, s); 
tempGuess.Re := Str2Num(s); 
if (tempGuess.Re <0 ) then 
begin 

sysbeep(IO); 

end; 

end; 

MaxIterText ; begin 

GetlText(h5, s); 

StringToNum(s, tempiter); 
if (tempiter <0 ) then 
begin 

sysbeep(IO); 

end; 

end; 

ToleranceText : begin 

GetlText(h6, s); 
tempTolerance := Str2Num(s); 
if (tempTolerance <= 0 ) then 
begin 

sysbeep(IO); 

end; 

end; 

end; 

until done; 

DisposDialog(theDialog); 

end; 

procedure GetCoeff; 
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{= This procedure provides the dialog box for the 
{= coefficients of the characteristic equation. 



=} 

=} 



{= = = 

type char_set = set of char; 
const 






CancelBtn = 1 ; 
OKBtn = 2; 

SOText =3; 

SUext =4; 

S2Text = 5; 

S3Text = 6; 

S4Text = 7; 

SSText = 8; 

S6Text = 9; 

STText = 10; 

SSText =11; 

S9Text = 12; 



SIOText = 13; 



var 

Idno.tempitemhit.sssindex.ssslength 
: Integer; 

h ; Array[1..13] of Handle; 
set_valid :char_set; 
op_set : char_set; 
sss: str255; 
sssc: char; 

noerrorcheck.leftparen: boolean; 

Procedure skipblank ; 
begin 

while ((sss[sssindex] = ’ ’) and (sssindex <= ssslength)) do 
sssindex := sssindex + 1 ; 



Procedure Checknext ( setofnext1,setofnext2:char_set); 

var opchar: char; 

begin 

opchar ;= sssc; 
sssindex := sssindex +1 ; 
sssc := sss[sssindex]; 
if (sssc = ’ ’) then begin 
skipblank; 

sssc := sss[sssindex]; 
if sssindex <= ssslength then begin 
if not (sssc in setofnextl) then 



end; 
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noerrorcheck:= False; 
end else 

if (opchar in (op_set + ('('])) then 
noerrorcheck := false 

end else 

if not (sssc in setofnext2) then 
noerrorcheck:= false 

end: 



Procedure CheckExpression(matchparen:boolean); 
begin 

leftparen:= true; 

while ((noerrorcheck and leftparen) and (sssindex < ssslength)) do 
case (sssc) of 

'O'. .'9': Checknext(op_set+['.',’)'],op_set+['0'.. '9', I 
'a'.'b': Checknext(op_set+ [')'], op_set+['a','b', ’)']); 

Checknext(['0'..’9‘.'a','b','n,['0‘..'9','a','b‘,'(']); 

Checknext(['0'..'9'],['0'..'9']): 

'(':begin 

Checknext([‘0’..’9’.’a’,'b','-','('],[’0’..’9’, 'a', 'b','-’, ’(']): 
CheckExpression (false); 
if sssc = ')' then 

if sssindex < ssslength then begin 
sssc := sssisssindexj; 

if sssfsssindex + 1] = ’ ’ then skipblank; 
if sssindex < ssslength then 

Checknext(op_set+[')'].op_set+[')']) 
else if not matchparen then noerrorcheck ;= false 
end else begin 

if not matchparen then noerrorcheck := false 
end 

else noerrorcheck := false; 
end: 



')’:begin 

if matchparen then noerrorcheck;= false; 
leftparen := false: 
end; 

' ';begin 

skipblank ; 

if sssindex = ssslength then begin 

if not (sss[sssindex] in ['0’..'9','a',’b']) then 
noerrorcheck:= false 

end else if sssindex < ssslength then sssc := sss[sssindex]; 
end 



end; 



if not noerrorcheck then sysbeep (30); 
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end: 

{begin 

if not InitDegreeStatus then 
begin 

SysBeep(30): 

AlertBoxI; 

end 

else} 

begin 

set_valid := ['O’. .'9', ’a’, 'b', 7','-’,’ '.''''.T,')','.’); 
op.set := 

if not InitDegreeStatus then 
begin 

SysBeep(30): 

AlertBoxI; 

end 

else 

begin 

Idno := 300 + InitDegree; 

theDialog := GetNewDialog(ldno, nil, Pointer (-1)); 
for n ;= 3 to InitDegree + 3 do 
GetDltem(theDialog, n, theType, h|n], r); 
done := False; 

ModalDialog(nil, itemHit); 



repeat 

tempitemhit := itemHit; 
case itemHit of 
CancelBtn : done := True; 

OKBtn : done := True; 

3..12 ; begin 

while tempitemhit = itemHit do 
ModalDialog(nil, tempitemhit); 

GetIText (h(itemHit],sss); 
noerrorcheck ;= true; 
sssindex;= 1 ; 
ssslength := Length (sss); 
sssc := sss[1]: 

if not (sss[1] in ['O'. .'9', 'a', 'b', then 
noerrorcheck := False: 

CheckExpression ( True); 

GetlText(h[itemHit],lnfixArray[lnitDegree + 3 - itemHit]) 
GetCoeffStatus :=True; 
itemhit := tempitemhit: 
if itemhit = OKBtn then done := true 
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end; 



end; 

until done; 

DisposDialog(theDialog) ; 
end; 
end; 

{end;} 



procedure PlotOneParameter; 

= = = = = ^ = = = = = = = = = = = = = = = = = = = = = = = = 

{= This procedure provides the dialog box for the =} 

{= plot data of the one parameter root locus method. =} 



const 

CancelBtn = 1 ; 

OKBtn = 2; 

LinearBtn = 3; 

LogarithmicBtn = 4; 

AutoBtn = 5; 

ManualBtn = 6; 

AMinGainText = 7; 

AMaxGainText = 8; 

XMinText = 9; 

XMaxText =10; 

YMinText =11; 

YMaxText = 12; 

PlotPointsText = 13; 

var 

saveSoundVol : Integer; 
radButton : array [3.. 6] of ControlHandle; 
h : Array[7..13] of Handle; 
tempAMinGain : Extended; 
tempAMaxGain : Extended; 
tempXMin : Extended; 

tempXMax : Extended; 

tempYMin : Extended; 

tempYMax : Extended; 

tempPoints : Longint; 

OkPlot :boolean; 
begin 

if not InitDegreeStatus or not GetCoeffStatus then 
begin 

SysBeep(30); 

Alert Box2; 
end 
else 
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begin 

GetSoundVol (saveSoundVol); 

SetSoundVol (1): 

theDialog := GetNewDialog(257, nil, Pointer (-1)); 
for n := LinearBtn to ManualBtn do 

GetDltem(theDialog,n,theType,Handle(radButton[n]),r): 
SetCtIValue (radButton[LinearBtn],1 ); 
linear := True; 

SetCtIValue (radButton[AutoBtn],1); 

AutoScale := True; 
for n := 7 to 13 do 

GetDltem(theDialog, n, theType, h[n], r); 
tempAMinGain := Str2Num(’0.T); 
tempAMaxGain := Str2Num(’10000'); 
lempXMin ;= Str2Num{'-10'); 
tempXMax := Str2Num('5’); 
tempYMin := Str2Num('-10'); 
tempYMax := Str2Num('10'); 

StringToNum(’50’, tempPoints); 
done := False; 

OkPlot := False; 
repeat 

ModalDialog(nil,itemHit); 
case itemHit of 
CancelBtn : done := True; 

OKBtn : begin 

done := True; 

OkPlot := True; 
end; 

LinearBtn : begin 

SetSoundVol (1); 

for n := LinearBtn to LogarithmicBtn do 

SetCtlValue(radButton[n],Ord(n = itemHit)); 
linear := True; 
end; 

LogarithmicBtn ; begin 

SetSoundVol (1); 

for n := LinearBtn to LogarithmicBtn do 

SetCtlValue(radButton[n],Ord(n = itemHit)); 
linear ;= False; 
end; 

AutoBtn : begin 

SetSoundVol (1); 

for n := AutoBtn to ManualBtn do 

SetCtlValue(radButton[n],Ord(n = itemHit)); 
AutoScale := True; 
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end; 



ManualBtn : begin 

SetSoundVol (1); 

for n := AutoBtn to ManualBtn do 

SetCtlValue(radButton[n),Ord(n = itemHit)); 
AutoScale := False: 
end; 

AMinGainText : begin 

GetlText(h[itemHit] , s); 
tempAMinGain := Str2Num(s); 
if (tempAMinGain < 0 ) or (tempAMinGain > 1 e7) then 
begin 

sysbeep(lO); 

end; 

end; 

AMaxGainText : begin 

GetlText(h[itemHit] , s); 
tempAMaxGain := Str2Num(s); 
if (tempAMaxGain < 0 ) or (tempAMaxGain > 1e7) then 
begin 

sysbeep(IO); 

end; 

end; 

XMinText ; begin 

GetlText(h[itemHit] , s); 
tempXMin := Str2Num(s); 
if (tempXMin <-1e7 ) or (tempXMin > 1e7) then 
begin 

sysbeep(IO); 

end; 

end; 

XMaxText : begin 

GetlText(h[itemHit) , s); 
tempXMax := Str2Num(s); 
if (tempXMax <-1e7 ) or (tempXMax > 1e7) then 
begin 

sysbeep(IO): 

end: 

end: 

YMinText : begin 

GetlText(h[itemHit] , s); 
tempYMin := Str2Num(s); 
if (tempYMin <-1e7 ) or (tempYMin > 1e7) then 
begin 

sysbeep(IO): 

end; 
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YMaxText 



end; 

: begin 

GetlText(h[itemHit] , s); 
tempYMax := Str2Num(s); 
if (tempYMax <-1 e7 ) or (tempYMax > 1 e7) then 
begin 

sysbeep(IO): 
end; 
end; 

PlotPointsText : begin 

GetlText(h[itemHit] , s); 

StringToNum(s, tempPoints); 
if (tempPoints <= 0) or (tempPoints > 150) then 
begin 

sysbeep(IO): 

end: 

end; 

end;{case end} 
until done; 

DisposDialog(theDialog) ; 

if OkPlot then 

begin 

AMinGain := tempAMinGain; 

AMaxGain := tempAMaxGain; 

XMn := tempXMin; 

XMx := tempXMax; 

YMn := tempYMin; 

YMx := tempYMax; 

Points := tempPoints; 

HideCursor; 

Rewrite(OutFile,'RootLocus1 .data’); 

GetRootl ; 

Close(OutFile): 

ShowCursor; 

MakeRegend; 

end;{if} 

end;{case} 

end;{PlotOneParameter} 
procedure PLotTwoParameter; 



{= = = = = = = = = = = = = = = = = = = = = = = = = = = = =} 

{= This procedure provides the dialog box for the =} 

{= plot data of the hvo parameter root locus method. =} 

{= = = = = = = = = = = = = = = = = = = = = = = } 



const 

CancelBtn = 1 ; 
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OKBtn = 2; 

LinearBtn = 3; 

LogarithmicBtn = 4; 

AMarkStBtn = 5; 

AMarkEdBtn = 6; 

AJustifyRBtn = 7; 

AJustifyLBtn = 8; 

BMarkStBtn = 9; 

BMarkEdBtn = 10; 

BJustifyRBln =11; 

BJuslifyLStn =12; 

AMinGainText = 13; 

AMaxGainText =14; 

BMinGainText = 15; 

BManGainText =16; 

LociText =17; 

PlotPointsText =18; 

XMinText = 19; 

XMaxTexl = 20; 

YMinText = 21 ; 

YMaxText = 22; 

var 

saveSoundVol.i : Inleger; 

radButlon : array [3. .12] of ControlHandle; 

h : Array[13..22j of Handle; 

tempAMinGain : Extended; 

tempAMaxGain : Extended; 

tempBMinGain ; Extended; 

tempBMaxGain ; Extended; 

tempXMin : Extended; 

tempXMax : Extended; 

tempYMin : Extended; 

tempYMax : Extended; 

tempStep, tempPoints : Longint; 

OkPlot : Boolean; 
begin 

if not InitDegreeStatus or not GetCoeffStatus then 
begin 

SysBeep{30); 

AlertBox2; 

end 

else 

begin 

GetSoundVol (saveSoundVol); 

SetSoundVol (1); 

theDialog := GetNewDialog(258, nil, Pointer (-1)); 
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for n := LinearBtn to BJustifyLBtn do 
GetDltem(theDialog,n,theType,Handle(radButton[n]),r); 
SetCtIValue (radButton[LinearBtn],1 ); 
linear := True; 

SetCtIValue (radButton[AMarkStBtn],1); 

AMarkStatus := False; 

SetCtIValue (radButton[AJustifyRBtn],1); 

ARJustification ;= True; 

SetCtIValue (radButton[BMarkStBtn],1); 

BMarkStatus ;= False; 

SetCtIValue (radButton[BJustifyRBtn],1); 

BRJustification := True; 
for n := 13 to 22 do 

GetDltem(theDialog, n, theType, h[n], r); 
tempAMinGain ;= Str2Num(’0.1'); 
tempAMaxGain := Str2Num{’10000’); 
tempBMinGain := Str2Num('0.1'); 
tempBMaxGain := Str2Num('10000'); 
tempXMin := Str2Num('-10'); 
tempXMax ;= Str2Num('5'); 
tempYMin := Str2NumC-10'); 
tempYMax := Str2Num('10'); 

StringToNum('4', tempStep); 

StringToNum('50', tempPoints); 
done := False; 

OkPlot := False; 
repeat 

ModalDialog(nil.itemHit): 
case item Hit of 
CancelBtn : done := True; 

OKBfn : begin 

done := True; 

OkPlot := True; 
end; 

LinearBtn : begin 

SetSoundVol (1); 

for n := LinearBtn to LogarithmicBtn do 

SetCtlValue(radButton[n],Ord(n = itemHit)); 
linear := True; 
end; 

LogarithmicBtn : begin 

SetSoundVol (1); 

for n := LinearBtn to LogarithmicBtn do 

SetCtlValue(radButton[n],Ord(n = itemHit)); 
linear ;= False; 
end; 
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AMarkStBtn : begin 

SetSoundVol (1); 

for n := AMarkStBtn to AMarkEdBtn do 

SetCtlValue{radButton(n],Ord(n = itemHit)); 
AMarkStatus := False; 
end; 

AMarkEdBtn : begin 

SetSoundVol (1); 

for n := AMarkStBtn to AMarkEdBtn do 

SetCliVa!ue(rad8utton[n],Ord{n = itemHit)); 
AMarkStatus := True; 
end; 

AJustifyRBtn : begin 

SetSoundVol (1); 

for n ;= AJustifyRBtn to AJustifyLBtn do 
SetCtlValue{radButton[n],Ord(n = itemHit)); 
ARJustification := True; 
end; 

AJustifyLBtn : begin 

SetSoundVol (1); 

for n := AJustifyRBtn to AJustifyLBtn do 
SetCtiVaJu8{radButton[n],Ord{n = itemHit)); 
ARJustification := False; 
end; 

BMarkStBtn : begin 

SetSoundVol (1); 

for n := BMarkStBtn to BMarkEdBtn do 

SetCtlValue{radButton[n].Ord(n = itemHit)); 
BMarkStatus := False; 
end; 

BMarkEdBtn : begin 

SetSoundVol (1); 

for n := BMarkStBtn to BMarkEdBtn do 

SetCl!Value(radButton[n],Ord{n = itemHit)); 
BMarkStatus := True; 
end; 

BJustifyRBtn ; begin 

SetSoundVol (1); 

for n := BJustifyRBtn to BJustifyLBtn do 
SetCtlValue(radButton(n],Ord(n = itemHit)); 
BR Justification ;= True; 
end; 

BJustifyLBtn : begin 
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SetSoundVol (1 ); 

for n := BJustifyRBtn to BJustifyLBtn do 
SetCtlValue(radButton[n],Ord(n = itemHit)); 
BRJustification := False; 
end; 

AMinGainText : begin 

GetlText(h[itemHit],s); 
tempAMinGain := Str2Num(s); 
if (tempAMinGain < 0 ) or (tempAMinGain > 1e7) then 
begin 

sysbeep(IO); 

end; 

end; 

AMaxGainText : begin 

GetlText(h[itemHit],s); 
tempAMaxGain := Str2Num(s); 

if (tempAMaxGain < 0 ) or (tempAMaxGain > 1 e7) then 
begin 

sysbeep(IO); 

end; 

end; 

BMinGainText : begin 

GetlText(h[itemHit],s); 
tempBMinGain := Str2Num(s); 

if (tempBMinGain < 0 ) or (tempBMinGain > 1e7) then 
begin 

sysbeep(IO); 

end; 

end; 

BManGainText : begin 

GetlText(h[itemHit],s); 
tempBMaxGain := Str2Num(s); 

if (tempBMaxGain < 0 ) or (tempBMaxGain > 1e7) then 
begin 

sysbeep(IO); 

end; 

end; 

XMinText : begin 

GetlText(h[itemHit],s); 
tempXMin := Str2Num(s); 
if (tempXMin <-1e7 ) or (tempXMin > 1e7) then 
begin 

sysbeep(IO); 

end; 

end; 

XMaxText : begin 
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GellText(h[ilemHitl,s); 
tempXMax := Str2Num(s); 
if (lempXMax <-1e7 ) or (lempXMax > 1e7) then 
begin 

sysbeep(IO): 

end; 

end; 

YMinTexl : begin 

GetIT ext(h[itemHlt).s); 
tempYMin > Str2Num(s); 
if (tempYMin <-1e7 ) or (lempYMin > 1e7) then 
begin 

sysbeep(IO); 

end; 

end; 

YMaxText : begin 

GetlText(h[ilemHit),s); 
tempYMax := Sfr2Num(s); 
if (tempYMax <-1e7 ) or (tempYMax > 1e7) then 
begin 

sysbeep(IO); 

end; 

end; 

LociText : begin 

GetIText(h[itemHit).s); 

StringToNum(s, tempStep); 
if (tempStep <= 0) or (tempStep > 15) then 
begin 

sysbeep(IO); 

end; 

end; 

PlotPointsText : begin 

GetlText(h(itemHit),s); 

StringToNum(s, tempPoints); 
if (tempPoints <= 0) or (tempPoints > 150) then 
begin 

sysbeep(IO); 

end; 

end; 

end;{case end} 
until done; 

DisposDialog(theDialog); 

if OkPlot then 

begin 

AMinGain := tempAMinGain; 

AMaxGain := tempAMaxGain; 

BMinGain := tempBMinGain; 
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BMaxGain ;= tempBMaxGain; 

XMn := tempXMin; 

XMx := tempXMax; 

YMn := tempYMin; 

YMx := tempYMax; 

Step := tempStep; 

Points := tempPoints; 

HideCursor; 

Rewrite(OutFile ,'RootLocus2 .data') ; 
GetRoot2; 

Close(OutFile); 

ShowCursor; 

MakeRegend; 

end;{if} 

end;{case} 

end;{PlotTwoParameter} 

begin 

end.{unit MyDialog} 
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unit TurboGraph(1 1000); 



{= = = = = = = = = = = = = = = = = = = = = = = } 

{= =} 

{= Turbo Pascal Numerical Methods Toolbox = } 

{= Copyright (C) 1987 Borland International =} 

{= =} 

{= This unit provides routines for displaying graphics. = } 

{= =} 

(==========================================} 



{$U RootsFinder} 

{$U SpecVar} 

interface 

uses 

MemTypes, QuickDraw, OSIntf, Toolintf, PackIntf.PasPrinter, SANE, MacPrint , 

RootsFinder, SpecVar; 

const 

MaxWorldsGIb = 10; { The maximum number of worlds that can be defined } 

MaxWindowsGIb = 10; { The maximum number of windows that can be defined } 

{MaxPlotGIb - 1500;}{ The maximum number of points in a plot array } 

XI Offset = 40; { The upper left corner of the axis } 

Y1 Offset =10; 

X20ffset = 30; { The lower right corner of the axis } 

Y20ffset = 90; 

type 

WrkString = Str255; { A general string type } 

WorldType = record { Used to store world coordinates } 

XI, Y1,X2, Y2 :real; 
end; 

WindowType = record { Used to store window information } 

XI, Y1, X2, Y2 : integer; { The windows screen coordinates } 

Header : WrkString; { Header for a window } 

W : WindowRecord; { Mac window record } 

P : WindowPtr; { Mac window pointer } 

H : PicHandle; { A handle to a picture } 
end; 

Worlds = array[1.. MaxWorldsGIb] of WorldType; { Holds world info } 
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Windows = array[1 ..MaxWindowsGIb] of WindowType;{ Holds window info } 
{ PlotArray = array[1..MaxPlotGlb, 1..2] of Extended;} 

var 

{ Coordinates of the currently active window } 

XMinGIb, YMinGIb, XMaxGlb, YMaxGlb : integer; 

{ Coordinates of the currently selected world } 

XIWIdGIb. X2WldGlb, YIWIdGIb, Y2WldGlb : real; 

{ World coordinate scaling factors } 

AxGlb, AyGIb, BxGlb, ByGIb : real; 

{ Coordinates of an axis if one is defined } 

XIRefGIb, X2RefGlb, YIRefGIb, Y2RefGlb : integer; 

{ The maximum world and window number defined } 

MaxWorldGIb, MaxWindowGIb : integer; 

{ The currently selected world and window } 

WorldNdxGlb, WindowNdxGlb : integer; 

{ Aspect ratio for a true circle } 

AspectFactor : real; 

AspectGIb : real; 

{ Flags if an axis is defined } 

AxisGIb : boolean; 

{ Flags if hatching is turned on } 

HatchGIb : boolean; 

{ Flags if clipping should be performed } 

ClippingGIb : boolean; 

{ The currently selected line style } 

LineStyleGIb : integer; 

{ The current foreground color } 

ForeColorGIb : integer; 

{ Holds the worlds defined by the user } 

World ; Worlds; 

{ Holds the windows defined by the user } 

Window : Windows; 
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procedure Error{ProcName : WrkString); 

{ Report that an error occurred in the procedure named "ProcName" } 

procedure SetForegroundColor(Color : integer); 

{ Set the foreground drawing color } 

procedure SetBackgroundColor(Color : integer): 

{ Set the background color } 

procedure DP(X, Y : integer); 

{ Plot a pixel at position (X, Y) } 

function PD(X, Y : integer) : boolean; 

{ Return true if the color of the pixel at (X, Y) matches ForeColorGIb } 

procedure DrawStraight(X1 , X2, Y : integer); 

{ Draw a horizontal line from X1,Y to X2,Y } 

procedure lnvertWindow(WinNum : integer); 

{ Invert the window referenced by WinNum } 

function RealToStr(R : real) : WrkString; 

{ Returns the string representation of the real R. } 

function lnlToStr(l : integer) : WrkString: 

{ Returns the string representation of the integer I. } 

procedure SetLineStyle(Style ; integer); 

{ Select the current line style } 

function GetLineStyle : integer; 

{ Returns the current line style } 

procedure DefineWorld(WorldNum : integer; X_1, Y_1, X_2, Y_2 : real); 

{ Defines a world coordinate system with a specific number } 

procedure SelectWorld(WorldNum : integer); 

{ Select the world associated with WorldNum } 

procedure DefineWindow(WinNum, XLo, YLo, XHi, YHi, WindowType : integer); 
{ Defines a window with a specific window number and window type } 

procedure SetVisibility(WinNum ; integer; Visible : boolean): 

{ Sets the visibility of a window to either TRUE or FALSE } 

procedure SelectWind(WinNum : integer; Visible : boolean); 

{ Selects a window as visible or invisible } 
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procedure DefineHeader(WinNum : integer; Hdr : WrkString); 

{ Define a header for a window } 

procedure RemoveHeader(WinNum : integer); 

{ Clears a header from a window } 

procedure ReDefineWindow{WinNum, X1, Y1, X2, Y2, WindowType : integer); 
{ Redefines the coordinates of an existing window } 

function WindowX(X : real) : integer; 

{ Converts an X world coordinate into a screen coordinate } 

function WindowY(Y : real) : integer; 

{ Converts a Y world coordinate into a screen coordinate } 

procedure ResetWorlds; 

{ Resets all worlds to the maximum dimensions of the screen } 

procedure InitGraphic; 

{ Initializes the graphics system } 

function Clip(var X1, Y1, X2, Y2 : integer) : boolean; 

{ Clips a line to the current window and returns TRUE if any part } 

{ of the line is still in the window after the clip operation } 

procedure DrawPointfXr, Yr : real); 

{ Draw a point in world coordinates } 

function PointDrawn(Xr, Yr : real) : boolean; 

{ Returns TRUE if the point at (Xr, Yr) is drawn } 

procedure DrawLine(X1, Y1, X2, Y2 : real); 

{ Draw a line in world coordinates } 

procedure DrawLineClipped (X1, Y1, X2, Y2 : integer); 

{ Draw a line clipped in screen coordinates } 

procedure DrawCrossDiag(X, Y, Scale ; integer); 

{ Draw a cross at screen coordinate (X, Y) with a scaling factor } 

procedure DrawWye(X, Y, Scale ; integer); 

{ Draw a Y at screen coordinate (X, Y) with a scaling factor } 

procedure DrawDiamond(X, Y, Scale : integer); 

{ Draw a cross at screen coordinate (X, Y) with a scaling factor } 
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procedure DrawCircleDirect(Xr, Yr, R : integer); 

{ Draw a circle at screen coordinate (Xr, Yr) with a radius R } 

procedure DrawCircle(X_R, Y_R, Xradius : real); 

{ Draw a circle at world coordinate (X_R, Y_R) with a radius R } 

procedure DrawCross(X1 , Y1, Scale : integer): 

{ Draw a cross at screen coordinate (X, Y) with a scaling factor } 

procedure DrawStar(X, Y, Scale : integer); 

{ Draw a star at screen coordinate (X, Y) with a scaling factor } 

procedure DrawSquareC(X1 , Y1, X2, Y2 : integer; Fill : boolean): 

{ Draw a square in screen coordinates with optional filling } 

procedure DrawSquare(X1, Y1, X2, Y2 : real; Fill : boolean); 

{ Draw a square in world coordinates with optional filling } 

procedure DrawAscii(X, Y : integer; Size, CharByte ; byte); 

{ Draw a character with ASCII code CharByte } 

procedure DrawText(X, Y, Scale : integer; Txt : Str255); 

{ Draw a string at screen coordinate (X, Y) with a scaling factor } 

procedure DrawTextW(X, Y : real; Scale : integer; Txt : WrkString); 

{ Draw a string at world coordinate (X, Y) with a scaling factor } 

procedure TextStyle(Face : Style): 

{ Face = (bold, italic, underline, outline, shadow, condense, extend) } 
procedure HardCopy(TopWin : boolean): 

{ Do a screen dump of the currently selected window (TopWin = TRUE) } 
{ or the whole screen (TopWin = FALSE) } 

procedure OpenPic(WinNum : integer; ShowPic : boolean): 

{ Open a picture for a specific window and only show the drawing ) 

{ if ShowPic is set to TRUE } 

procedure DrawPic(WinNum : integer); 

{ Draw the picture associated with a window } 

procedure ErasePic(WinNum : integer); 

{ Erase the picture associated with a window } 

procedure ClearWindow(WinNum : integer): 

{ Clear the content portion of a window } 
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function Vv'hereX : integer; 

{ Returns the X cursor position } 

function WhereY : integer; 

{ Returns the Y cursor position } 

procedure SetWindow(X1 , Y1, X2, Y2 : integer); 

{ Creates an invisible windo\w inset from the currently selected one } 

procedure FindWorld(l : integer; A : PlotArray; NPoints : integer); 

{ Finds a world to fit a polygon defined in A } 

procedure FindWorld1(l : integer; XMn,YMn,XMx,YMx:ExTended ); 

{ Finds a world to fit a polygon defined in A } 

procedure DrawAxis(Footer1, Footer2 : WrkString; Arrows : boolean); 
{ Draws an axis with Footers and optional arrows on the axis } 

procedure ResetAxis; 

{ Sets AxisGIb to TRUE } 

procedure DrawPolygon(A : PlotArray; First, NPoints, Line, Scale, 
Lines : integer; CrossHairs : boolean); 

{ Draws a polygon defined in A with "NPoints" points, line style "Line” } 
{ and optional Lines from the axis to the points and cross hairs. } 

procedure Hatch(X_1, Y_1, X_2, Y_2, Delta : real); 

{ Hatch a bar in a histogram } 

procedure DrawHistogram(A :PlotArray; NPoints : integer; 

Hatching : boolean; HatchStyle : integer); 

{ Draw a histogram defined in A with "NPoints" points and optional } 

{ hatching with a given hatch style } 

implementation 

procedure Error{(ProcName : WrkString)}; 

{ Report that an error occurred in the procedure named "ProcName" } 
begin 

DrawStringCERROR in ’); 

DrawString(ProcName); 

DrawStringC Press the Button to exit.’): 
repeat until Button; 

Halt; 

end; { Error } 
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procedure SetForegroundColor{(Color : integer)}; 

{ Set the foreground drawing color } 
begin 

case Color of 

0 : begin 

ForeColor(BlackColor) ; 

ForeColorGIb := 0; 
end; 

1 ; begin 

ForeColor(WhiteColor) ; 

ForeColorGIb := 1 ; 
end; 
end; 

end; { SetForegroundColor } 

procedure SetBackgroundColor{(Color : integer)}; 

{ Set the background color } 
begin 

case Color of 

0 : BackColor(BlackColor); 

1 : BackColor(WhiteColor); 
end; 

end; { SetBackgroundColor } 

procedure DP{(X, Y ; integer)}; 

{ Plot a pixel at position (X, Y) } 
begin 

MoveTo(X, Y); 

LineTo(X, Y); 
end; { DP } 

function PD{(X, Y : integer) : boolean}; 

{ Return true if the color of the pixel at (X, Y) matches ForeColorGIb } 
var 

BlackPixel ; boolean; 

PixelColor : integer; 
begin 

BlackPixel ;= GetPixel(X, Y); 
if BlackPixel then 
PixelColor := 0 
else 

PixelColor := 1 ; 

PD := PixelColor = ForeColorGIb; 
end; { PD } 

procedure DrawStraight{(X1 , X2, Y : integer)}; 

{ Draw a horizontal line from X1 ,Y to X2,Y } 
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begin 

MoveTo(X1, Y): 

LineTo(X2, Y); 
end; { DrawStraight } 

procedure lnvertWindow{(WinNum : integer)}; 

{ Invert the vtrindow referenced by WinNum } 
begin 

if WinNum in [1..MaxWindowGlb] then 

lnvertRect(Window[WinNum].W.Port.PortRect) 

else 

Error(’lnvert Window’); 
end; { InvertWindow } 

function RealToStr{(R : real) : WrkString}; 

{ Returns the string representation of the real R. } 
var 

Int, Frac : Longint; 

S1, S2 ; Str255; 

Negitive : boolean; 
begin 

51 := ”; 

52 := 

if R < 0.0 then 
Negitive := TRUE 
else 

Negitive := FALSE; 

R ;= ABS(R); 

Int := Trunc(R); 

Frac := Round(100.0 * (R - Int)); 
NumToString(lnt, S1); 

NumToString(Frac, S2); 
if Length(S2) = 1 then 
S2 := S2 + 'O’: 

S2 := S1 + + S2; 

if Negitive then 
RealToStr := + S2 

else 

RealToStr := ” + S2; 
end; { RealToStr } 

function lntToStr{(l : integer) : WrkString); 

{ Returns the string representation of the integer I. } 
var 

Form : DecForm; 

Str : DecStr; 
begin 
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Form .Style '.= FixedDecimal; 

Form. Digits := 0; 

Num2Str(Form, I, Str); 

IntToStr := Str; 
end; { IntToStr } 

procedure SetLineStyle{(Style : integer)}; 
var 

LineStyle : Pattern; 
begin 

case Style of 

0 : LineStyle := Black; 

1 : LineStyle := White; 

2 : LineStyle := Gray; 

3 : LineStyle := LtGray; 

4 : LineStyle ;= DkGray; 
otherwise 

Style := 0; 

LineStyle := Black; 
end; 

LineStyleGIb := Style; 

PenPat(LineStyle); 
end; { SetLineStyle } 

function GetLineStyle{ : integer}; 
begin 

GetLinestyle ;= LineStyleGIb; 
end; { GetLineStyle } 

procedure DefineWorld{(WorldNum : integer; X_1, Y_1, X_2, Y_2 ; real)}; 
begin 

if ((X_1 <> X_2) and (Y_1 <> Y_2)) and 
(WorldNum in [1 ..MaxWorldsGIbj) then 
with World[WorldNum] do 
begin 

X1 := X_1 ; 

Y1 := Y_1; 

X2 := X_2; 

Y2 := Y_2; 

if WorldNum > MaxWorldGIb then 
MaxWorldGIb := WorldNum; 
end 

else if WorldNum in [1..MaxWorldsGlb] then 
Error(’DefineWorld #1’) 
else 

Error('DefineWorld #2'); 
end; { DefineWorld } 
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procedure SelectWorld{(WorldNum : integer)}; 
begin 

if (WorldNum in [1..MaxWorldGlb]) then 
with WorldfWorldNum] do 
begin 

WorldNdxGlb := WorldNum; 

XIWIdGIb := X1; 

YIWIdGIb := Y1; 

X2WldGlb := X2; 

Y2WldGlb ;= Y2; 
end 
else 

Error(’SelectWorld’): 
end; { SelectWorld } 

procedure DefineWindow{(WinNum, XLo, YLo, XHi, YHi, WindowType : integer; Visible 

boolean)}; 

var 

BoundsRect : Rect; 

Title : Str255; 

RefCon : Longint; 

Visible ; boolean; 

GoAwayFlag : boolean; 

begin 

if WinNum in [1 ..MaxWindowsGIb] then 
begin 

if WinNum > MaxWindowGIb then 
MaxWindowGIb := WinNum; 
with BoundsRect do 
begin 

Left := XLo; 

Top := YLo; 

Right := XHi; 

Bottom := YHi; 
end; 

Window[WinNum].X1 := XLo; 

Window[WinNum].Y1 := YLo; 

Window[WinNum].X2 := XHi; 

Window[WinNum].Y2 ;= YHi; 

Title := 

GoAwayFlag := TRUE; 

Visible := FALSE; 

RefCon := WinNum; 
if WindowType = documentProc then 
WindowType := 16 * WindowType + 8; { Add Zoom box } 
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Window[WinNum].P := NewWindow(@Window[WinNum].W, BoundsRect, 
Title, Visible, WindowType, 

POINTER(-I), GoAwayFlag, RefCon); 
end 
else 

Error(’DefineWindow'): 
end; { DefineWindow } 

procedure SetVisibility{(WinNum : integer; Visible : boolean)}; 
begin 

ShowHide(Window[WinNum).P, Visible); 
end; { SetVisibility } 

procedure SelectWind{(WinNum : integer; Visible : boolean)}; 
begin 

if (WinNum in [1..MaxWindowGlb]) then 
begin 

SelectWindow(Window[WinNum].P}; 
if Visible then 

ShowWindow(Window[WinNum].P); 
SetPort(@Window[WinNum].W.Port); 
with Window[WinNum] do 
begin 

WindowNdxGlb := WinNum; 

XIRefGib := W.Port.PortRect.Left; 

YIRefGib := W.Porf.PortRect.Top; 

X2RefGlb := W.Port.PortRect.Right; 

Y2RefGlb ;= W.Port.PortRect. Bottom; 

BxGlb := (X2RefGlb - XIRefGib - 16) / (X2WldGlb - XIWIdGIb); 
ByGIb := (Y2RefGlb - YIRefGib - 16) / (Y2WldGlb - YIWIdGIb); 
AxGlb := XIRefGib - XIWIdGIb * BxGlb; 

AyGIb := YIRefGib - YIWIdGIb * ByGIb; 

AxisGIb := FALSE; 
end; 
end 
else 

Error('SelectWind’) ; 
end; { SelectWind } 

procedure DefineHeader{(WinNum : integer; Hdr : WrkString)}; 
begin 

Window[WinNum]. Header := Hdr; 

SetWTitle(Window[WinNum].P, Hdr); 
end; { DefineHeader } 

procedure RemoveHeader{(WinNum : integer)}; 
begin 
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DefineHeader(WinNum, "); 
end; { RemoveHeader } 

procedure ReDefineWindow{(WinNum, X1, Y1, X2, Y2, WindowType : integer)}; 
begin 

if (WinNum in [1 ..MaxWindowsGIb]) then 
begin 

CloseWindow(Window[WinNum].P); 

Window[WinNum].P := NIL; 

DefineWindow(WinNum, XI, Y1, X2, Y2, WindowType): 
SelectWind(WinNum, FALSE); 

DefineHeader(WinNum, Window[WinNum]. Header); 
end 
else 

Error('ReDefineWindow'); 
end; { ReDefineWindow } 

function WindowX{(X : real) : integer}; 
var 

Temp : real; 
begin 

Temp := AxGlb + BxGlb * X; 
if Temp > Maxint then 
WindowX := Maxint 
else if Temp < -32767 then 
WindowX ;= -32767 
else 

WindowX := trunc(Temp); 
end; { WindowX } 

function WindowY{(Y : real) : integer}; 
var 

Temp : real; 
begin 

Temp := Ay Gib + ByGIb * Y; 
if Temp > Maxint then 
WindowY := Maxint 
else if Temp < -32767 then 
WindowY := -32767 
else 

WindowY := trunc(Temp); 
end; { WindowY } 

procedure ResetWorlds; 
var 

I : integer: 
begin 
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for I := 1 to MaxWorldsGIb do 
DefineWorld(l, XMinGIb, YMinGIb, XMaxGlb, YMaxGlb); 
SelectWorld(l): 
end; { ResetWorlds } 

procedure InitGraphic; 
var 

Index : integer; 

BoundsRect : Rect; 

Title : Str255; 

Visible : boolean; 

RefCon : Longint; 

GoAwayFlag : boolean; 

WindowType : integer; 
begin 

XMinGIb := screenBits.bounds.Left; 

YMinGIb := screenBits.bounds.Top; 

XMaxGlb := screenBits.bounds. Right; 

YMaxGlb := screenBits.bounds. Bottom; 
for Index := 1 to MaxWindowsGIb do 
begin 

Window[lndex].P := NIL; 

Window(!ndex].H := NIL; 
end; 

ResetWorlds; 

MaxWorldGIb := 0; 

MaxWindowGIb := 0; 

WindowNdxGlb := 0; 

WorldNdxGlb := 0; 

AspectFactor := 0.44; 

AspectGIb := ABS(AspectFactor) * AspectFactor; 

AxisGIb := false; 

HatchGIb := false; 

ClippingGIb := true; 

SetLineStyle(O); { Solid Black lines } 
end; { InitGraphic } 

function Ciip{(var X1, Y1, X2, Y2 : integer) : boolean}; 
var 

1x1, iyl, 1x2, Iy2, Dummy, XILoc, X2Loc : integer; 
ClipLoc : boolean; 

Temp : real; 

function lnside(X, Xxl, Xx2 : integer) : integer; 
begin 

Inside := 0; 
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if X < Xx1 then 
Inside := -1 
else if X > Xx2 then 
Inside ;= 1 ; 
end; { Inside } 

begin { Clip } 

Clip := true; 

ClipLoc := true; 
if ClippingGIb then 
begin 

XI Loc := X1 RefGIb; 

X2LOC := X2RefGlb; 

1x1 ;= lnside(X1, XI Loc, X2Loc); 
lyl := lnside(Y1, Y1 RefGIb, Y2RefGlb); 

1x2 := lnside(X2, XI Loc, X2Loc); 

Iy2 := lnside(Y2, Y1 RefGIb, Y2RefGlb); 
if (1x1 or 1x2 or lyl or Iy2) <> 0 then 
begin 

if XI <> X2 then 
begin 

if 1x1 <>0 then 
begin 

if 1x1 < 0 then 
Dummy := XI Loc 
else 

Dummy := X2Loc; 
if Y2 <> Y1 then 
begin 

Temp := (Y2 - Y1) / (X2 - XI) * (Dummy - XI); 
if Temp > Maxint then 
Temp := Maxint 
else if Temp < -32767 then 
Temp := -32767; 

Y1 := Y1 + trunc(Temp); 
end; 

XI := Dummy; 
end; 

if (1x2 <> 0) and (XI <> X2) then 
begin 

if 1x2 < 0 then 
Dummy := XI Loc 
else 

Dummy := X2Loc; 
if Y2 <> Y1 then 
begin 

Temp := (Y2 - Y1) / (X2 - XI) * (Dummy - XI); 
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if Temp > Maxint then 
Temp := Maxint 
else if Temp < -32767 then 
Temp := -32767; 

Y2 := Y1 + trunc(Temp): 
end; 

X2 := Dummy; 
end; 

Iy1 := lnside(Y1, YIRefGIb, Y2RefGlb); 

Iy2 := !nside(Y2, YIRefGIb, Y2RefGlb); 
end; 

if Y1 <> Y2 then 
begin 

if Iy1 <> 0 then 
begin 

if Iy1 < 0 then 
Dummy ;= YIRefGIb 
else 

Dummy := Y2RefGlb; 
if X1 <> X2 then 
begin 

Temp ;= (X2 - X1) / (Y2 - Y1) * (Dummy - Y1); 
if Temp > Maxint then 
Temp := Maxint 
else if Temp < -32767 then 
Temp ;= -32767; 

X1 := X1 + trunc(Temp); 
end; 

Y1 := Dummy; 
end; 

if Iy2 <> 0 then 
begin 

if Iy2 < 0 then 
Dummy := YIRefGIb 
else 

Dummy := Y2RefGlb; 
if X1 <> X2 then 
begin 

Temp (X2 - X1) / (Y2 - Y1) * (Dummy - Y1); 
if Temp > Maxint then 
Temp := Maxint 
else if Temp < -32767 then 
Temp := -32767; 

X2 := X1 + trunc(Temp); 
end; 

Y2 := Dummy; 
end; 
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end: 

lyl := lnside(Yl , Y1 RefGIb, Y2RefGlb); 

Iy2 := lnside(Y2, Y1 RefGIb, Y2RefGlb); 
if (Iy1 <> 0) or (Iy2 <> 0) then 
ClipLoc := false; 
if ClipLoc then 
begin 

Ixl := lnside(X1, XILoc, X2Loc); 

1x2 := lnside(X2, XILoc, X2Loc); 
if (1x2 <> 0) or (1x1 <> 0) then 
ClipLoc := false; 
end: 

Clip := ClipLoc: 
end; 
end; 

end; { Clip } 

procedure DrawPoint{(Xr, Yr : real)}; 
var 

X, Y : integer: 
begin 

X := WindowX(Xr): 

Y := WindowY(Yr): 

DP(X, Y); 
end; { DrawPoint } 

function PointDrawn{(Xr, Yr : real) : boolean}; 
begin 

PointDrawn := PD(WindowX(Xr), WindowY(Yr)); 
end; { PointDrawn } 

procedure DrawLine{(X1, Y1, X2, Y2 : real)}; 
begin 

MoveTo(WindowX(X1), WindowY(YI)); 
LineTo(WindowX(X2), WindowY(Y2)); 
end; { DrawLine } 

procedure DrawLineClipped{(X1 , Y1, X2, Y2 : integer)}; 
begin 

if Clip(X1, Y1, X2, Y2) then 
begin 

MoveTo(Xl, Y1); 

LineTo(X2, Y2); 

{ end 
else 
begin 

Moveto(X1 - 3,Y1 - 3); 
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TextSize(9); 

DrawString(Ds): 

LabelSet := False;} 
end; 

end; { DrawLineClipped } 

procedure DrawCrossDiag{(X, Y, Scale : integer)}; 
begin 

DrawLineClipped(X - Scale, Y + Scale, X + Scale + 1, Y - Scale - 1); 
DrawLineClipped(X - Scale, Y - Scale, X + Scale + 1, Y + Scale + 1); 
end; ( OrawCrossDiag } 

procedure DrawWye{(X, Y, Scale : integer)}; 
begin 

DrawLineClipped(X - Scale, Y - Scale, X, Y); 

DrawLineClipped(X + Scale, Y - Scale, X, Y); 

DrawLineClipped(X, Y, X, Y + Scale); 
end; { DrawWye } 

procedure DrawDiamond{(X, Y, Scale : integer)}; 
begin 

DrawLineClipped(X - Scale, Y, X, Y - Scale - 1); 

DrawLineClipped(X, Y - Scale + 1, X + Scale, Y + 1); 
DrawUneC!ipp8d(X + Scale, Y + 1, X, Y + Scale); 

DrawLineCiipped(X, Y + Scale, X - Scale, Y); 
end; { DrawDiamond } 

procedure DrawCircleDirect{(Xr, Yr, R : integer)}; 
type 

Circ = array[1..14] of integer; 
var 

Xkl , Xk2, Ykl , Yk2, Xpl , Yp1 , Xp2, Yp2 : integer; 

Xfact, Yfact : real; 

I : integer; 

X : Circ; 



procedure initX; 
begin 



X[1] 


:= 0; 


X[2] 


:= 121 


X[3] 


:= 239 


X[4] 


:= 355 


X[5] 


:= 465 


X[6] 


:= 568 


X[7] 


:= 663 
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X[8] := 749; 

X(9] := 823: 

X(10] := 885; 

X[11] := 935; 

X[12] := 971; 

X[13] := 993; 

X[14] := 1000; 
end; { InitX } 

begin { DrawCircleDirect } 

InitX; 

Xfact := abs(R); 

Yfact ;= Xfact * AspectGIb; 
if Xfact > 0.0 then 
begin 

Xk1 ;= trunc(X[1] * Xfact + 0.5); 

Ykl := trunc(X[14] * Yfact + 0.5); 
for I := 2 to 14 do 
begin 

Xk2 := trunc(X(l] * Xfact + 0.5); 

Yk2 := trunc(X[14 - I + 1] * Yfact + 0.5); 

Xpl := Xr - Xkl ; 

Ypl := Yr + Ykl ; 

Xp2 := Xr - Xk2; 

Yp2 ;= Yr + Yk2; 

DrawLine(Xp1 , Ypl, Xp2, Yp2); 

Xpl := Xr + Xkl ; 

Xp2 := Xr + Xk2; 

DrawLine(Xp1, Ypl, Xp2, Yp2); 

Ypl := Yr - Ykl ; 

Yp2 := Yr - Yk2; 

DrawLine(Xp1, Ypl + 1, Xp2, Yp2 + 1); 

Xpl := Xr - Xkl ; 

Xp2 := Xr - Xk2; 

DrawLine(Xp1 , Ypl + 1, Xp2, Yp2 + 1); 

Xkl := Xk2; 

Ykl := Yk2; 
end; 
end 
else 

DP(Xr, Yr): 

end; { DrawCircleDirect } 

procedure DrawCircle{(X_R, Y_R, Xradius : real)}; 
begin 

DrawCircleDirect(WindowX(X_R), WindowY(Y_R), trunc(Xradius)): 
end; { DrawCircle } 
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procedure DrawCross{(X1, Y1, Scale : integer)}; 
begin 

DrawLineClipped(X1 - Scale, Y1, X1 + Scale + 2, Y1); 
DrawLineClipped(X1 , Y1 - Scale, X1, Y1 "+ Scale + 1); 
end; { DrawCross } 

procedure DrawStar{(X, Y, Scale : integer)}; 
begin 

DrawUneClipped(X - Scale, Y + Scale, X + Scale + 1, Y - Scale - 1); 
DrawLlneC!ipped(X - Scale, Y - Scale, X + Scale + 1, Y + Scale + 1); 
DrawLineCiipped{X - Scale - 2, Y, X + Scale + 4, Y); 
end; { DrawStar } 

procedure DrawSquareC{(X1 , Y1, X2, Y2 : integer; Fill : boolean)}; 
var 

I ; integer; 

procedure DSC(X1, X2, Y : integer); 
begin 

DrawStraight(X1 , X2, Y); 
end; { DSC } 

begin { DrawSquareC } 
if not Fill then 
begin 

DrawLineClipped(X1, Y1, X2, Y1); 

DrawLineClipped{X2, Y1, X2, Y2); 

DrawLineClipped{X1, Y2, X2, Y2); 

DrawLineClipped{X1, Y2, XI, Y1); 
end 
else 

for I := Y2 to Y1 do 
DSC(X1, X2, I); 
end; { DrawSquareC } 

procedure DrawSquare{(X1 , Y1, X2, Y2 : real; Fill ; boolean)}; 
var 

I, X1Loc, Y1Loc, X2Loc, Y2Loc : integer; 

DirectModeLoc : boolean; 

procedure DS(X1, X2, Y : integer); 
begin 

if LineStyleGIb = 0 then 
DrawStraight(X1, X2, Y) 
else 

DrawLine(X1, Y, X2, Y); 
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end; { DS } 



procedure DSC(X1, X2, Y : integer); 
begin 

DS(X1, X2, Y); 
end; { DSC } 

procedure DrawSqr(X1, Y1, X2, Y2 : integer; Fill : boolean); 
var 

I : integer; 
begin 

if not Fill then 
begin 

DS(X1, X2, Y1); 

DrawLine(X2, Y1, X2, Y2); 

DS(X1, X2, Y2); 

DrawLine(X1, Y2, X1, Y1); 
end 
else 

for I := Y1 to Y2 do 
DS(X1, X2, I); 
end; { DrawSqr } 

begin { DrawSquare } 

XIL 0 C := WindowX(XI); 

YIL 0 C := WindowY(YI); 

X2LOC := WindowX(X2); 

Y2LOC := WindowY(Y2); 
if not Fill then 
begin 

DSC(X1 Loc, X2Loc, Y1 Loc); 

DrawLineClipped(X2Loc, YILoc, X2Loc, Y2Loc); 
DSC(X1Loc, X2LOC, Y2Loc); 

DrawLineClipped(X1Loc, Y2Loc, XIL 0 C, YILoc); 
end 
else 

for I := Y1 Loc to Y2Loc do 
DSC(X1Loc, X2LOC, I); 
end; { DrawSquare } 

procedure DrawAscii{(X, Y : integer; Size, CharByte : byte)}; 
begin 

MoveTo(X, Y); 

TextSize(Size * 12); 

DrawChar(Chr(CharByte)); 
end; { DrawAscii } 
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procedure DrawText{(X, Y, Scale ; integer; Txt : WrkString)}; 
var 

Index : integer; 

EscStr : boolean; 

StringLen : integer; 

AsciiValue ; integer; 

SymbolScale ; integer; 

SymbolCode : integer; 
begin 

Index := 1 ; 

EscStr := FALSE; 

StringLen := Length(Txt); 

while (Index <= StringLen) and (not EscStr) do 

begin 

if Txt[lndex] = #27 then 
EscStr ;= TRUE; 

Index := Index + 1 ; 
end; 

if not EscStr then 
begin 

MoveTo(X, Y); 

TextSi 2 e(Scale * 12); 

DrawString(Txt); 

end 

else 

begin 

Index := 1; 

while Index <= StringLen do 
begin 

AsciiValue := Ord(Txt[lndex]); 
if AsciiValue = 27 then 
begin 

SymbolScale := Scale; 

Index := Index + 1 ; 
if Index <= StringLen then 
begin 

SymbolCode := Ord{Txt[lndex]) - 48; 

if (Index + 2 <= StringLen) and (Ord(Txt[lndex + 1]) = 64) then 
begin 

SymbolCode := Ord(Txt[lndex]) - 48; 

Index := Index + 2; 
end; 

case SymbolCode of 

1 : DrawCross(X + SymbolScale, Y + Scale, SymbolScale); 

2 : DrawCrossDiag(X + SymbolScale, Y + Scale, SymbolScale); 
3,4 : DrawSquareC(X, Y + (SymbolScale shl 1 ) - 1 , 

X + (SymbolScale shl 1), Y - 1, (SymbolCode = 4)); 
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5 : begin 

DrawDiamond(X + trunc(1.5 * SymbolScale), 

Y + SymbolScale - 1, SymbolScale + 1); 

X := X + SymbolScale; 
end; 

6 : DrawWye(X + SymbolScale, Y + SymbolScale - 1, SymbolScale); 

7 : begin 

DrawStar(X + SymbolScale shl 1, Y + SymbolScale - 1, SymbolScale); 
X := X + SymbolScale shl 1 ; 
end; 

8 : DrawCircleDirect(X + SymbolScale, Y + (SymbolScale shr 1), 

SymbolScale + 1 ); 

end; 

X := X + 3 * SymbolScale; 

SymbolScale := Scale; 
end; 
end 
else 

DrawAscii(X, Y, Scale, Ascii Value); 

Index := Index + 1 ; 
end; 
end; 

end; { DrawText } 

procedure DrawTextW{(X, Y : real; Scale : integer; Txt : WrkString)}; 
begin 

DrawText(WindowX(X), WindowY(Y), Scale, Txt); 
end; { DrawTextW } 

procedure TextStyle{(Face : Style)}; 

{ Face = (bold, italic, underline, outline, shadow, condense, extend) } 
begin 

TextFace(Face); 
end; { TextStyle } 

procedure HardCopy{(TopWin : boolean)}; 
begin 

PrDrvrOpen; 
if TopWin then 
{ Print the top folder. } 

PrCtlCall(iPrEvtCtl, LPrEvtTop, 0, LScreenBits) 
else 

{ Print the whole screen. } 

PrCtlCall(iPrEvtCtl, LPrEvtAII, 0, LScreenBits); 

PrDrvrClose; 
end; { HardCopy } 
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procedure OpenPic{(WinNum : integer; ShowPic : boolean)}; 
begin 

if Window[WinNum].H <> NIL then 
begin 

KillPicture(Window[WinNum].H); 

Window[WinNum].H := NIL; 
end; 

RectRgn(Window(WinNum).W.Port.clipRgn, ScreenBits.bounds); 
Window[WinNum].H := OpenPicture(Window[WinNum].W.Port.PortRect); 
if ShowPic then 
ShowPen 
end; { OpenPic } 

procedure DrawPic{(WinNum : integer)}; 
var 

PictRect : Rect; 
begin 

PictRect := Window[WinNum].W.Port.PortRect; 
with PictRect do 
begin 

if ((Bottom - Top) < 200) OR ((Right - Left) < 200) then 
begin 

Right := Right - 16; { so we don’t overwrite the grow region } 

Bottom := Bottom - 16;{ on a small window. } 

end; 
end; 

DrawPictu re(Window[WinNum]. H, PictRect) ; 
end; { DrawPic } 

procedure ErasePic{(WinNum : integer)}; 
begin 

if Window[WinNum].H <> NIL then 
begin 

KillPicture(Window[WinNum].H); 

Window[WinNum].H := NIL; 
end; 

end; { ErasePic } 

procedure ClearWindow{(WinNum : integer)}; 
begin 

EraseRect(Window[WinN urn} .W. Port. PortRect); 
end; { ClearWindow } 

function WhereX{ : integer}; 
var 

Pt : Point; 
begin 
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GetPen(Pt): 

WhereX ;= Pt.H; 
end; { WhereX } 

function WhereY{ ; integer}; 
var 

Pt : Point; 
begin 

GetPen(Pt); 

WhereY := Pt.V; 
end; { WhereY } 

procedure SetWindow{(Xl , Y1, X2, Y2 : integer)}; 
begin 

X1 RefGIb := X1 ; 

Y1 RefGIb := Y1 ; 

X2RefGlb := X2; 

Y2RefGlb := Y2; 

BxGlb := (X2 - X1) / (X2WldGlb - XIWIdGIb); 

ByGIb := (Y2 - Y1) / (Y2WldGlb - YIWIdGIb); 

AxGlb := X1 - XIWIdGIb * BxGlb; 

AyGIb := Y1 - YIWIdGIb * ByGIb; 

AxisGIb := FALSE; 
end; { SetWindow } 

procedure FindWorld{(l : integer; A : PlotArray; NPoints : integer)}; 
var 

J : integer; 

Xmax, Ymax, Xmin, Ymin, Xmid, Ymid, Xdiff, Ydiff : real; 
begin 

NPoints := abs(NPoints); 
if NPoints > 2 then 
if I in [1..MaxWorldsGlb] then 
begin 

Xmax := A[1, 1]; 

Ymax := A[1, 2]; 

Xmin := Xmax; 

Ymin := Ymax; 

for J := 2 to NPoints do 

begin 

if A[J, 1] > Xmax then 
Xmax ;= A(J, 1] 
else 

if A[J, 1] < Xmin then 
Xmin := A[J, 1j; 
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if A[J, 2] > Ymax then 
Ymax := A(J, 2] 
else 

if A(J, 2] < Ymin then 
Ymin := A[J, 2]; 

end; 

Xmin := Round(Xmin): 

Ymin := Round(Ymin) - 0.5; 

Xmax ;= Round(Xmax); 

Ymax Round(Ymax) + 0.5; 

DefineWor!d(l, Xmin, Ymin, Xmax, Ymax); 

SelectWorld(l); 

end 

else 

Error(’FindWorld #1’) 
else 

Error(’FindWorld # 2’); 
end; { FindWorld } 

procedure FindWorld1{(l : integer; A : PlotArray; NPoints : integer)}; 
var 

J : integer; 

Xmax, Ymax, Xmin, Ymin, Xmid, Ymid, Xdiff, Ydiff : real; 
begin 

Xmin := XMn; 

Ymin := YMn{Round(YMn)} ; 

Xmax := XMx{Round(XMx)}; 

Ymax := YMx{Round(YMx)} ; 

DefineWorld(l, Xmin, Ymin, Xmax, Ymax); 

SelectWorld(l); 
end; { FindWorld } 

procedure DrawAxis{(Footer1 , Footer2 : WrkString; Arrows : boolean)}; 
var 

LineStyleLoc, XkO, YkO, Xkl, Ykl, Xk2, Yk2, 

MaxExponentX, MaxExponentY, TickPoint : integer; 

TickSmall, TickLarge, Max, Min, Tick, Offset, Diff : real; 

function Log(X : real) : real; 

{ Base 10 logarithm of X. } 
begin 

Log := Ln(X) / Ln(IO.O); 
end; { Log } 

function ALog(X : real) : real; 
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{ Ten raised to the X power. } 
begin 

ALog := Exp(X * Ln(10.0)); 
end; { Alog } 

function Frac(R : real) : real; 

{ Return the fractional part of the real number R. } 
begin 

Frac := R - Int(R); 
end; { Frac } 

procedure Ticks(NTicks : integer; Max, Min : real; 
var TickSmall, TickLarge : real); 

{ NTicks : The approximate number of tick marks in the interval. } 
{ Max, Min : World coordinates of the axis extremes. } 

{ TickSmall, TickLarge :Tick mark intervals, in world coordinates .} 
var 

TickLog : array[1..4] of real; 

I, J, ChA : integer; 

Delta, XTicks, LogTicks, Mant, MinDiff, Diff : real; 
begin 

TickLog[1] ;= 0.0; 

TickLog[2] := Log(2.0); 

TickLog [3] := Log (5.0); 

TickLog [4] := 1.0; 

XTicks := NTicks; 

Delta := Max - Min; 

XTicks := Delta / XTicks; 

LogTicks ;= Log (XTicks); 

ChA := trunc( LogTicks); 
if LogTicks < 0.0 then 
ChA := ChA - 1 ; 

MinDiff := 1 .0; 

Mant := LogTicks - ChA; { Fractional part of logarithm } 

for I := 1 to 4 do 

begin 

Diff := Abs(Mant - TickLog(l]); 
if (Diff < MinDiff) then 
begin 

MinDiff ;= Diff; 

J := I; 
end; 
end; 

LogTicks := ChA + TickLog[J]; { Logarithm of tick mark } 

TickLarge := A LOg( LogTicks); { The tick mark } 

{ Find the small tick marks, that are two tick scales smaller } 
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if J > 2 then 
begin 
J := J - 2; 

LogTicks := ChA + TickLog[J]; 
end 
else 
begin 
J := J + 1 : 

LogTicks := ChA + TickLog[J] - 1 ; 
end; 

TickSmall := ALog(LogTicks); 
end; { Ticks } 

function StringNumber(X1 : real; MaxExponent : integer) : WrkString; 
begin 

StringNumber := RealToStr(X1 * Exp(-MaxExponent * Ln(IO.O))); 
end; { StringNumber } 

function GetExponent(X1 : real) : integer; 
begin 

GetExponent := 0; 
if X1 <> 0.0 then 
if ABS(X1) >= 1.0 then 

GetExponent := trunc(Ln(ABS(X1)) / Ln(10.0)) 

else 

GetExponent := -trunc(ABS(Ln(ABS(X1))) / Ln(IO.O) + 1.0); 
end; { GetExponent } 

procedure DrawNum(X1, Y1, MaxExponent : integer; Number : real): 
var 

StrNumber : WrkString; 
begin 

TextSize(9): 

StrNumber := StringNumber(Number, MaxExponent): 

Y1 ;= Y1 - 3; 

MoveTo(X1, Y1): 

DrawString(StrNumber): 

TextSize(12); 
end; { DrawNum } 

procedure DrawExponent(X1 , Y1, MaxExponent : integer); 
var 

NumStr : WrkString; 
begin 

MoveTo(X1, Y1): 

TextSize{9): 

DrawChar('x'): 
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DrawCharC '); 

DrawChar(T): 

DrawChar(’O’); 

XI := WhereX + 1; 

Y1 := WhereY - 3; 

MoveTo(X1, Y1): 

NumStr := IntToStr(MaxExponent); 

TextSize(7): 

DrawString{NumStr): 

TextSize(12); 
end; { DrawExponent } 

begin { DrawAxis } 

LineStyleLoc := LinestyleGIb; 

SetLineStyle(O); { Black } 

XkO := XIRefGlb + X1 Offset; 

YkO := Y2RefGlb - Y20ffset; 

Xkl := XkO; 

Yk1 := YIRefGIb + Y1 Offset; 

Xk2 := X2RefGlb - X20ffset; 

Yk2 := YkO; 

MoveTo(XK0, YKO); { Draw the Y axis with optional Arrows } 
LineTo(XK1, YK1); 
if Arrows then 
begin 

MoveTo(XkO, Yk1); 

LineTo(XkO - 4, Yk1 + 4); 

MoveTo(XkO, Yk1); 

LineTo(XkO + 4, Yk1 + 4); 

DP(XkO, Yk1 - 1); 
end; 

MoveTo(XkO, YkO); { Draw the X axis with optional Arrows } 
LineTo(Xk2 + 1, Yk2); 
if Arrows then 
begin 

MoveTo(Xk2, Yk2); 

LineTo(Xk2 - 4, Yk2 - 4); 

MoveTo(Xk2, Yk2); 

LineTo(Xk2 - 4, Yk2 + 4); 
end; 

if Footer! <> " then { Draw the 1 st footer below the X axis } 
begin 

MoveTo(XkO, YkO + 45); 
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TextSize(9); 

DrawString(Footerl): 

end; 

if Footer2 <> ” then { Draw the 2nd footer below the X axis } 
begin 

MoveTo(XkO, YkO + 65); 

TextSize(9); 

DrawString(Footer2) ; 
end; 

if (ABS(YkO - Yk1) >= 35) and (ABS(Xk2 - Xk1) >= 150) then 
begin 

if ABS(Y2WldGlb) > ABS(YIWIdGlb) then 
MaxExponentY := GetExponent(Y2WldGlb) 
else 

MaxExponentY := GetExponent(YIWIdGlb); 

if MaxExponentY <> 0 then { Draw the power of ten on top of Y axis } 
DrawExponent(Xk1 - 30, Ykl + 2, MaxExponentY); 

TickPoint := YkO; 
if Y1 WIdGIb > Y2WldGlb then 
begin 

Max ;= YIVVIdGib; 

Min := Y2WidGlb: 
end 
else 
begin 

Max ;= Y2WldGlb; 

Min := YIWIdGIb; 
end; 

{ Using the Max and Min values, this procedure call calculates } 

{ large and small Tick Marks in world coordinates. } 

Ticks(5, Max, Min, TickSmall, TickLarge); 

Offset ;= Min / TickLarge; 

Offset := Offset - Frac(Offset); 

Tick := Offset * TickLarge; 
if Tick < Min then 
Tick := Tick + TickLarge; 

{ Tick is the world coordinate at which the tick mark is to be drawn } 
Diff := Max - Min; 

{ Plot large tick marks and Numeric labels } 
while Tick <= Max do 
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begin 

TickPoint := YkO - Trunc((YkO - Ykl) * (Tick - Min) / Diff); 
MoveTo(XkO, TickPoint); 

LineTo(XkO - 4, TickPoint); 

DrawNum(X1 Offset - 30, TickPoint + 7, MaxExponentY, Tick); 
Tick := Tick + TickLarge; 
end; 

{ The same repeated for the small tick marks, } 

{ only without axis numbering. } 

Offset := Min / TickSmall; 

Offset ;= Offset - Frac(Offset); 

Tick := Offset * TickSmall; 
if Tick < Min then 
Tick := Tick + TickSmall; 
while (Tick + 0.01) < Max do 
begin 

TickPoint := YkO - Trunc((YkO - Ykl) * (Tick - Min) / Diff); 
MoveTo(XkO, TickPoint); 

LineTo(XkO - 2, TickPoint); 

Tick := Tick + TickSmall; 
end; 

if ABS(X2WldGlb) > ABS(X1 WIdGIb) then 
MaxExponentX := GetExponent(X2WldGlb) 
else 

MaxExponentX := GetExponent(X1 WIdGIb); 
if MaxExponentX <> 0 then { Draw power of ten label on X axis } 
DrawExponent(Xk2 - 25, YkO + 28, MaxExponentX); 

{ This is the same as for the Y axis, but the window } 

{ and world are appropriate for the X axis. } 

TickPoint := XkO; 
if XI WIdGIb > X2WldGlb then 
begin 

Max := XI WIdGIb; 

Min := X2WldGlb; 
end 
else 
begin 

Max := X2WldGlb; 

Min ;= XI WIdGIb; 
end; 

Ticks(5, Max, Min, TickSmall, TickLarge); 

Offset := Min / TickLarge; 

Offset := Offset - Frac(Offset); 

Tick := Offset * TickLarge; 
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if Tick < Min then 
Tick := Tick + TickLarge; 

Diff := Max - Min; 
while Tick <= Max do 
begin 

TickPoint := XkO + Trunc((Xk2 - XkO) * (Tick - Min) / Diff); 
MoveTo(TickPoint, YkO); 

LineTo(TickPoint, YkO + 4); 

DrawNum(TickPoint - 14, YkO + 20, MaxExponentX, Tick); 
Tick ;= Tick + TickLarge; 
end; 

Offset := Min / TickSmall; 

Offset := Offset - Frac(Offset); 

Tick := Offset * TickSmall; 
if Tick < Min then 
Tick := Tick + TickSmall; 
while (Tick + 0.01) < Max do 
begin 

TickPoint := XkO + Trunc((Xk2 - XkO) * (Tick - Min) / Diff); 
MoveTo(TickPoint, YkO); 

LineTo(TickPoint, YkO + 2); 

Tick := Tick + TickSmall; 
end; 
end; 

SetLineStyle(LineStyleLoc); 

AxisGlb := TRUE; 
end; { Draw Axis } 

procedure ResetAxis; 
begin 

AxisGlb := true; 
end; { ResetAxis } 



procedure DrawPolygon{(A : PlotArray; First, NPoints, Line, Scale, 
Lines ; integer; Crosshairs : boolean)}; 



var 

i, XI, X2, Y1, Y2, XOffset, YOffset, 

XIRefLoc, YIRefLoc, X2RefLoc, Y2RefLoc, 

Delta Y, XOsI, XOs2, YOsI , YOs2 : integer; 

AutoClip, DirectModeLoc, PlotLine, PlotSymbol, Flipped : boolean; 
XILoc, YILoc, X2Loc, Y2Loc : integer; 

Temp : real; 

LineStyleLoc2 : integer; 

DrawPt : Boolean; 



procedure DrawPointClipped(X, Y : integer); 
begin 
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if (X1 > XIRefGIb) and (X2 < X2RefGlb) then 
if (Y1 > Y1 RefGIb) and (Y2 < Y2RefGlb) then 
DP(X, Y): 

end; { DrawPointClipped } 

procedure Drawltem(X, Y : integer); 
var 

LineStyleLoc : integer; 
begin 

LineStyleLoc := LineStyleGIb; 

SetLineStyle(O); { Black } 
case Line of 

2 : DrawCrossDiag(X, Y, Scale); 

3, 4 : DrawSquareC(X - Scale, Y + Scale, X + Scale, Y - Scale, (Line = 4)) 

5 : DrawDiamond{X, Y, Scale + 1); 

6 : DrawWye(X, Y, Scale +1); 

1 : DrawCross(X, Y, Scale); 

8 : DrawCircleDirect(X, Y, Scale + 1); 

9 : begin 

PlotLine := false; 
if AutoClip then 
DrawPointClipped(X, Y) 
else 

DP(X, Y); 
end; 

7 : DrawPoint(X, Y){, Scale)}; 

end; 

SetLineStyle(LineStyleLoc); 
end; { Drawitem } 

begin { DrawPolygon } 
if AxisGIb then 
Flipped := FALSE 
else 
begin 

Flipped := TRUE; 

Temp := World[WorldNdxGlb].Y1 ; 

World[WorldNdxGlb].Y1 := World[WorldNdxGlb].Y2; 
World(WorldNdxGlb].Y2 := Temp; 

SelectWorld(WorldNdxGlb) ; 

SelectWind(WindowNdxGlb, TRUE); 
end; 

if abs(NPoints - First) >= 2 then 
begin 

AutoClip := (NPoints < 0); 

NPoints := abs(NPoints); 

XOs1 := 1 ; 



149 



X0s2 ;= 1 ; 

YOsI :=6; 

YOs2 := 6; 
if AxisGIb then 
begin 

XOs1 :=X1 Offset: 

XOs2 := X20ffset; 

YOs1 :=Y1 Offset; 

YOs2 := Y20ffset: 

if {(X2RefG!b - XOs2 - XIRefGIb + XOs1) > (XOs1 + XOs2)) and 
((Y2RefGib - YOs2 - YIRefGIb + YOsI) > (YOs1 + YOs2)) then 
begin 

X1 RefLoc := X1 RefGIb; 

X1 := X1 RefGIb + XOsI ; 

Y1RefLoc:= YIRefGIb; 

Y1 := Y1 RefGIb + YOsI ; 

X2RefLoc := X2RefGlb; 

X2 := X2RefGlb - XOs2; 

Y2RefLoc := Y2RefGlb; 

Y2 := Y2RefGlb - YOs2; 

SetWindow(X1, Y1, X2, Y2); 

AxisGIb := TRUE; 
end; 
end; 

PiotLine := {Line >= 0); 

PlotSymbo! := (Line <> 0); 

Line := abs(Line); 

Scale := abs(Scale); 
if Lines < 0 then 

DeltaY := Trunc(1.0 / (abs(YIWIdGlb) + abs(Y2WldGlb)) * 

abs(YIWIdGlb) * abs(Y2RefGlb - YIRefGIb)) + 1 
else 

DeltaY := 0; 
if (NPoints < 2) then 

Error(’DrawPolygon #1’) 
else 
begin 

if CrossHairs then 
begin 

LineStyleLoc2 ;= LineStyleGIb; 

SetLineStyle(3): { Light Gray } 

MoveTo(X1 RefGIb, YIRefGIb + Y2RefGlb - WindowY(O.O)); 
LineTo(X2RefGlb, YIRefGIb + Y2RefGlb - WindowY(O.O)); 
MoveTo(WindowX(0.0), YIRefGIb); 

LineTo(WindowX(0.0), Y2RefGlb); 

SetLineStyle(LineStyleLoc2): 

end; 
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X1 := WindowX(A[First, 1]); 

Y1 := Y1 RefGIb + Y2RefGlb - WindowY(A[First, 2]) ; 
Drawltem(X1, Y1); 
if Abs(Lines) = 1 then 
if AutoClip then 

DrawLineClipped(X1, Y2RefGlb - DeltaY, XI, Y1) 
else 
begin 

MoveTo(X1, Y2RefGlb - DeltaY); 

LineTo{X1, Y1); 
end; 

DrawPt := True; 
for l:= First + 1 to N Points do 
begin 

X2 := WindowX(A[l, 1]); 

Y2 := Y2RefGlb + Y1 RefGIb - WindowY(A[l, 2]); 
Drawltem(X2, Y2); 
if StepA then begin 
if not AMarkStatus then 
begin 

if (A[l, 2] > 0) and DrawPt then {Clip(Xl,Y1,X2.Y2)} 
begin 

if ARJustification then 
begin 

MoveTo(X2+15, Y2+5); 

TextSize(9): 

DrawString(Ds): 

end 

else 

begin 

MoveTo(X2-65, Y2); 

TextSize(9); 

DrawString(Ds); 

end; 

DrawPt := False; 
end; 
end 
else 
begin 

if (I > (NPoints - InitDegree)) and (A[l, 2] > 0) and DrawPt 
begin 

if ARJustification then 
begin 

MoveTo(X2+15, Y2+5); 

TextSize(9); 

DrawString(Ds); 

end 



then 
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else 

begin 

MoveTo(X2-65, Y2); 

TextSize(9); 

DrawString(Ds); 

end; 

DrawPt := False; 
end; 
end; 

end;{StepA} 
if not StepA then begin 
if not BMarkStatus then 
begin 

if (All, 2] > 0) and DrawPt then {Clip(X1,Y1,X2,Y2)} 
begin 

if BRJustification then 
begin 

MoveTo(X2+15, Y2+5); 

TextSize(9); 

DrawString(Ds); 

end 

else 

begin 

MoveTo(X2-65, Y2); 

TextSize(9); 

DrawString(Ds); 

end; 

DrawPt := False; 
end; 
end 
else 
begin 

if (I > (NPoints - InitDegree)) and (A[l, 2] > 0) and DrawPt 
begin 

if BRJustification then 
begin 

MoveTo(X2+15, Y2+5); 

TextSize(9); 

DrawString(Ds); 

end 

else 

begin 

MoveTo(X2-65, Y2); 

TextSize(9); 

DrawString(Ds); 

end; 

DrawPt := False; 



then 
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end; 

end; 

end;{not StepA} 

if Abs( Lines) = 1 then 
if AutoClip then 

DrawUneClipped(X2, Y2RefGlb - DeltaY, X2, Y2) 
else 
begin 

MoveTo(X2, Y2RefGlb - DeltaY); 

UneTo(X2, Y2); 
end; 

if PlotLine then 
if AutoClip then 

DrawLineClipped(X1 , Y1 , X2, Y2) 
else 
begin 

MoveTo(X1, Y1); 

LineTo(X2, Y2); 
end; 

XI := X2; 

Y1 := Y2; 
end; 
end: 

if AxisGIb then 
begin 

SetWindow(X1 RefLoc, YIRefLoc, X2RefLoc, Y2RefLoc); 
AxisGIb ;= false; 
end; 
end 
else 

Error('DrawPolygon # 1'); 
if Flipped then 
begin 

Temp := World(WorldNdxGlb].Y1 ; 

World[WorldNdxGlb].Y1 := World(WorldNdxGlb].Y2; 
World[WorldNdxGlb].Y2 := Temp; 

SelectWorld( WorldNdxGlb) ; 

SelectWind(WindowNdxGlb, TRUE); 
end; 

end; { DrawPolygon } 

procedure Hatch{(X_1, Y_1, X_2, Y_2, Delta ; real)}; 
var 

X1, Y1, X2, Y2 : integer; 

procedure HatchDirect(X1 , Y1, X2, Y2, Delta : integer); 
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var 

I, Yst, Yen, Count : integer; 

XtRefLoc, X2RefLoc, YlRefLoc, Y2RefLoc : integer; 
ClippingLoc : boolean; 

XI D, Y1 D, X2D, Y2D : integer; *' 

begin { Hatch Direct } 
if Delta <> 0 then 
begin 

HatchGIb := true; 

ClippingLoc := ClippingGIb; 

ClippingGIb := true; 

XtRefLoc :=X1RefGlb; 

XI RefGIb := XI ; 

X2RefLoc := X2RefGlb; 

X2RefGlb := X2; 

YlRefLoc := Y1 RefGIb; 

Y1 RefGIb := Y1 ; 

Y2RefLoc := Y2RefGlb; 

Y2RefGlb := Y2; 

Yst := Y1 + Delta; 

Yen := Y1 - X2 + XI + Delta; 

if Delta < 0 then 

begin 

Delta := -Delta; 

! := Yst; 

Yst := Yen; 

Yen ;= I; 
end; 

Count := (Y2 - Y1 -t- X2 - XI -i- X2 - XI) div Delta; 

for I ;= 1 to Count-1 do 

begin 

XI D := XI; 

Y1 D := Yst; 

X2D := X2; 

Y2D := Yen; 

if Clip(X1D, Y1D, X2D, Y2D) then 
begin 

MoveTo(X1D, Y1D); 

LineTo(X2D, Y2D); 
end; 

Yst := Yst + Delta; 

Yen := Yen + Delta; 
end; 

ClippingGIb := ClippingLoc; 

HatchGIb := false; 

XI RefGIb := XtRefLoc; 
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X2RefGlb := X2RefLoc; 

YIRefGIb := YIRefLoc; 

Y2RefGlb := Y2RefLoc; 
end; 

end; { Hatch Direct } 
begin { Hatch } 

HatchDirect(trunc(X_1), trunc(Y_1), trunc(X_2), trunc(Y_2), trunc(Delta)) 
end; { Hatch } 

procedure DrawHistogram{(A :PlotArray; NPoints : integer; 

Hatching : boolean; HatchStyle : integer)}; 



var 

X1, X2, Y2, NPixels, Delta, NDiff, YRef, LineStyleLoc, I ; integer; 
Fract, S, Y ; real; 

DirectModeLoc, Negative : boolean; 

XIL 0 C, YIL 0 C, X2Loc, Y2Loc : integer; 

XtRefLoc, YIRefLoc, X2RefLoc, Y2RefLoc, YAxis : integer; 

Temp : real; 

begin { DrawHistogram } 
if ABS(NPoints) >= 2 then 
begin 

LineStyleLoc := LinestyleGIb; 

SetLineStyle(O); { Black } 
if AxisGIb then 
begin 

XtRefLoc := Window[WindowNdxGlb].X1 ; 

YIRefLoc := Window[WindowNdxGlb].Y1 ; 

X2RefLoc := Window[WindowNdxGlb].X2; 

Y2RefLoc := Window[WindowNdxGlb].Y2; 
SetWindow(XlRefGlb + XI Offset, YIRefGIb + Y1 Offset, 
X2RefGlb - X20ffset, Y2RefGlb - Y20ffset); 

AxisGIb := TRUE; 
end; 

Negative := NPoints < 0; 

NPoints := ABS(NPoints); 

NPixels := X2RefGlb - XtRefGIb; 

Delta := NPixels div NPoints; 

NDiff := NPixels - Delta * NPoints; 

Fract := NDiff / NPoints; 

S := -Fract; 

XI := XtRefGIb; 

Temp := Y2RefGlb + Yt RefGIb - AyGIb; 
if Temp > Maxint then 
Temp := Maxint 
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else 

if Temp < -32767 then 
Temp := -32767; 

YRef := trunc(Temp); 
if Negative then 

DrawStraight(X1, X2RefGlb, YRef); 

YAxis := Y1 RefGIb; 
if BYGIb > 0 then 
YAxis := Y2RefGlb; 
for f := 1 to NPoints do 
begin 

X2 := X1 + Delta: 

Y := A(l, 2]; 

if not Negative then 

Y := ABS(Y): 

Temp := AyGIb -i- ByGIb * Y; 
if Temp > Maxint then 
Temp ;= Maxint 
else 

if Temp < -32767 then 
Temp := -32767; 

Y2 ;= Y2RefGlb -i- Y1 RefGIb - trunc(Temp); 
if not Negative then 
begin 

MoveTo(X1 , YAxis): 

LineTo(X1, Y2); 

MoveTo(X1, Y2); 

LineTo(X2, Y2); 

MoveTo{X2, Y2): 

LineTo(X2, YAxis); 
if Hatching then 
if Odd(l) then 

Hatch(X1, Y2, X2, YAxis, HatchStyle) 
else 

Hatch(X1. Y2, X2, YAxis, -HatchStyle): 
end 
else 
begin 

MoveTo(X1, YReO: 

LineTo(X1, Y2); 

MoveTo(X1, Y2): 

LineTo(X2, Y2); 

MoveTo(X2, Y2); 

LineTo(X2, YReO: 
if Hatching then 
if YRef - Y2 < 0 then 
if Odd(l) then 
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Hatch(X1, YRef, X2. Y2, HatchStyle) 
else 

Hatch{X1, YRef, X2. Y2, -HatchStyle) 
else 

if Odd(l) then 

Hatch(X1 , Y2, X2,YRef, HatchStyle) 
else 

Hatch(X1, Y2. X2, YRef, -HatchStyle): 

end; 

X1 := X2; 
end; 

if AxisGIb then 
begin 

SetWindow(X1 RefLoc, Y1 RefLoc, X2RefLoc, Y2RefLoc); 
AxisGIb := FALSE; 
end; 

SetLineStyle(LineStyleLoc): 

end 

else 

Error{’DrawHistogram'): 
end; { DrawHistogram } 

begin 

end. { TurboGraph } 
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* This is the resource file that defines the menus and icons for MacRootLocus. 
MacRootLocus.Rsrc 
TYPEDUCX3 
,256 (36) 

CE Parameter Dialog 
70 100 300 412 
Visible NoGoAway 
1 
0 

256 

TYPE DITL 
,256 (36) 

11 

Btnitem Enabled 
185 240 210 300 
Cancel 

Btnitem Enabled 
150 240 175 300 
CK 

EditTextItem Enabled 
60 260 75 280 

EditTextItem Enabled 
95 165 110 205 
1 

EditTextItem Enabled 
120 165135 205 
100 

EditTextItem Enabled 
145 165160 205 
IE-6 

StatText Disabled 
20 15 40 250 

Characteristic Equation Parameter 

StatText Disabled 
60 40 75 240 
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Degree of the polynomial 

StatText Disabled 
95 40 110 150 
InitGuess 

StatText Disabled 
120 40135 150 
Maxiter 

StatText Disabled 
145 40160150 
Tolerance 



TYPEDLjOG 
,257 ( 36 ) 

One Parameter Root Locus Plot Data 
50 86 315 426 
Visible NoGoAway 
1 
0 

257 

TYPE DITL 
,257 ( 36 ) 

21 

Btnitem Enabled 
35 270 55 330 
Cancel 

Btnitem Enabled 
10 270 30 330 
PLOT 

Radioitem Enabled 
95 40 110 265 
Linear Point Interval 

Radioitem Enabled 
115 40 130 265 
Logarithmic Point Interval 

Radioitem Enabled 
145 25 160 150 
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Auto Scale Axis 



Radbltem Enabled 
145 165 160 320 
Manual Scale Axis 

EditTextIlem Enabled 

65 110 80160 

0.1 

EditTextIlem Enabled 
65 265 80 315 
10000 

EditTextItem Enabled 
170 85 185 125 
-10 

EditTextItem Enabled 
170 225 185 265 
5 

EditTextItem Enabled 
195 85 210 125 
-10 

EditTextItem Enabled 
195 225 210 265 
10 

EditTextItem Enabled 
230 135 245 165 
50 

StatText Disabled 
20 10 40 255 

One Parameter Root Locus Plot Data 

StatText Disabled 
65 25 80 100 
AMin Gain 

StatText Disabled 
65 180 80 260 
AMax Gain 

StatText Disabled 
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170 25 185 80 
XMin 

StatText Disabled 
170 165 185 220 
XMax 

StatText Disabled 
195 25 210 80 
YMin 

StatText Disabled 
195165 210 220 
YMax 

StatText Disabled 
230 25 245125 
Points To Plot 



TYPEDLDG 
.258 ( 36 ) 

Two Parameter Root Locus Plot Data 
30 86 330 426 
Visible NoGoAway 
1 
0 

258 

TYPE DITL 
,258 ( 36 ) 

35 

Btnitem Enabled 
35 270 55 330 
Cancel 

Btnitem Enabled 
10 270 30 330 
Plot 

Radioitem Enabled 
120 40135 265 
Linear Point Interval 

Radioitem Enabled 
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140 40 155 265 
Logarithmic Point Interval 

Radioltem Enabled 
250 130 265 180 
Start 

Radioltem Enabled 
250 185 265 230 
End 

Radioltem Enabled 
250 235 265 288 
Right 

Radioltem Enabled 
250 293 265 338 
Left 

Radioltem Enabled 
270 130 285 180 
Start 

Radioltem Enabled 
270 185 285 230 
End 

Radioltem Enabled 
270 235 285 288 
Right 

Radioltem Enabled 
270 293 285 338 
Left 

EditTextItem Enabled 

65 110 80 160 

0.1 

EditTextItem Enabled 
65 265 80 315 
10000 

EditTextItem Enabled 

90 110 105 160 

0.1 
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EditTextItem Enabled 
90 265 105 315 
10000 

EditTextItem Enabled 

170135185165 

5 

EditTextItem Enabled 
170 285185 315 
50 

EditTextItem Enabled 
195 85 210 125 
-10 

EditTextItem Enabled 
195 225 210 265 
5 

EditTextItem Enabled 
220 85 235125 
-10 

EditTextItem Enabled 
220 225 235 265 
10 

StatText Disabled 
20 10 40 255 

Two Parameter Root Locus Plot Data 

StatText Disabled 
65 25 80 105 
AMin Gain 

StatText Disabled 
65 180 80 260 
AMax Gain 

StatText Disabled 
90 25 105 105 
BMin Gain 

StatText Disabled 
90 180 105 260 
BMax Gain 
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StatText Disabled 
195 25 210 80 
X Min 

StatText Disabled 
195 165 210 220 
XMax 

StatText Disabled 
220 25 235 80 
YMin 

StatText Disabled 
220 165 235 220 
YMax 

StatText Disabled 
170 25 185 125 
How Many Loci 

StatText Disabled 
170 185 185 275 
Points To Plot 

StatText Disabled 
250 25 265 125 
AMark Point 

StatText Disabled 
270 25 285 125 
BMark Point 

TYPEDLOG 

.301 ( 36 ) 

Characteristic Equation Coefficient Data 
105 56 260 456 
Visible NoGoAway 
1 
0 

301 

TYPE DITL 
,301 
7 
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Btnitem Enabled 
15 320 35 380 
CK 

Btnitem Enabled 
50 320 70 380 
Cancel 

EditTextItem Enabled 
110 145 130 255 

EditTextItem Enabled 
110 270130 380 

StateTextItem Enabled 
20 15 40 305 

Characteristic Equation Coefficient Data 

StateTextItem Enabled 
85 148 105 208 
S“1 

StateTextItem Enabled 
85 273 105 333 
S**0 



TYPEDLDG 
,302 (36) 

Characteristic Equation Coefficient Data 
105 56 260 456 
Visible NoGoAway 
1 
0 

302 

TYPE DITL 
,302 
9 

Btnitem Enabled 
15 320 35 380 
CK 

Btnitem Enabled 
50 320 70 380 



Cancel 



EditTextItem Enabled 
110 20 130 130 

EditTextItem Enabled 
110145130 255 

EditTextItem Enabled 
110 270 130 380 

StaleTextItem Enabled 
2015 40 305 

Characteristic Equation Coefficient Data 

StateTextItem Enabled 

85 23 105 83 

S**2 

StateTextItem Enabled 
85 148 105 208 
S**1 

StateTextItem Enabled 
85 273 105 333 
S“0 

TYPEDUCX3 
,303 (36) 

Characteristic Equation Coefficient Data 
75 56 285 456 
Visible NoGoAway 
1 
0 

303 

TYPE DITL 
,303 
11 

Btnitem Enabled 
15 320 35 380 
CK 

Btnitem Enabled 
50 320 70 380 



Cancel 



EditTextItem Enabled 
110 270 130 380 

EditTextItem Enabled 
165 20 185130 

EditTextItem Enabled 
165 145185 255 

EditTextItem Enabled 
165 270 185 380 

StateTextItem Enabled 
20 15 40 305 

Characteristic Equation Coefficient Data 

StateTextItem Enabled 
85 273 105 333 
S**3 

StateTextItem Enabled 
140 23 160 83 
S**2 

StateTextItem Enabled 
140 148 160 208 
S**1 

StateTextItem Enabled 
140 273 160 333 
S**0 

TYPEDLOG 
,304 (36) 

Characteristic Equation Coefficient Data 
75 56 285 456 
Visible NoGoAway 
1 
0 

304 

TYPE DITL 
,304 
13 
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Btnitem Enabled 
15 320 35 380 
CK 



Btnitem Enabled 
50 320 70 380 
Cancel 

EditTextItem Enabled 
110 145 130 255 

EditTextItem Enabled 
110 270 130 380 

EditTextItem Enabled 
165 20 185 130 

EditTextItem Enabled 
165 145 185 255 

EditTextItem Enabled 
165 270185 380 

StaleTexiltem Enabled 
20 15 40 305 

Characteristic Equation Coefficient Data 

StateTextItem Enabled 
85 148 105 208 
S**4 

StateTextItem Enabled 
85 273 105 333 
S**3 

StateTextItem Enabled 
140 23 16*) 83 
S**2 

StateTextItem Enabled 
140 148 160 208 
S**1 

StateTextItem Enabled 
140 273160 333 
S**0 
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PiPEDLOG 



,305 ( 36 ) 

Characteristic Equation Coefficient Data 
75 56 285 456 
Visible NoGoAway 
1 
0 

305 

TYPE DITL 
,305 
15 

Btnitem Enabled 
15 320 35 380 
CK 

Btnitem Enabled 
50 320 70 380 
Cancel 

EditTextItem Enabled 
110 20 130 130 

EditTextItem Enabled 
110 145 130 255 

EditTextItem Enabled 
110 270130 380 

EditTextItem Enabled 
165 20 185 130 

EditTextItem Enabled 
165145185 255 

EditTextItem Enabled 
165 270 185 380 

StateTextItem Enabled 
20 15 40 305 

Characteristic Equation Coefficient Data 

StateTextItem Enabled 
85 23 105 83 
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S‘*5 



StateTextItem Enabled 
85 148 105 208 
S*M 

StateTextItem Enabled 
85 273105 333 
S**3 

StateTextItem Enabled 
140 23160 83 
S**2 

StateTextItem Enabled 
140148 160 208 
S“1 

StateTextItem Enabled 
140 273 160 333 
S**0 

TYPEDLOG 
,306 (36) 

Characteristic Equation Coefficient Data 
45 56 310 456 
Visible NoGoAway 
1 
0 

306 

TYPE DITL 
,306 
17 

Btnitem Enabled 
15 320 35 380 
CK 

Btnitem Enabled 
50 320 70 380 
Cancel 

EditTextItem Enabled 
110 270 130 380 
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EditTextlter: Enabled 
165 20 185 130 

EditTextItem Enabled 
165 145185 255 

EditTextItem Enabled 
165 270185 380 

EditTextItem Enabled 
220 20 240 130 

EditTextItem Enabled 
220 145 240 255 

EditTextItem Enabled 
220 270 240 380 

StateTextItem Enabled 
2015 40 305 

Characteristic Equation Coefficient Data 

StateTextItem Enabled 
85 273 105 333 
S*‘6 



StateTextItem Enabled 
140 23 160 83 
S**5 



StateTextItem Enabled 
140 148 160 208 
S*M 



StateTextItem Enabled 
140 273 160 333 
S“3 



StateTextItem Enabled 
195 23 215 83 
S**2 



StateTextItem Enabled 
195148 215 208 
S**1 



StateTextItem Enabled 



195 273 215 333 
S“0 



TYPEDLOG 
.307 (36) 

Characteristic Equation Coefficient Data 
45 56 310 456 
Visible NoGoAway 
1 
0 

307 

TYPE DITL 
,307 
19 

Btnitem Enabled 
15 320 35 380 
CK 

Btnitem Enabled 
50 320 70 380 
Cancel 

EditTextItem Enabled 
110 145 130 255 

EditTextItem Enabled 
110 270 130 380 

EditTextItem Enabled 
165 20 185 130 

EditTextItem Enabled 
165 145 185 255 

EditTextItem Enabled 
165 270 185 380 

EditTextItem Enabled 
220 20 240 130 

EditTextItem Enabled 
220 145 240 255 

EditTextItem Enabled 
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220 270 240 380 



StateTextItem Enabled 
20 15 40 305 

Characteristic Equation Coefficient Data 

StateTextItem Enabled 
85 148 105 208 
S**7 

StateTextItem Enabled 
85 273 105 333 
S*‘6 

StateTextItem Enabled 
140 23 160 83 
S**5 

StateTextItem Enabled 
140 148 160 208 
S*M 

StateTextItem Enabled 
140 273 160 333 
S**3 

StateTextItem Enabled 
195 23 215 83 
S**2 

StateTextItem Enabled 
195 148 215 208 
S**1 



StateTextItem Enabled 
195 273 215 333 
S**0 

TYPEDLDG 
.308 (36) 

Characteristic Equation Coefficient Data 
45 56 310 456 
Visible NoGoAway 
1 
0 

308 
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TYPE DITL 
,308 

21 

Btnitem Enabled 
15 320 35 380 
CK 

Btnitem Enabled 
50 320 70 380 
Cancel 

EditTextItem Enabled 
11020 130 130 

EditTextItem Enabled 
110 145 130 255 

EditTextItem Enabled 
110 270 130 380 

EditTextItem Enabled 
165 20 185 130 

EditTextItem Enabled 
165 145 185 255 

EditTextItem Enabled 
165 270 185 380 

EditTextItem Enabled 
220 20 240 130 

EditTextItem Enabled 
220 145 240 255 

EditTextItem Enabled 
220 270 240 380 

StateTextItem Enabled 
20 15 40 305 

Characteristic Equation Coefficient Data 

StateTextItem Enabled 
85 23 105 83 
S“8 
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StateTextItem Enabled 
85 148 105 208 

S**7 



StateTextItem Enabled 
85 273105 333 
S**6 

StateTextItem Enabled 
140 23 160 83 
S**5 

StateTextItem Enabled 
140 148 160 208 
S**4 



StateTextItem Enabled 
140 273 160 333 
S**3 



StateTextItem Enabled 
195 23 215 83 
S**2 

StateTextItem Enabled 
195 148 215 208 
S**1 

StateTextItem Enabled 
195 273 215 333 
S**0 

TYPEDLDG 
,309 (36) 

Characteristic Equation Coefficient Data 
40 56 315 456 
Visible NoGoAway 
1 
0 

309 

TYPE DITL 
,309 
23 
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Btnitem Enabled 
15 320 35 380 
CK 

Btnitem Enabled 
50 320 70 380 
Cancel 

EditTextItem Enabled 
75 20 95130 

EditTextItem Enabled 
130 20 150 130 

EditTextItem Enabled 
130 145 150 255 

EditTextItem Enabled 
130 270 150 380 

EditTextItem Enabled 
185 20 205 130 

EditTextItem Enabled 
185 145 205 255 

EditTextItem Enabled 
185 270 205 380 

EditTextItem Enabled 
240 20 260 130 

EditTextItem Enabled 
240 145 260 255 

EditTextItem Enabled 
240 270 260 380 

StateTextItem Enabled 
20 15 40 305 

Characteristic Equation Coefficient Data 

StateTextItem Enabled 

50 23 70 83 

S**9 

StateTextItem Enabled 
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105 23 125 63 
S**8 

StateTextItem Enabled 
105 148 125 208 
S**7 

StateTextItem Enabled 
105 273 125 333 
S“6 

StateTextItem Enabled 
160 23 180 83 
S**5 

StateTextItem Enabled 
160 148 180 208 
S*M 

StateTextItem Enabled 
160 273 180 333 
S**3 

StateTextItem Enabled 
215 23 235 83 
S**2 



StateTextItem Enabled 
215 148 235 208 
S**1 



StateTextItem Enabled 
215 273 235 333 
S‘*0 

TYPEDUOG 

,310(36) 

Characteristic Equation Coefficient Data 
40 56 315 456 
Visible NoGoAway 
1 
0 

310 

TYPE DITL 
,310 
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25 



Btnitem Enabled 
15 320 35 380 

o< 

Btnitem Enabled 
50 320 70 380 
Cancel 

EditTextItem Enabled 
75 20 95 130 

EditTextItem Enabled 
75145 95 255 

EditTextItem Enabled 
130 20 150 130 

EditTextItem Enabled 
130 145 150 255 

EditTextItem Enabled 
130 270 150 380 

EditTextItem Enabled 
185 20 205130 

EditTextItem Enabled 
185 145 205 255 

EditTextItem Enabled 
185 270 205 380 

EditTextItem Enabled 
240 20 260 130 

EditTextItem Enabled 
240 145 260 255 

EditTextItem Enabled 
240 270 260 380 

StateTextItem Enabled 
20 15 40 305 

Characteristic Equation Coefficient Data 
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StateTextItem Enabled 

50 23 70 83 

S**10 

StateTextItem Enabled 
50 148 70 208 
S**9 

StateTextItem Enabled 
105 23 125 83 
S**8 

StateTextItem Enabled 
105 148 125 208 
S**7 



StateTextItem Enabled 
105 273 125 333 
S**6 



StateTextItem Enabled 
160 23 180 83 
S“5 



StateTextItem Enabled 
160 148 180 208 
S**4 



StateTextItem Enabled 
160 273 180 333 
S**3 

StateTextItem Enabled 
215 23 235 83 
S**2 

StateTextItem Enabled 
215 148 235 208 
S**1 



StateTextItem Enabled 
215 273 235 333 
S**0 



TYPE MENU 

,1000 
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\14 

About MacRootLocus . 
(- 

,1001 

File 

EQ Parameter/E 
Get Coeff^G 
Print Screen /S 
Print Winck)w /R 
QuiVQ 

,1002 

Edit 

Undo /U 
(* 

Cut /X 
Copy /C 
Paste /V 
Clear 

,1003 

Plot 

One Parameter/0 
Two Parameter/T 

,1004 

Help 

EQ Parameter/A 
Get Coeff/F 
One Parameter/N 
Two Parameter/W 
Print Out/P 



LIST OF REFERENCES 



1. Apple Computer, Inz., Macintosh System Software User’s Guide, pp. 
123-154,1988. 

2. Graham, N., Introduction to Computer Science: A Structured Approach, 2nd 
ed., pp. 239 — 249, West Publishing Company, 1982 

3. Ralston, A., A First Course in Numerical Analysis, pp. 368 — 371, McGraw- 
Hill Book Co., 1978. 

4. Borland International, Inc., Turbo Pascal Tutor (a Self-study Guide to 
Turbo Pascal on the Macintosh), 1987. 

5. Borland International, Inc., Turbo Pascal for the Mac (User’s Guide and 
Reference Manual), 1986. 



181 



INITIAL DISTRIBUTION LIST 



NO. Copies 



1. Defence Technical Information Center 2 

Cameron Station 

Arlington, Virgini 22304—6145 

2. Library, Code 0142 2 

Naval Postgraduate School 

Monterey, CA 93943-5002 

3. Chairman, Code 62 1 

Department of Electrical and Computer Engineering 

Naval Postgraduate School 
Monterey CA 93943-5000 

4. Professor George J. Thaler, Code 62Tr 5 

Naval Postgraduate School 

Monterey CA 93943—5000 

5. Professor Harold A. Titus, Code 67Ts 1 

Naval Postgraduate School 

Monterey CA 93943—5000 

6. Professor Sellung Kwak, Code 52Kw 1 

Naval Postgraduate School 

Monterey CA 93943—5000 

7. Professor James W. Bamham 1 

Assoc. Head. Dept of Mech Engineering & Mechanics 

College of Engineering Drexel Univ. 

Philadelphia PA, 19104 

8. Major Ko, Sung Hoon 5 

286—38 SooYu 4 Dong DoBong Gu Seoul 132 

Seoul Korea 

9. Park, Seung Chin 1 

3401 N. Columbus #16j 

Tucson AZ. 85712 

10. Kim Yoo Chang 1 

460 W. Forest ave #903 

Detroit Michigan 48201 



182 



1 



11. Lcdr Kenneth Mac Donald 
1524 Dolphin Court 
Orange Park Florida 32673 



12. Lcdr Roy Wood 1 

1522 Oak Knoll Road 

Virginia Beach, Virginia 23414 

13. Lcdr Hwang, Jung Sub 1 

SMC1209 Naval Postgraduate School 

Monterey, CA 93943 

14. Seo Yong Seok 1 

SMC1448 Naval Postgraduate School 

Monterey, CA 93943 

15. Hong-on Kim 1 

SMC2665 Naval Postgraduate School 

Monterey, CA 93943 

1C. Kang, Mung Hung 1 

SMC1375 Naval Postgraduate School 
Monterey, CA 93943 

17. Kwon, Hui Man 1 

SMC1375 Naval Postgraduate School 

Monterey, CA 93943 

18. Yang, Young Hyee 1 

SMC2440 Naval Postgraduate School 

Monterey, CA 93943 



183 



ip! 



Thesis 

K716248 Ko 

c.l MacRootLocus, a CAD 

design tool for feed- 
back control systems. 



Thesis 
K716248 Ko 

c.’ MacRootLocus, a CAD 

design tool for feed- 
back control systems. 



